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Heat Transfer From Transverse and Yawed 


LIONEL V. BALDWIN 
VIRGIL A. SANDBORN 
JAMES C. LAURENCE 


Aeronautical Research Scientists, 
Lewis Research Center, National 
Aeronautics and Space Administration, 
Cleveland, Ohio 


Cylinders in Continuum, Slip, and Free 
Molecule Air Flows 


A general Nusselt number correlation is presented for transverse cylinders in subsonic 
and supersonic air flows where dissociation is negligible. New and existing data in 
the following experimental range have been correlated: 


Mach M = 0.001 to 6.0, 


Reynolds Nx. = 0.02 to 300,000, Knudsen Nx, = 4 X 10~* to 37. 

In subsonic flow, heat transfer from a cylinder yawed at an angle to the air velocity is 
predicted by the transverse cylinder Nusselt number correlation when the normal velocity 
component is used as the characteristic velocity. Finally, recovery temperature data from 
cylinders are divided into three regimes by a Knudsen number criterion. 


Introduction 


OR MANY YEARS PAPERs reporting experimental data 
of heat transfer from cylinders have appeared in the literature. 
This can probably be attributed not only to the relative ease with 
which electrically heated wires lend themselves to such an experi- 
ment, but also to the practical importance of cylinders as sensitive 
elements for instruments. For example, thermocouple response 
characteristics are important for many control applications. 
Similarly, hot-wire anemometers are valuable tools for studying 
both the mean and turbulent properties of air flows in aerodynamic 
research; and at very low density, hot wires find application as 
pressure gages of the Pirani-type. Despite these familiar applica- 
tions, engineers do not generally realize that the accumulated 
data on heat transfer from cylinders are the most complete in- 
formation available today for the prediction of heat transfer in 
rarefied air. The utility of individual contributions is greatly en- 
hanced when these are viewed as part of a summary. But un- 
fortunately, there has been little effort to correlate recent experi- 
mental data over a wide range of flow conditions. This paper 
presents the results of such an attempt with special emphasis on a 
general correlation of heat-transfer data from transverse and 


1This paper has been written as a compilation of two papers, 
“Heat Transfer From Cylinders in Subsonic Slip Flow,’ presented 
at the 1958 Annual Meeting (Paper No. 58—A-202), and “Heat 
Transfer From Yawed Cylinders in Rarefied Air Flows,"’ presented 
at the Heat Transfer Conference, 1959 (Paper No. 59—HT-5). Both 
papers were contributed by the Heat Transfer Division of Tue 
AMERICAN Society oF MecHanicat ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, April 
28, 1959. 
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yawed cylinders in rarefied air flows. The discussion of this 
Nusselt number correlation stresses the primary dependence of 
the data on the flow characteristics, which vary from free mole- 
cule flow to continuum, and from speeds near those of natural 
convection to supersonic. Secondary effects, such as the de- 
pendence of the heat-transfer coefficient on cylinder and air tem- 
perature, are reviewed only briefly. Since data on the recovery 
temperature of cylinders give an insight to the boundary cri- 
terion for free molecule flow, these results are also presented 
herein. 


Dimensionless Groups of Correlation 


A significant portion of the experimental results which are to 
be presented were obtained in the last two years at the Lewis 
Research Center. The data from these experiments are corre- 
lated by dimensionless groups such as Nusselt, Reynolds, Mach, 
or Knudsen number. Like previous investigators, we were con- 
fronted with the problem of evaluating the air density, viscosity, 
and thermal conductivity for these dimensionless groups. The 
temperature at which to evaluate these air properties is the one 
that gives the “best correlation.’’ Certainly the forms that 
were chosen in the past could be justified, but we hope that the 
success of the proposed Nusselt-number Reynolds-number corre- 
lation presents a good case for the method used in this paper. 

The convective heat-loss rate defines the heat-transfer coefficient 
h: 


The heat-transfer coefficient is defined in terms of the recovery 
temperature 7. The heat-transfer data are presented in terms 


49.02 +/7, speed of sound, ft/sec 
exponential, e*, dimensionless 


dimensionless function of s; f(s) 


>) | 


D = cylinder diameter, ft 


8? 3? 3? 
dimensionless function of s; g(s) = 39 | exp — = x + 2 + 
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of a Nusselt number based on air thermal conductivity (Hilsen- 
rath [1]}*) at total air temperature 


Nu = — (2) 


The Reynolds iwmber for the cylinder is defined by the free- 
stream density and velocity, but with air viscosity at total air 
temperature: 


pU..D 


The Knudsen number is based on the mean free path at free- 
stream conditions and cylinder diameter: 


Nre = 


(3) 


Ns = D (4) 
Kennard, 1938 [2], Equation (126b), gives the mean free path in 
terms of the viscosity of a rarefied gas composed of hard elastic 
spheres having a Maxwellian velocity distribution; 4 = 0.499%pX. 
As is usually done, this relation together with experimental 
viscosity is taken as the definition of the mean free path. The 
Knudsen number is therefore related to the ratio of the Mach 
number to Reynolds number, since there is a simple relation 
between the acoustic speed and the mean molecular speed 6 (e.g., 
Stalder, Goodwin, and Creager [3)). . 


ry M 


(5) 
2 


= 
The ratio of viscosity at static temperature to that at total tem- 
perature is required by the definition of the Reynolds number 
which is used herein. All viscosities have been taken from Hilsen- 
rath [1]. 
The Mach number is based on free-stream velocity and the 
acoustic velocity at static air temperature: 


ua (6) 
a, 

All the data presented have been converted to this basis. Where 
there has been a choice of measurements of different wire tem- 
peratures, the data nearest (7',, — 7',.) = 500 F were used for 
reasons which will be given later. 


? Numbers in brackets designate References at end of paper. 


Nusselt-Number Correlation for Transverse Cylinders 

The results of individual investigators are grouped according to 
experimental range for ease in the graphical presentation of re- 
sults. In this presentation we are primarily concerned with 
demonstrating the variation of heat-transfer characteristics over 
a very wide range of flows. 

A summary of some subsonic continuum data is given in Fig. 1. 
Hilpert’s data [4] are the backbone of current engineering correla- 
tions (McAdams [5]). The results of Schmidt and Wenner [6] 
for large cylinders confirm the results of Hilpert at the high 
Reynolds numbers. The data of Collis and Williams [7] and 
Cole and Roshko [8] extend the familiar correlation for con- 
tinuum flow by two decades to Nre = 0.01. These data were 
obtained with small diameter wires having length to diameter 
ratios greater than 2000. In Fig. 1, the results of Cole and 
Roshko fall slightly below the Collis and Williams data. A 
possible explanation is that Cole and Roshko reported Nusselt 
numbers which were extrapolated to AT’ = 0. The Collis and 
Williams data which are used here were obtained at AT’ ~ 420 F, 
In this experimental range, the heat-transfer coefficient increases 
slightly as cylinder temperature is increased. 

The effects of natural convection were observed by Collis and 
Williams as the Reynolds number decreases (e.g., in Fig. 1, 
0.00021-in. cylinder Nx. < 0.2 or the 0.0000354-in. cylinder Nr. 
< 0.03). Cole and Roshko report that decreasing the aspect 


ratio of the cylinder increases the Nusselt number (after allow- . 


ance for conduction losses at the ends). This observation is 
limited to very low Nre results at nearly zero Mach number be- 
cause no effect of aspect ratio is discernible in the rarefied flow 
results to be presented in the following discussion. 

As the air velocity increases, the Mach number can be ex- 
pected to be a flow parameter along with Nr. Scadron and 
Warshawsky [9] varied cylinder diameter and air density for sub- 
sonic Mach numbers from 0.10 to 0.90. Their experimental pro- 
cedure differed from that used by most of the investigators men- 
tioned here in that they measured transient thermal response 
rather than steady-state he atloss. Therefore the Nusselt num- 
bers that ‘they calculate from the time-constant of step-current 
response are for very low cylinder temperatures. In terms of our 
equation (3), their results were correlated within 7.4 per cent 
average deviation of a single observation by: 


Nyu = 0.431 (i 0.20 M2) (7) 


g(s sin ®) = dimensionless function of s sin ®; substitute 
s sin ® for s in g(s) formula 
9. = unit conversion factor, lbyft/lbr sec* 
h = average convective heat-transfer coefficient for 
whole cylinder, Btu/(sec) (sq ft)(deg R) 
Io(z),Ji(z) = modified Bessel function of first kind; (0) zero and 
(1) first order, respectively, dimensionless 
J = conversion factor, ft-lbr/Btu 
k, = air thermal conductivity at total air temperature, 
Btu/(sec) (ft) (deg R) 


l = cylinder length, ft 
£ = turbulent integral scale, ft 
M = free-stream Mach number, dimensionless 
Nx» = free-stream Knudsen number, dimensionless 
Nyu = Nusselt number, dimensionless 
Nre = Reynolds number, dimensionless 
p = static pressure, lbr/sq ft 
q- = over-all convective heat-transfer rate, Btu/sec 
R = gas constant for air, ft-lbr/Iby-deg R 
s = \ = M, molecular speed ratio, dimensionless 
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T = static or free-stream air temperature, deg R 
T, = recovery or equilibrium cylinder temperature, deg 
R 
T..2 = theoretical recovery temperature of an infinitely 
long cylinder in a free-molecule flow of air, deg R 


T, = total air temperature, deg R 
T,, = cylinder temperature, deg R 
AT = T,, — T,, temperature difference, deg R 
U.. = free-stream air speed, ft/sec 
§ = ye", mean molecular speed, ft/sec 
Vm = V 2g.RT, most probable molecular speed, ft/sec 
a@ = accommodation coefficient, dimensionless 
Y = ratio of specific heats, 1.40 for air, dimensionless 
» = T,/T,, recovery temperature ratio, deg R 
\ = mean free path of air at free-stream conditions, ft 
M, = air viscosity at total air temperature, lby/ft-sec 
mw = 3.14159 
® = angle the free-stream velocity vector subtends with 
the cylinder axis, radians 
p = air density at free-stream conditions, lby/cu ft 
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Fig. 1 Summary of subsonic continuum heat transfer from tr ylind 


The 0.05 Mach line calculated from equation (7) is shown dashed 
in Fig. 1 in the range of their work 250 < Nre < 30,000. At least 
some of the difference in the dashed line and the summary curve 
is traceable to the secondary effect of cylinder temperature on 
the heat-transfer coefficient. The important point of the data of 
Scadron and Warshawsky is the relatively small effect of Mach 
number on heat transfer at constant Reynolds number that was 
observed [equation (7)] as M increased to 0.90. This Mach 
number effect is practically nil beiow M = 0.30 and it decreases 
the Nusselt number by about 15 per cent at M = 0.90. Naturally, 
more data are required before firm conclusions can be drawn, 
but until more data become available, an allowance of the order 
of V 1 + 0.20 M? might be used to account for the effect of com- 
pressibility in subsonic continuum flows. 

The remainder of this section is devoted to heat-transfer data 
from cylinders in rarefied air flows. The data shown in Fig. 2 
extend from the slip region to free molecule flow and are subsonic 
throughout. The Mach number is a parameter shown by solid 
line curves. The early results of Laurence and Landes [10], 
seemed anomalous when the Mach number was shown to be a 
parameter in the range of 0.10 < M < 0.30 at Nre approaching 
100. Subsequent experiments showed that this ‘“Mach-number 
effect’’ was real and the anomalous features faded when it was 
noted that an equally effective correlating parameter would be 
the Knudsen number. Curves of constant Knudsen number 
are shown in Fig. 2 as dashed lines. The heat-loss characteristics 
shown in this figure differ markedly from the familiar NR. corre- 
lation of continuum flow. The significant features of Fig. 2 are 
discussed as part of the summary correlation at the end of this 
section. For the moment, recall that the three flow parameters 
(Nre, M, Nxa) are related by equation (5); so only two moduli 
can be regarded as independent. Therefore either Nyu = f(Nre, 
M) or Nxu = f(Nre, Nxa) completely specifies the correlation. 
Where it is convenient, both of these correlation forms are given 
for ready reference (as in Fig. 2). 

The data of Baldwin [11, 12] and Cybulski and Baldwin [13] 
were obtained with tungsten wires about 0.00020 in. in diameter. 
The aspect ratio was about 350 in the earlier work and 1/D = 
2500 for the more recent measurements. The good agreement at 
Nxa < 0.01 in Fig. 2 is evidence that the aspect ratio phenomenon 
discussed by Cole and Roshko, does not extend into the rarefied 
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flow region. Although the data of Laurence and Landes and 
Baldwin are not in good agreement, both sets of data exhibit the 
same trends as the flow parameters are varied. 

In Fig. 2, the data of Cybulski and Baldwin at Nx. = 4.70 
borders on the free molecule region, which will shortly be shown 
to exist for Nx. > 5.0 from recovery temperature data. The 
following equation gives the heat-transfer coefficient h of a trans- 
verse cylinder of infinite length in a rarefied, diatomic gas flow 
(Stalder, Goodwin, and Creager [3]) 


om pv,,ag (8) 


(8) 
Re-expressing equation (8) in terms of the Nusselt and Knudsen 
numbers: 
Ru, ) g(8) 
= —— —— — (allM 9 


At a total air temperature of 85 F, the dimensionless constant in 
parentheses is 0.0322. g/(s) is a function of the Mach number (or 
molecular speed ratio); which is defined in the Nomenclature. 
At Mach numbers greater than about 2.4, 


g(s) = 6 
2 
Combining equations (10) and (5) with equation (9) and evaluat- 


ing the gas properties for air streams at a total temperature of 85 
F, we get 


(M > 2.4) (10) 


= 0.193 aNre 


If the accommodation coefficient a of the metal-air system is 
known, equations (9) and (11) predict the heat transfer. 

The data of Cybulski and Baldwin were found to be in best 
agreement with equation (9) when a = 0.57. This agreement is 
shown in Fig. 2 for Nxa > 5.0. The theory correctly predicts 
both the first power dependence at constant M and also the large 
separation of the constant M parametric curves. 

The data of Spangenberg [14] are shown in Fig. 3. Spangen- 
berg varied cylinder diameter from 0.00005 to 0.00152 in. and 
worked with three wire materials: tungsten, platinum, and a 
platinum 10 per cent rhodium alloy. He reported more data than 
have been plotted here, but the portion that is not shown confirmed 


(M > 2.4) (11) 
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and overlapped the data in Fig. 3. The subsonic data in Fig. 3 
confirm all the trends noted in Fig. 2. The agreement between 
Figs. 2 and 3 is good for M > 0.60, but the low M results of 
Spangenberg are somewhat higher in Nu than the Baldwin data. 
Nevertheless, the similarity between Figs. 2 and 3 is considered 
a positive result, although the discrepancies at M < 0.6 are larger 
than would be expected for these experiments. The supersonic 
data of Fig. 3 group in a narrow band, so for 1.05 < M < 1.90, 
the Nusselt number is adequately given by specifying the 
Reynolds number alone. 

Some additional experimental data are given in Fig. 4. The 
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Fig. 3 Heat transfer from transverse cylinders in slip flow 


results of Lowell [15] and Winovich and Stine [16] agree well in 
the transonic range. The supersonic data of Winovich and Stine 
groups in a single Nx. band as did that of Spangenberg in Fig. 3. 
However, Lowell reports that the supersonic Nw. data have a 
systematic variation with Mach numbers greater than 1.0. 
Note that the supersonic Nx, data of Lowell fall over the transonic 
results. However, the trend is relatively small and if it truly 
exists, there is no evidence of it in the data of either Spangenberg 
or Winovich and Stine. 

A summary plot of all available supersonic data, including some 
not previously published, is given in Fig. 5. The data group in 
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a narrow band with a scatter of the order of +15 per cent. The 
very low Reynolds number data of Stalder. Goodwin, and Creager 
[3], Wong [17], and Weltmann and Kuhns [18] were obtained in 
nitrogen rather than air, however, this is thought to have no ef- 
fect on the heat transfer. While free-moiecular theory predicts 
that the data should be represented by a single curve for M > 2.4 
[see equation (11)], it was not foreseen that the data in the 
transition to free molecule and slip regions should also be a single 
curve. The data are limited to Nae < 1000, however, there 
is no reason to suspect that the correlation will differ greatly at 
higher Nr as long as dissociation does not occur. 

The two curves given as Figs. 1 and 5 are the boundaries of the 
convective heat-transfer correlation for transverse cylinders which 
is given as Fig. 6. The continuum curve at very low Mach num- 
ber (Fig. 1) is well established over a large range of Reynolds 
numbers. The mid-range between the incompressible, continuum 
flow data and the supersonic band of data is covered by the sub- 
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Fig. 4 Heat transfer from transverse cylinders in transonic slip flow 


sonic rarefied flows such as shown in Figs. 2 and 3. The data of 
Baldwin and Cybulski were chosen in preparing the suggested 
correlation curves of Fig. 6, because they cover the entire subsonic 
range of rarefied flows. Future work may alter these sub- 
sonic curves in the rarefied regime, but there is no doubt that the 
present curves correctly predict the trends and are sufficiently 
accurate for most engineering work. As discussed in the section 

n “Secondary Effects,” future fundamental data on the ac- 
commodation coefficient should also help clarify the subsonic 
portion of Fig. 6. 

The subsonic flow data in the transition region from slip to 
free-molecular flow show the gradual change in slope of the con- 
stant Mach number curves necessary to approach the first power 
Reynolds-Nusselt number relation required for free-molecular 
flow. As the Reynolds number decreases, the heat transfer de- 
pends less on the air veiocity (as compared to the continuum 
curve), and instead the Knudsen number becomes the governing 
parameter in subsonic flows. 

The supersonic curve appears to fair in very closely with the 
subsonic continuum data at high Reynolds numbers. In the ab- 
sence of supersonic data at Reynolds numbers greater than 10° it 
would appear necessary to assume that the subsonic curve is ade- 
quate to predict heat transfer in supersonic flows. 

Recently, heat-transfer measurements for spheres have been 
correlated by Schaaf and Chambri [22], in the same manner as 
used for Fig. 6. The results for spheres are quite similar to Fig. 6, 
although the data for spheres are quite meager by comparison. 


Secondary Effects on the Nusselt Number Correlation 


This brief review will not adequately describe all the observed 
effects, because in some cases the basic parameters have not been 
isolated. However, the state of the art is better today than even 
a few years ago, and future research should clarify some of the 
phenomena which occur. 
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It has long been recognized that the turbulence of the approach 
air stream can affect the observed heat-transfer characteristics. A 
recent review by Van der Hegge Zijnen [23] seems to clarify this 
problem in subsonic flows. He characterizes the turbulent field by 
the turbulent intensity and the ratio of the turbulent scale to 
cylinder diameter £/D. He shows that, at very small or very 
large £/D, the effect of turbulence on heat transfer is relatively 
small for all intensities, but for £/D ~ 1.5, an increase in turbu- 
lent intensity has a very large effect which is especially pro- 
nounced at large Nre. Concerning the data presented herein, 
it seems safe to neglect the effect of turbulence on the data from 
fine wires (say Nr. < 100). However, the turbulence must be 
regarded as an uncontrolled parameter for most of the results 
referred to here in the continuum region. 

Other secondary effects are usually associated with the wire 
temperature. For continuum flow, an increase in body tempera- 
ture 7’,, causes an increase in the heat-transfer coefficient h. This 
effect led to the use of “film temperatures’’ for evaluating air 
properties in engineering correlations (McAdams, also Douglas 
and Churchill [24]). In rarefied air flows, the observations are 
more complex (e.g., Baldwin). In subsonic flows at Nx. > 0.02, 
the heat-transfer coefficient h decreases with increasing cylinder 
temperature for 33 < AT < 200 F, but A assumed a constant 
value for 200 < AT < 500 F. The results were dependent on 
Nixa (or Nx.) and to a lesser extent on M. In rarefied flows at 
supersonic speed, the heat-transfer coefficient also decreases with 
rising 7',,, but this dependence on the flow parameters is different 
as shown by Winovich and Stine. From a practical viewpoint, 
this means that the “film temperature” type of correlation will not 
minimize the effect of body temperature from data in rarefied air 
flow, but rather will accentuate the systematic spread. 
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The dependence of the heat-transfer coefficient on body tem- 
perature reverses as Reynolds numbers corresponding to rarefied 
flow are approached. Very few fundamental data exist on the 
variation of the accommodation coefficient a with metal tem- 
perature, but a recent summary reported by Estermann [25] 
shows that @ exhibits a marked variation with metal temperature 
in the neighborhood of 100 to 600 F. This variation is one possible 
explanation of the reversal noted. As more data become available 
on @, it might be fruitful to reanalyze some of the data from fine 
wires at various 7. It is an experimental fact that data at very 
low AT’ scatter more than for higher operating temperatures. 
Furthermore, according to Estermann, there is less likelihood that 
a is strongly dependent on 7’, for T,, ~ 500 F. The accommoda- 
tion coefficient must be regarded as an uncontrolled parameter in 
the data reported herein and it is likely to have increasing impor- 
tance as Nx, rises above 0.01. 

Finally, the data presented were all obtained at essentially am- 
bient total temperature. There is some evidence that air tempera- 
ture has a nearly negligible effect on the Nusselt-number correla- 
tion in the form used here for temperatures as high as 3000 F 
(Baldwin, also Glawe and Johnson [26]). However, very high 
Mach number (high temperature) flows where air dissociation 


becomes important are obviously out of the range of the data in 
this review. 


Nusselt-Number Correlation for Yawed Cylinders 


Experimental data for heat transfer from yawed cylinders are 
scant compared with the transverse orientation. This discussion 
of heat transfer from yawed cylinders is limited to subsonic free- 
stream velocities in the continuum and slip flow regimes. A 
theoretical solution for the heat transfer from a yawed cylinder 
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in a free molecule flow is now available for all Mach numbers. 

For a transverse cylinder in subsonic continuum flow, Nxu was 
found to be a function of Nx. with perhaps a 15 per cent variation 
at constant Nr. due to a Mach number effect at 0.56 < M < 1.0 
(Fig. 1 and also Scadron and Warshawsky [9]). Reshotko and 
Beckwith [27] give a solution of the compressible laminar 
boundary layer equations at the stagnation-line of yawed cylin- 
ders. For Mach numbers from zero to about 0.5, they found that 
the effect of yaw can almost be eliminated by giving Nwu as a 
function of Nae sin ®, where ® is the angle the free-stream 
velocity vector subtends with the cylinder axis (® = 90 deg for 
transverse case). Thus for angles from 90 to about 20 deg, use 
of the summary correlation (Fig. 6) with Nx. sin ® should give a 
good engineering estimate for Nx. sin ® > 400 in subsonic flows. 

The anemometer data of Schubauer and Klebanoff [28] are 
sometimes cited as experimental evidence for the validity of Nr. 
sin ® as the correlating parameter for yawed cylinders. Actually, 
these authors’ contribution was to show that U sin ® correlated 
heat loss data obtained at various velocities and angles. That is, 
for their application, they did not need to vary D or p systemati- 
cally. Data were subsequently obtained with fine wires at vari- 
ous Nre by Spangenberg and by Sandborn and Laurence [29] but, 
as could be expected from the transverse cylinder discussion, an 
additional flow parameter is required along with Nre sin ® due 
to slip effects. The importance of Nxs for transverse orientation 
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in slip flow suggests that the heat-transfer results for yawed 
cylinders should correlate using Nx, and Npe sin ® as parameters. 
By equation (5) it can be seen that M sin ® is the related modu- 
lus. 

The data of Sandborn and Laurence have been reanalyzed by 
Baldwin, Sandborn, and Laurence [30] to test the use of Nre sin 
® and Nx, (or M sin ®) for correlation. The results are shown in 
Fig. 7. The points are keyed to values of the free stream Mach 
number and the angle of yaw. Some of the data at M = 0.70 
and 0.90 fail to align with the transverse cylinder data (® = 90 
deg). However, the general agreement of the majority of the 
yawed results with ® = 90 deg and the grouping of the yawed wire 
data for M sin ® = 0.10, shows that this form of correlation is 
an adequate assumption for most engineering purposes in the sub- 
sonic slip flow regime. 

Fig. 8 shows the subsonic data of Spangenberg. The data of 
Figs. 7 and 8 show that either M sin ® or Nx» is adequate to 
correlate the data for angles between 90 and 20 degrees. 

Although there appear to be no experimental data for yawed 
cylinders in very rarefied flow, the change in orientation only 
slightly modifies the free molecule flow analysis. The transverse 
cylinder problem is detailed by Stalder, Goodwin, and Creager 
{[31]; and Talbot [32] outlined the yawed cylinder solution. The 
similarity is clear when the result is written as: 


g(s sin ®) yp, 


Nyu = 0.0322 (all M) 


(12) 
or for M > 2.4, 
= 0.193 aNre sin ® 


(M > 2.4) (13) 
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The Mach number function g(s sin ®) is the same as in equation 
(9) except that it is a function of M sin ®, These forms include 
the transverse orientation as a special case, and equations (12) 
and (13) are applicable over the entire range of yaw angles, 90 to 
Odeg. For axial flow past the cylinder, the heat transfer becomes 
nearly independent of the free-stream Mach number. A solution 
of the form, 

Nwu = 0,303 


- (® = 0 deg) 
My 


(14) 
is equivalent to the solution for a flat plate at zero angle of attack 
as given by Oppenheim [33]. Note that since Nx, and Nx, each 
appear to the first power, the characteristic length may be can- 
celed from the equation. As will be discussed shortly, equations 
(12) to (14) can be expected to hold for Nx» > 5.0. 


Recovery Temperature of Transverse Cylinders 


The adiabatic or recovery temperature of transverse cylinders 
in continuum flows has been found to depend primarily on the 
Mach number of the air stream. Typical curves for subsonic 
flows are given in Fig. 9. The solid curve which reaches a 
minimum at M = 0.75 is from Simmons [34]; the shaded area 
indicates the spread of data about the mean. Simmons used long 
thermocouple wires stretched across an airjet and he varied wire 
diameter, air density, and velocity at 20,000 > Nr. > 1000. The 
abrupt minimum at M ~ 0.75 was observed only with cylinders 
having a length to diameter ratio greater than about 50. Note 
also the increase in the data spread near this hump. If a pair of 
support prongs were moved into position near the long thermo- 
couple wire (not in contact with the wire, but in front of it to form 
the usual H profile of probe supported couples) the hump was no 
longer observed. In the latter case, the recovery temperature 
was found to be more like the dashed curve of Fig. 9, which was 
taken from Scadron and Warshawsky. The data of Scadron and 
Warshawsky are from 30,000 > Nr. > 250, but the thermocouples 
were probe supported. It seems probable that the disappearance 
of the distinct minimum at M ~ 0.75 is associated with an inter- 
ference of the wire supports with the flow over the cylinder at 
aspect ratios less than 50. 
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RECOVERY TEMPERATURE RATIO, 


The preceding discussion is intended to show that, although 
the recovery temperature data are correlated by Mach number 
alone, there may be applications where the cylinder aspect ratio 
is a significant factor. However, the data given in Fig. 10 were 
all obtained with fine wires having length to diameter ratios 
greater than 100, and these curves have no sharp minimums. In 
Fig. 10, the Reynolds number range is lower and extends well 
into the slip regime (1000 > Na. > 0.10). From the general 
agreement between the data of Figs. 9 and 10, we conclude that 
within experimental scatter the recovery temperature data are 
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almost entirely free of Reynolds number dependence for at least 
30,000 > Nr. > 0.10. 

The recovery temperatures of continuum and slip flow are all 
less than the total air temperature. However, one of the interest- 
ing predictions of free-molecular-flow analysis is that the recovery 
temperature of a transverse cylinder will exceed the total tem- 
perature of the gas stream in which it is located. Thus the re- 
covery temperature of an infinitely long cylinder T,,,.. when placed 
in a diatomic gas flow is given by Stalder, Goodwin, and Creager 
as: 

T f(s) 
go — 


T, g(s) (15) 


That is, theory for very rarefied flow predicts that 7 will be a 
function only of Mach number. Equation (16) is plotted in Fig. 
11, where the experimental evidence supporting this prediction is 
also given. The Knudsen number of each data point is also 
shown. It is clear that there is a transition region between the 
continuum and slip flow behavior of the curve of Lowell in Fig. 11 
and the free molecule flow limitingcurve. In this transition region, 
7 depends on both M and Nx». Rough limits on this transition 
from slip to fully developed free molecule flow can be set at 5.0 > 
Nxa > 0.2. 

Fig. 12 is a crossplot of the data of Fig. 11 (on which the 
measured recovery temperature is compared to the prediction of 
free-molecule-flow analysis). Unlike inferences drawn from Nus- 
selt number data, the recovery temperature results enable us to 
test free-molecule-flow analysis independent of the accommoda- 
tion coefficient. So, the general grouping of both the subsonic and 
supersonic data at Nx. > 5.0 in Fig. 12 is excellent evidence for 
the existence of fully developed free-molecule flow. 


Summary and Conclusions 


A general Nusselt-number correlation is presented as Fig. 6 for 
transverse cylinders in subsonic and supersonic air flows where 
air dissociation is negligible. This empirical correlation covers 
flow conditions varying from continuum to the free molecule 
regime. It is based on data obtained in the following experi- 
mental range: Mach M = 0.001 to 6.0, Reynolds Nx, 0.02 to 
300,000, Knudsen Nx, 4 X 10~ to 37. The familiar Nusselt- 
number cortelation in subsonic continuum flow is a function only 
of Reynolds number. This curve has been extended by two 
decades at low Reynolds numbers to Nae = 0.01. The first 
deviation from the continuum correlation occurs at about Nre = 
400, and for Nr. < 400, an additional flow parameter along with 
Nnre is required to specify the Nusselt number. This additional 
parameter may be taken to be either the Knudsen number or the 
Mach number, because there is a relation between the three 
groups (Nre, Nxa,M). Asareminder of the physical significance 
of the additional parameter, it is preferable to use the Knudsen 
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number since the need for the additional modulus is due to the 
beginning of rarefied flow phenomena. However, for graphical 
presentation it is more convenient to specify the correlation by 
Nwu = f(N re, M) throughout the rarefied flow region. Thus, for 
all Nre < 400, the Nusselt number is given as a continuous func- 
tion of Reynolds number by a family of Mach number curves for 
0.05 <M < 1,0. These parametric curves have two outstanding 
characteristics: (a) The slope of the constant M curves on the 
plot of log Nwu versus log Ne gradually increases from about 0.5 
at Nx. = 0.001 to the 1.0 values predicted by free-molecule- 
flow theory at Nx. = 5.0; (6) the spread of these parametric 
curves is greatest in free-molecule flow (Nx, > 5.0) and the 
M curves gradually converge to a common point on the continuum 
Nr. curve at Np. > 400. The other important feature of the gen- 
eral Nx, correlation iu rarefied flows is the grouping of all super- 
sonic-flow data into a single band, so that Nyu is essentially a 
function of Nae only for all Nre below 1000. Furthermore, it 
appears safe to assume that the correlation given as Nyw = 
(Ne) for all supersonic M will persist into continuum flow above 
Nre = 1000, but lack of data prohibits firm conclusions. 

In subsonic flow, heat transfer from a cylinder yawed at an 
angle ® to the air velocity U.. is predicted by the transverse 
cylinder Nusselt number correlation when the normal velocity 
component U,. sin ® is used as the characteristic velocity. That 
is, both Reynolds and Mach numbers must be based on U,, sin ® 
even for slow subsonic flow, so that the yawed wire results are 
compared at constant Knudsen number. In continuum and slip 
flows, Nx. data at angles from 90 deg (transverse) to about 20 deg 
are correlated in this manner, while free-molecule-flow theory pre- 
dicts this procedure is valid for this regime over the entire range, 
90 to zero degrees. 

Secondary dependence of the Nusselt number on cylinder and 
air temperatures causes up to 20 per cent deviations from the 
general correlation. These effects are complicated, but appar- 
ently systematic, and a brief review is included in this paper. 

Finally, recovery temperature data from cylinders are divided 
into three parts by a Knudsen number criterion: Nx. < 0.2, the 


. 
ratio of recovery to total air temperature 7 = ( r ) is independ- 
t 


ent of Nx» and 7 decreases from 1.0 to 0.95 as Mach number in- 
creases from 0 to 4.0; 0.2 < Na < 5.0, 7 is a function of both 
Nx» and M in this transition region to fully developed free- 
molecule flow; Na > 5.0, 9 predicted by free-molecule-flow 
theory is observed, and 7 increases from 1.0 to 1.167 as M in- 
creases from 0 to 7.0 again independent of Ns. 
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Magnetohydrodynamic Effects Upon Heat 
Transfer for Laminar Flow 
Across a Flat Plate 


Forced-convection heat transfer for laminar flow of electrically conducting fluids across 
a flat plate is considered for a magnetic field of constant inductance acting normal to 
the free stream velocity and fixed relative to the plate. The boundary condition on the 
surface of the plate is taken to be either a constant temperature or constant heat flux, and 
solutions are presented for the following cases: (a) Fluids having a Prandtl number of 
unity for which both Joule heating and frictional heating are accounted for; (b) fluids 
having moderate and large Prandtl numbers for negligible Joule and frictional heating; 
and (c) fluids having low Prandtl numbers for negligible frictional heating. 


Introduction 


ONSIDERABLE INTEREST has recently developed in 
regard to forced-convection heat transfer involving electrically 
conducting fluids in the presence of magnetic fields. The primary 
motivation of this interest has been the prospect of reducing aero- 
dynamic heat transfer [1, 2, 3, 4, 5].!_ Rossow [1] has considered 
a fluid with a Prandtl number of unity flowing across a flat plate 
for which a magnetic field of constant inductance acts normal to 
the free stream velocity. The magnetic field was considered fixed 
relative to the plate or relative to the fluid, and solutions were ob- 
tained for the plate temperature equal to the free stream tem- 
perature and for an adiabatic plate. 

Although the fluids of most common concern have been con- 
ducting gases, it is not hard to envision practical heat-transfer 
problems involving other types of electrically conducting fluids, 
such as salt water or liquid metals. 


The purpose of the present investigation is to extend certain 
of Rossow’s results for laminar flow across a flat plate to include 
more general boundary conditions on the plate surface as well 
as more general types of fluids. The magnetic field is considered 
to be fixed relative to the plate and to act normal to the free stream 
velocity with constant magnetic inductance. It is further con- 
sidered that the plate surface has either a constant temperature 
or constant heat flux, and solutions are presented for the following 
cases: 


1 Fluids having a Prandtl number of unity for which both 
Joule heating and frictional heating are taken into account. 

2 Fluids having moderate and large Prandtl numbers (ap- 
proximately greater than one) for negligible Joule and frictional 
heating. 

3 Fluids having low Prandtl numbers for negligible frictional 
heating. 


The first case would be approximately representative of many 
gases, and is essentially the same problem considered by Rossow 
with the exception that the boundary conditions are more general. 
The second case would correspond to ordinary types of conducting 
liquids, such as salt water, for which both Joule heating and fric- 
tional heating are generally negligible. The third case applies 
to liquid metals such that frictional heating is negligible whereas 


1 Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division of THe American 
Society or MecHAnicaL EnGIngeers and presented at the ASME- 
AIChE Heat Transfer Conference, Storrs, Conn., August 9-12, 1959. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, May 4, 
1959. Paper No. 59—HT-14. 


menclature 


B. = magnetic induction vector (mag- 
nitude Bo) 

specific heat at constant pressure 

dimensionless stream functions 

total enthalpy, c,t + u?/2 

dimensionless enthalpy function, 
see Equations (14) and (25) 

thermal conductivity 

magnetic parameter, 7By"/pu. 

pressure 

Peclet number, 

Prandtl number, v/a 

local surface heat: flux 

Reynolds number, u,.z/v 

Laplace transform variable 

temperature 

dimensionless temperature dif- 
ference, see Equations (33) and 
(42) 
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PP 


velocity vector (magnitude U) 
velocity component in z-direction 
velocity component in y-direction 
distance co-ordinate measured 
along plate from leading edge 


distance co-ordinate measured 
normal to plate 


= (8.Pr/12)n* 

= thermal diffusivity 

= fo"(0) = 0.332 

= fi"(0) = —0.894 

= gamma function 

= temperature difference, t — t, 
Laplace transform of 0 
kinematic viscosity 


= SPPR a 


Subscripts 
w = conditions at the plate surface 
o = free stream conditions 
Superscripts 
’ = differentiation with respect to 7 
* = conditions for no magnetic field pres- 
ent 


density 

electrical conductivity 

time 

hydrodynamic potential function 

dissipation function 

unit hydrodynamic 
function, ®/u, 


potential 


WV = hydrodynamic stream function 
y = unit hydrodynamic stream func- 


tion, V/u, 
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Joule heating may be significant. For cases 1 and 2 only first- 
order heat-transfer results will be considered. 


Basic Equations 


The present analysis is based upon the following assumptions: 1.25 , 
1 The flow is incompressible, steady, and two-dimensional. Pa 
2 The properties of the fluid, including electrical conductivity, 


are constant. 
3 Both the induced magnetic field and the charge density are 

negligible. 
In accord with these assumptions the equations expressing the 


conservation of mass, momentum, and energy for laminar flow 
may be written, respectively, as [1] 


) and f, 

N 


7 
div =0 / 
0.50 
D 
DU wg +2 x Be) / 
Dr p p 
Dt 1 DP / 
= + — — + — & +— (OX (3) 0.25 
Dr pe, Dr Cy pc, 


where ®, is the dissipation function [6]. The last terms on the 
right side of Equations (2) and (3) represent, respectively, 


the induced magnetic force and Joule heat generation within 3 6 
the fluid due to the resulting electric current. ny [3 

‘From Rossow [1] the boundary layer form of Equation (2) for 
flow across a flat plate becomes Fig. 1 Zero-order and first-order stream function derivatives from 


references [1] and [6] 


4 
baile oy oy? p (4) Correspondingly, from continuity 


and from this equation Rossow has shown that the zero-order and 


first-order terms for the velocity components u and v are? oy oz p 
= + +...) (5) such that 
oBe? 
v= —y (9) 
+ maf,’ +...] 
— The stream and potential functions V and ® thus become 
V2 fo + 3mafi +...) (6) oBy? 
V= f (—vdz + udy) = ugy — zy 
where the primes denote differentiation with respect to 7. The mF 
function is the dimensionless Blasius stream function [6], oB? 
and f,() is given by (udz + = — “2p — 
1 3 
+ — + So"fi = fo’ (7) Further, defining unit stream and potential functions as = 
WV/u, and @ = P/u,, then 
with the boundary conditions y =y — may (10) 
fi=fi' =0; 7 =0 
fi =-1; = — y’) (11) 
$ 2 


Tabulated values of f,(y) and its first two derivatives are given 
in reference [1],* and the derivatives fy’ and f,’ are shown in Fig. 1. 
It is of further interest to determine the velocity components 


for the potential flow occurring outside the region of viscous effects Prandt] Number of Unity 
(momentum boundary layer). For this case Equation (4) be- 


These potential flow results will be employed when considering 
fluids with low Prandt! number. 


eames Taking into account both Joule heating and frictional heating 
the boundary-layer form of Equation (3) reduces to 
du _ 
ot ot du\? Bo? 
and ox oy oy \ oy 
Bo? 
u=U, - m4 z (8) However, it is more convenient in this case to employ the total 


— | energy equation rather than Equation (12). Upon multiplying 
8 Rossow has aleo considered second-order terms. Equation (4) by u and adding the result to Equation (12), the 
3 In reference [1 the first-order stream function is denoted by fo(n). | total energy equation for Pr = 1 becomes [1] 
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v 
or 


where A, is the total enthalpy given by 
hs = Cot + u?/2 


Constant Surface Temperature. For the case of a constant surface 
temperature on the plate, a dimensionless enthalpy function will 
be defined as 


and Equation (13) becomes 
oH OH 


= 15 
v (15) 


The zero-order and first-order terms of the enthalpy function may 
be expressed as 


H = Hy(n) + mzH,(n) +... (16) 


Upon substituting Equations (5), (6), and (16) into Equation (15), 
and collecting like powers of mz, the equations for Hy and H, be- 
come 


1 
H,” + fo'H,’ = 0 (17) 


1 3 
H,’ + — fo'H, 2 (18) 


with the boundary conditions 
H.=1, Hi =0; 7 =0 
H, = 0, H, = 0; n= @ 


The solution of Equation (17) has been shown by Pohlhausen [6] 
to be 


Hy = 1 — fi’ (19) 
and correspondingly Equation (18) becomes 


H," + SoH,’ — fo'H, Sifo” 


If a transformation is considered of the form 
H, = 2fo'fr’ (21) 


it may readily be shown that the function g() is analogous to the 
temperature function considered by Rossow [1] for the case of 
equal plate and free stream temperatures, and it follows that 
g'(0) = —0.352. 

From Equations (19) and (21) 


H,'(0) = —0.332 
H,'(0) = 0.352 
such that 


ot * 
( = -(0. {0.332 — 0.352(mz) + ...] 
Oy yo 2c, 


Further 
ot 
(=) 
oy y=0 
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and the expression for the heat flux from the surface of the plate 
becomes 


q = 0.332k (0. [1 — 1.06(mz) +...] (22) 


This may be written in terms of a heat-flux ratio as 


= 1 — 1.06(mz) +... (23) 
q | 


where q* is the heat flux with no magnetic field present given by 
k 

* = 0.332 {| — }) (6, — (Re)'” 24 

q 33 ( )( ye ) (24) 


which is the well-known result obtained by Pohlhausen [6}. 

One may note that the ratio g/g* is independent of the tem- 
perature difference 0,, and thus Rossow’s results for 0, = 0 are 
identical with Equation (23). This invariance with 6, is not ob- 
tained for Prandtl] numbers other than unity, except for cases in 
which Joule and frictional heating are negligible. 

Constant Surface Flux. For the case of a constant surface heat 
flux the dimensionless enthalpy function will be defined as 


Again letting 

H = H,(n) + mzHi(n) +... 
then from Equations (5), (6), and (13) 


Hy” + fell! — fr'H, = 0 (26) 


1 3 1 3 
H,” + 2 SoH,’ fi'Ho 2 (27) 


with the boundary conditions 
H,’ = 1, H,' = 0; 7 =0 
Hy, = 0, H, = 0; 7 = @ 


Numerical solution of Equations (26) and (27) has been carried 
out on a Datatron computer, giving the results 


H,(0) = —2.18 
H,(0) = —1.59 


H(n = 0) = —2.18 — 1.59(mz) +... 


such that the expression for the surface temperature difference 
becomes 


= 
2c 


+ 2.18 [1 + 0.729(mz) + ...] (28) 


This may be written in terms of a temperature difference ratio as 
0. Ue */2c, 
=1+ 0.7: 29 
— wa*/2c, + 0.729(mz) + (29) 


where 0, * is the surface temperature difference with no magnetic 
field present given by 


6,* = 2.18 (“) (Re) ~'/2 (30) 
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It may be noted from Equation (28) that the adiabatic wall 


To 1, = 0; 9= 0 
temperature is independent of the magnetic field, which is the 


same result found by Rossow [1]. T%=0; © 
From Equations (32), Equation (35) may be rewritten as 
Moderate and High Prandtl Numbers a 
Neglecting both Joule heating and frictional heating, the T.” + 7 BoPrn*T.’ = 0 
boundary layer form of Equation (3) becomes 


oz oy oy? Pr\‘/? P. 
1 (82) [ — POPE es) at (37) 
It will now be convenient to employ an approximation in regard 0 


to the velocity components given by Equations (5) and (6). For 

nonmagnetic flows (mz = 0) Fage and Falkner [7], as well as 

others, have considered the stream function derivative fo’ to be a 1 

linear function of 7. This approximation is reasonable for fluids aahict 6'/°T (4/3) 

with large Prandtl numbers, since for a large Prandtl number the 

thickness of the thermal boundary layer is much less than that Correspondingly, Equation (36) becomes 

of the momentum boundary layer, and thermal gradients are con- 1 

sequently restricted to the region for which fo’ is approximately 17,” + — @)Prn*7;' — BoPrnT: 

linear with 7. For low Prandtl numbers this would not be the 4 

case. However, from the heat-transfer results of Fage and Falk- 3 BoPr\'/ B.Pr 

ner for nonmagnetic flow across an isothermal flat plate, it is exp (- 

found that for a Prandtl number as low as 0.77 the linear ap- 

proximation is in error by only 2.3 per cent. and it may readily be verified that the solution of this equation 
In the present analysis it will be assumed that both fo’ and f,’ is 

are linear functions of 7, such that 


where 


fe = Bun? (5) new (- 
(32) From Equations (37) and (38) 
1 
2 
2 Bin T,'(0) = —¥; 


where = fo"(0) = 0.332 and = f,;"(0) = —0.894 from ref- 
erences [6] and [1], respectively. As may be seen from Fig. 1 T,'(0) = _3 vy (4) (2) : 
7 


the linear approximation is more restrictive in regard to f,' than 
fo’. However, it will be shown later that first-order heat-transfer 
results obtained by this means are approximately valid for Prandtl such that 
numbers as low as one.‘ ar B.Pr\'/* 
Constant Surface Temperature. For a constant temperature on the (27) = -79, ( ) 2 


surface of the plate, a dimensionless temperature difference will oy 
3 


2 
be defined as 


t—t. 0 
T= — = (33) 
ty — te w and the heat flux from the surface of the plate is given by 
and Equation (31) beeomes 3 
q = we 3 ime) +...] 
or or or (34) 2 vz 7 Bo 
oz oy dy? 


(39) 


Upon letting 


This may be written in terms of a heat-flux ratio as 
T = To(n) + mzT\(n) +... q 


= 1—1.15(mz) +... (40) 
the equations for the zero-order and first-order temperature func- ° 
tions, 79 and 7, are found from Equations (5), (6), and (34) to 3) . 
where g* = 0.338 (Pr)'/*(Re)'/2 (41) 
To” + 2 Prf.To’ = 0 (35) One may note that the ratio ¢/q* is independent of Prandtl num- 
_ ber for the present case of moderate and large Prandtl numbers. 
1 3 Constant Surface Flux. For the case of a constant surface heat 
Ty’ + > PrfTi’ — Prf’T: = — PrfiTo’ (36) flux the dimensionless temperature difference will be defined as 
(tte 1) 
with the boundary conditions T= bo — 9) 4] eo (42) 


q/k vz 


‘From reference [1] the linear approximation would not appear t . 
applicable to terms higher than the first-order. Again letting 
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T = To(m) + m2Ti(n) +... 
then from Equations (5), (6), and (31) 


1 1 
To” + PrfeTo’ — Prf'Ts = 0 (43) 


1 3 l 3 
Ti" + 2 PrfoT 1’ — Prf.’T, = Prf,'T> — PrfiTo’ (44) 


with the boundary conditions 
= 1, =0; =0 
T™=0, =0; n= 
From Equations (32), Equation (43) becomes 
1 1 
+ 4 BoPrn?T 2 BoPrnTo =0 
Letting z = (8)Pr/12)n’, this equation reduces to Kummer’s dif- 


ferential equation [8, 9], such that the solution may be written in 
terms of the Kummer function F(a, c, z) as 


5 4 
(62) F (3. 3° | (45) 


gis 
T(5/3) 


where "= Al 
Correspondingly, Equation (44) becomes 
1 3 
T," + BoPry?T,’ — BoPrnT; 
1 3 
= 2 - B:Prn?T 


and the solution of this equation may also be written in terms of 
Kummer functions as 


(46) 


(62)'/*F (3,4, 


(4/3) 


From Equations (45) and (46) 
/s 
0)=- 
T(0) v2 ( in) 


and it follows that the surface temperature difference is given by 
2 q ‘Ue. 
= % k vz 
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= 1 + 0.789(mz) +. (48) 


where 
2.16(% ) (49) 


Again it may be noted that for the present case magnetic effects 
are indeperdent of Prandtl number. 


Low Prandtl Numbers 
For negligible frictional heating Equation (3) may “ written 
approximately as 
ot ot ( ot ) o By? 


= + byt U (50) 
The method that will be employed to solve this equation is similar 
to that presented by Grosh and Cess [10] for nonmagnetic flows. 
For fluids with very low Prandtl numbers the thickness of the 
thermal boundary layer is much greater than that of the momen- 
tum boundary layer. This in turn indicates that viscous effects 
are limited to a small portion of the thermal boundary layer, and 
to a first approximation inviscid potential flow may be assumed 
such that the velocity components u and v occurring in Equation 
(50) are given by Equations (8) and (9). A discussion of the ap- 
plicability of this assumption for nonmagnetic flow across an iso- 
thermal flat plate is given in reference [10], for which good agree- 
ment was found to exist between inviscid flow and viscous flow 
solutions within the Prandtl number range of typical liquid 
metals. 

Following reference [10], Boussinesq’s transformation [11] will 
be employed. Thus Equation (50) may be written in terms of 
the independent variable ¢ and y as 


20 ( | 20 
* oy 

where @ and y are given by Equations (10) and (11). From the 
usual boundary layer assumptions this equation reduces to 

ag Cy 
for which the boundary condition with respect to @ may be taken 


as 9 = Ofor} = 0. Applying a Laplace transform with respect 
to ¢, then 


Cp 


Ue (51) 


and Equation (51) transforms to 


Uns 
dy? - - -(™=') (52) 
such that the solution may be carried out for specified boundary 
conditions with respect to y. 

Constant Surface Temperature. For a constant temperature on 


the surface of the plate, the boundary conditions with respect to y 
become 


y=0 
00 


and it follows that the transformed boundary conditions are 
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dy 
The solution of Equation (52) satisfying these conditions is 
9, mu,.* MU 
- -lepi -§ + 
8 a c,8 
Correspondingly 
€3) 
dy ymo a 


and taking the inverse it follows that 


| 
[1 -2( cms) (1 4 ma) 


To obtain an expression for the heat flux from the surface of 
the plate 


08 dy 00 
q= oy = —k(l — mz) (2). 


3 


Ue? 1 
[1 2(2 ) (mz) (: mz) | (53) 


Again expressing the results in terms of a heat-flux ratio 


q 1 — mz 


= =— 


\ 
1 = mr 
2 


q* = 0.564 (°) (Pe)'/* (55) 
z 


where 


The ratio g/q* is shown in Fig. 2. For most cases of practical 
interest the effect of Joule heating, as represented by the parame- 
ter (u..*/O,¢,), will be negligible. 

Constant Surface Flux. For the case of a constant surface heat 


flux 
( ) 
Oy you k 
and since 


then in terms of ¢ and W the temperature gradient at the surface 
becomes 
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Fig.2 Variation of the heat-flux ratio for fluids with low Prandtl numbers 


(3) 

oy k V/1 — 

Letting Q(s) represent the Laplace transform of this quantity, 
the boundary conditions for Equation (52) become 


—— 0; 


The solution of Equation (52) satisfying these conditions is 


2 
5- ep (-v + 
a cs 


from which it follows 


1 
G(s, 0) = Qe) + — (56) 


Recalling that 
q 1 
L£-(Q(s)} 
k — 


then with the aid of the convolution integral the inverse of Equa- 
tion (56) becomes 


Evaluating the integral in this equation, the surface tempera- 
ture difference is given by 


ku 1 — mz 


1 
+ —— (mz) mz) (57) 


Cp 


( 1 + Vmz(2 — =| 
In 
1 — mz 


+ (mz) (: mz) (58) 
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Fig. 3. Variation of the temperature difference ratio for fluids with low 
Prandtl numbers 


where 


0," 1.128 =) 


(59) 


The temperature difference ratio 0,/0,,* is shown in Fig. 3. In 
most cases the effect of Joule heating will again be negligible. 


Discussioa of Results 


Mathematically, the solutions for moderate and large Prandtl 
numbers represent asymptotic solutions for Pr — @, whereas 
the low Prandtl number results correspond to Pr —~ 0. It is of 
interest to compare these solutions as well as those for a Prandtl 
number of unity to determine the effect of Prandtl number upon 
convective heat transfer in a magnetic field. In order to put the 
results for all three cases on a comparable basis, both Joule and 
frictional heating will be neglected; i.e., the parameters (u?,/ 
O,c,) and (u,.*/O0,,*c,) will be set equal to zero. Further, Equa- 
tions (54) and (58) will be expanded in series form and only the 
first two terms retained. For a constant surface temperature this 
gives 


Pr—o: — = 1 — 1.15(mz) +.. 


4 = 1 — 1.06(mz) +... 
q 


Pr+0: 4 = 1-0.75(mz) +... 
q 
whereas for a constant surface heat flux 


= 1+ 0:780(m2) +... 


= 1 + 0.720(m2) +. 


o: 
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Pr—0: = 1+ 0.417(mz) +. 
8, * 

It may be noted that magnetic effects increase with increasing 
Prandtl number. However, this variation is not great, and for all 
practical purposes the first-order terms may be considered in- 
dependent of Prandtl number for Prandtl numbers greater than 
approximately one. It should be mentioned, however, that this 
would probably not be the case for higher-order terms. 

In addition to altering the momentum and thermal transfer 
processes, the application of a magnetic field may also alter the 
flow regime. For example, separation can occur as a result of 
the retarding magnetic force, and from reference [1] a second- 
order approximation for, the laminar separation point is found to 
be mz ~ 0.5. A second effect upon the flow regime concerns 
transition from laminar to turbulent flow, and from the stability 
analysis presented by Rossow [12] it is found that the transition 
Reynolds number is reduced through application of a transverse 
magnetic field fixed relative to the plate. 

In the present analysis it has been assumed that the induced 
magnetic field is negligible. In most applications this is a plausible 
assumption providing the electrical conductivity of the fluid under 
consideration is not large. However, for fluids with high electrical 
conductivities, such as liquid metals, the induced magnetic field 
may in some instances be appreciable. This effect is discussed by 
Cowling [13]. 
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Conical Turbulent Boundary Layer 
Experiments and a Correlation 
With Flat Plate Data 


Measurements in the turbulent boundary layer of unyawed cones for Mach numbers 
from 1 to 6 are presented. Specifically, the first measurements of total temperature 
profiles, directly determined skin friction, and local heat flux in the turbulent boundary 
layer of cone models are given. In addition, these experimental data are shown to 
justify the calculation of turbulent friction and heat transfer on unyawed cones by 
means of an incompressible plate friction law and simple auxiliary relations. 


W. S. BRADFIELD 


Senior Staff Scientist, Convair, 
Scientific Research Laboratory, 
San Diego, Calif. 


Experimental Methods 


HE EXPERIMENTAL INVESTIGATION was performed 
in the wind tunnel facilities of the University of Minnesota Rose- 
mount Aeronautical Laboratory. In Fig. 1 is shown a photo- 
graph of the 15-deg cone model with boundary layer probes in- 
stalled in the M,, = 3, 6 X 9-in. test section. An impact tem- 
perature probe [1]! was mounted on the upper surface of the cone 


1 Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division of Tae American 
Socrery or Mecuanicat Enoineers, and presented at the ASME- 
AIChE Heat Transfer Conference, Storrs, Conn., August 9-12, 
1959. 

Norse: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. Manuscript received at ASME Headquarters, April 28, 
1959. Paper No. 59—HT-6. 


Nomenclature: 


and an impact pressure probe on the lower surface at the meridi- 
onal plane. The probes were adjustable from outside the tunnel 
during operation by means of push rods projecting through the 
model skin as shown. The probe tips were approximately 5 
inches downstream from the tip of the model. 

The impact pressure probe was made of '/;-in. OD stainless 
steel tubing. The tip opening was 1.5 mils* high and 17 mils 
wide with outside dimensions of 3.2 X 20 mils. The probe was 
designed for fast response [2], the response time being about 15 
seconds. 

The locations of the probes relative to the surface of the model 
were measured with a 20-power stereoscopic microscope. 

The surface of the model was instrumented with static pressure 
taps and thermocouples [3]. Static pressure orifices were 14 mils 


21 mil = 0.001 in. 


area of surface 
local skin friction coefficient, r/q; 
average or total skin friction coefficient 


1 
f CA')dA’ 
A A 
c, = specific heat at constant pressure 
h = heat-transfer coefficient (quantity of heat per unit time 
per unit area per temperature difference) 
= thermal conductivity 
Mach number without subscript indicates local value of 
Mach number in undisturbed boundary layer 
one thousandth of an inch .. . roughly 1/40 of a millime- 
ter 
exponent in turbulent velocity profile 
c 
Prandtl number = 7 based on thermal equilibrium at 
the wetted surface 
pu? 


St = Stanton number = 


based on conditions just out- 

side boundary layer 

temperature 

velocity in the z-direction 

distance downstream from boundary layer origin on 
cone or plate 

distance perpendicular to the wetted surface 


y= 
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" 


boundary layer thickness determined by taking y = 6 
when u = 0.995% 


= boundary layer displacement thickness, 
0 Prva 
Y = ratio of specific heats c,/c, 
6 = boundary layer momentum defect thickness, 


8 
f ly 
o Pith | 


vy = kinematic viscosity u/p 
viscosity 

p = density 


= boundary layer recovery factor {2 

T. 
Subscripts 
© = upstream conditions ee 
1 = conditions locally at outer edge of boundary layer 4 
0 = stagnation conditions ce: 
e = properties based on recovery temperature (surface tem- I 
perature at condition of zero heat transfer) a 


incompressible 

either “‘transition’’ or “thermal boundary layer’ as ap- 
propriate 

conditions on surface of model 

as in R,, characteristic length used in calculating R, 

Superscripts 

* = indicates values corresponding to equation (9) 
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Fig. 1 


Installation of impact pressure and impact temperature probes on 15-deg total angle cone 
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in diameter and opened into '/j-in. OD brass tubing pressure 
leads. Orifice fittings were soldered into the model skin. Iron- 
constantan thermocouples were installed by soft soldering in 
holes drilled through the '/,-in-thick model skin. 

A combination of an upstream directed air jet issuing from the 
cone tip and a roughness trip (shown in Fig. 1) was used to insure 
the existence of a turbulent boundary layer. 

Direct measurements of skin friction were made using the 
balance system [4, 5] shown as Fig. 2. A segment (truncated 
cone) of the conical surface was supported from the central shaft 
on 6 steel balls in such a way that it was free to “float’’ with a 
clearance of 1 mil afforded at the leading edge and trailing edge of 
the floating surface. Lateral adjustment of the floating element 
was achieved by the selective adjustment of an internal parallel 
inclined wedge system. A conservative estimate of the maximum 
floating element misalignment with the conical surface was */; 
mil. However, the average boundary layer thickness for the skin 
friction measurements was approximately 120 mils giving a ratio 
of boundary layer thickness to surface step of 400. During the 
testing no waves could be detected as originating from either the 
front or rear crack when observed with a double pass schlieren 
system at high sensitivity. 

Determinations of local values of heat flux were made both from 
electrically heated segmented models [6] and from a “heat sink”’ 
cone [7]. 

Accuracy of Measurements. Random errors in temperature ratio, 
or in flow parameters determined by probing due to uncertainty 
in pressure and temperature measurement are:? 7',/7) = +1 per 
cent; momentum thickness +5 per cent to +7 per cent; dis- 
placement thickness +3 per cent to +5 per cent. 

For measurements on the model with the balance system, the 
experimental uncertainties are as follows: Reynolds number +3 
* See reference [7] for details. 
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Fig. 3 Mach number and total temperature profiles in the turbulent boundary layer of 15-deg total angle cone 


per cent; friction coefficient +3.4 per cent; Mach number 
+1.5 per cent. 

In connection with the direct determinations of skin friction 
(carried out in the blowdown wind tunnel), the pressure force on 
the floating element was determined by measurements of pressure 
on both the exterior and interior surfaces of the floating element. 
The pressure force was always less than 20 per cent of the friction 
drag for these measurements. 

The accuracy of heat-transfer measurements was determined 
from the uncertainties in the terms in the heat-balance equation. 
They are: (7, — T,) + 2 deg F; slope + 4%; AW/AA + 
1%. Thus the uncertainty inh was +5%. The corresponding 
uncertainty in Stanton number was +8%. 

Mach Number and Total Temperature Profiles. The first measure- 
ment of a cone flow turbulent thermal boundary layer total tem- 
perature profile was accomplished for this investigation. Simul- 
taneously, an impact pressure profile was obtained. The measure- 
ments were obtained at M; = 2.8 and R, = 4 X 10°. The results 
are shown as Fig. 3. 

As expected for a gas (Prandtl number less than unity), local 
total enthalpy values above the freestream value are observed 
near the outer edge of the boundary layer. The surface, in 
adiabatic equilibrium with the boundary layer flow, stabilized at 
T./Tu\y-0 = 0.928. This temperature ratio is one per cent lower 
than the recovery factor, according to the relation ¢ = (Pr)'/*, 
would indicate. 

Comparison [7] of the present data with existing two dimen- 
sional, nearly adiabatic data [8, 9], indicates that the cone 
boundary layer total temperature profile is similar to two-dimen- 
sional profiles in the y/6, versus To(y)/T presentation. Further, 
indications are that similarity with respect to Reynolds number 
variations is also obtained over a moderate Reynolds number range 
at M,; = 3. 
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Turbulent Boundary Layer Velocity Profiles. Velocity profiles may 
be calculated from the Mach profile number and free stream total 
temperature. A set of turbulent velocity profiles from a cone at 
M, = 3.7 is presented as Fig. 4. The data were obtained from a 
fifteen-degree total angle cone at zero yaw and in thermal equi- 
librium with the flow. A combination roughness (sand) and air 
trip [4, 5] was used to produce fully turbulent boundary layer 
for three of the five sets of data. One set was obtained with 
sand trip only and one with air trip only. 

Just as in the case of existing flat-plate data over a relatively 
wide range of Reynolds number and Mach number, the cone 
data are well represented by '/;-power profiles for the outer 
region in combination with linear profiles assumed for the laminar 
sublayer. 

An Engineering Prediction of Turbulent Skin Friction and Heat Transfer 
on a Cone at Zero Yaw at Supersonic Speeds. It was indicated by 
the results of the present investigation [7] that velocity profiles 
from the turbulent boundary layer of both plate and cone at the 
same M;, are congruent when presented as u/u; versus y/#. In 
other words 


a(y/0)™ (1) 


for both plate and cone where a has the same dependence on 
Mach number for both cases over the range of Mach number in- 
vestigated. It was found that n = 7 for the moderate range of 
Reynolds numbers investigated. 

If it is now assumed that a friction law of form C, = B(R5)* 
holds, this law together with equation (1) results in 


-2 .4 -6 


u/uy 
Fig. 4 Turbulent velocity profiles on a cone 
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for both cases. It is further assumed that B, a parameter depend- 
ing on Mach number, has the same dependence on Mach number 
for both cases. It follows that the friction coefficients will be 


related by 
C, +1 
@ 3 
(3) 


Here and in what follows, a bar over a quantity indicates a 
two-dimensional (plate or axial flow cylinder) boundary layer 
quantity. 

For the plate, 


= 


while for the cone 


Cc 
(5) 
2 dz x 
Again, as a consequence of the congruence of velocity and tem- 
perature profiles for plate and cone, the momentum defect thick- 
ness for both cases has the same dependence on M;. That is, 


= 2) (6) 


Using equations (2) and (6), equations (4) and (5) can be inte- 
grated. These results together with equation (3) give 


2 
E are 


(7) 
at a specified Mi, n, u/v;, and z = %. This development, which 
explicitly considers variable Mach number, yields precisely the 
same result as that obtained by Gazley [10] for the case of con- 
stant properties. 

It is general practice to obtain heat-transfer coefficients from 
skin friction coefficients for the turbulent boundary layer by use 
of the Reynolds analogy (see, for example, reference [11], p. 107). 
For compressible flow plate experiments, a form of the Reynolds 
analogy modified to account for property value variation has 
proved useful [12]. This relation, proposed by A. P. Colburn 
[18], for incompressible flow, is 


n+1 


c 
St = (8) 
where the Prandtl number (Pr*) is evaluated at the reference 
temperature 
T* = 7, + 05(T, — + 0.27, — (9) 
proposed by Eckert [12]. 


If it is assumed that equation (8) holds for cone flow as well as 
for plate flow, then 


(10) 
if T* = T*. 


For n = 7, equation (7) yields for local values of Stanton num- 
ber and skin friction 
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for T* = T* and R, = R, independent of Reynolds number and 
Mach number over the range covered by the experimental data 
discussed here. Note that for 5 < n < 10, equation (7) shows 
1.23 < C,/C, < 1.13. Thus the effect of variations of n on the 
transformation relation is small in this range. 

For obtaining the relation between mean skin friction and heat 
transfer on plates and cones, the relations 


2 
Cp = — | and Cy=— | 
Jo z* Jo 


are employed. It follows from the above that for average fric- 
tion and heat-transfer coefficients 


n+l 


Cy Bs n+1 


forn = 7. The effect of variations of n on the transformation re- 
lation equation (12) is negligible; for 5 <n < 10, 1.065 < C,/Cp 
< 1.035. 

Correlations of Direct Measurements of Turbulent Skin Friction and Heat 
Transfer. The results of the measurements at M; = 3.7 are shown 
as Fig. 5 in comparison with the data of Coles [13] on a flat plate 
in the form C, versus (R, — &,,) where the dashed line and the 
line through the cone data are drawn parallel to the faired line 
representing Coles’ data. R,, is the Reynolds number based on 
the distance between the measuring station and the beginning 
of transition as determined from schlieren observations. The 
dashed line is higher than the plate data by 18 per cent in accord 
with equation (11). The cone data are about 10 per cent higher 
than predicted by equation (11). 

Fig. 6 is adapted from Fig. 2 of reference [14]. It shows the 
relation between the heat-transfer data (at constant surface tem- 
perature) from two cones and a secant ogive and from the flat 
plate data of Pappas [15]. The mean line of the cone flow data is 
higher than the corresponding mean line of the flat plate data and 
parallel to the flat plate data mean line. The lines were obtained 
by least squares procedure from the data points. As predicted 
by equation (11), the mean line of the cone data is about 18 per 
cent higher than the plate data mean line. 

Fig. 7 shows a comparison between measured heat transfer and 
transformed skin friction measurements on cones at M, = 3.7. 
Stanton numbers are calculated from measured skin friction co- 
efficients by the Colburn relation, equation (8), where the 
Prandtl number is evaluated at the reference temperature pro- 
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Fig. 5 Comparison of direct measurements of turbulent skin friction on 
a cone and a plate at M; = 3.7 
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posed by Eckert [12], equation (9). The agreement is satisfac- 
tory. 

Fig. 8 shows a comparison between heat-transfer results and 
skin friction results for cones and plates over the present Mach 
number range. The data are evaluated at R, = Kk, = 3 X 10+. 
In order to provide a comparison of results on a T* = 7* basis, 
all Stanton number data were adjusted to 7T,, = 7, using the 
theory of Van Driest [17] (the same procedure was used in com- 
paring the data of Fig. 7). The maximum adjustment required 
amounted to 5 per cent of the measured value. The average ad- 
justment required was 3 per cent of the measured values. 

In nondimensionalizing the Stanton numbers and skin friction 
coefficients, incompressible flat plate theory was used in both 
cases. Values of Stanton number (incompressible) were calcu- 
lated utilizing equations (8) and (9). The reference temperature 
value used to calculate St; was calculated using the wind tunnel 
adiabatic wall temperature at each Mach number. 

The vertical bars through many of the data points in Fig. 8 
indicate the maximum scatter of the data at that Mach number. 
The symbol location in each case represents the mean value of the 
data. Where no vertical bar appears, the scatter of the data is 
represented by the symbol size. 

The solid curve is the experimentally determined mean curve 
through the plate data points. Obviously, this relation should 
not be expected to hold for flow around cones at subsonic speeds 
or for cones with detached shocks. 

The data correlations indicated by Figs. 5 through 8 are re- 
garded as justification for the calculation of turbulent friction 
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Fig. 7 Comparison of heat transfer on a cone with skin friction on a 
cone at M; = 3.7 
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Fig. 8 A comporison between heat-transfer and skin-friction measurements on plates, axial flow cylinders, and cones 


and heat transfer on unyawed cones at supersonic speeds by 
means of an incompressible plate friction law and the simple 
auxiliary relations given as equations (8), (9), (11), and (12). 

As a summary of this procedure, suppose it is desired to calcu- 
late the local heat flux at the surface of a cone at Mach number 
M; and R, at the outside edge of the boundary layer. Starting 
with the Schultz-Grunow local skin friction law 


0.370 


where R, = R,, the corresponding incompressible plate Stanton 
number §t, is calculated by equation (8). The compressible flow 
flat plate St is obtained by multiplying St, by St/St, from Fig. 8. 
The cone Stanton number St is obtained by increasing St by 18 
per cent in accordance with equation (11). The same procedure 
is followed for average values except the relation between cone and 
plate values is given by equation (12). 

Summary and Conclusions. The first measurements of the total 
temperature distribution through a turbulent conical boundary 
layer are presented. 

A comparison of directly measured local skin friction values on 
the cone with those on the plate at constant Mach number indi- 
cates that cone skin friction is higher but that it has the same de- 
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pendence on Reynolds number as plate skin friction (Fig. 5). A 
comparison of local heat-transfer values on cone and plate 
similarly shows higher heat transfer on a cone at the same M, 
and Reynolds number based on property values just outside the 
boundary layer (Fig. 6). 

It is confirmed that the modified Reynolds analogy of Colburn, 
equation (5), together with the reference temperature of Eck- 
ert, equation (6), permits the calculation of cone heat-transfer 
values from known cone skin-friction values (Figs. 7 and 8). 

Finally, an analysis is presented which predicts that skin-fric- 
tion and heat-transfer values on unyawed cone in supersonic flow 
should be approximately a constant times the corresponding flat- 
plate values and this prediction is confirmed by experiment. 
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The Influence of Free Stream Turbulence 


R. A. SEBAN 

Professor of Mechanical Engineering, 
University of California, Berkeley, 
Calif. Mem. ASME 


on the Local Heat Transfer From Cylinders 


Local heat-transfer coefficients and recovery factors are presented for three different cyl- 


inders in a two-dimensional subsonic air flow, with emphasis on the effect of screen-pro- 
duced turbulence on these quantities. The increase in turbulent intensity so realized 
produced larger local heat-transfer coefficients, in a way dependent upon the location on 
the cylinders, through a direct increase in the heat transfer to the laminar boundary 
layer, through an earlier transition to turbulence, or through an alteration in the charac- 
ter of the separated flow. Alternatively, recovery factors were affected less, being in- 
variant with respect to the turbulent intensity for attached boundary layer flow, but 
demonstrating large changes in those separated flow regions for which increased free 
stream turbulence produced substantial changes in the nature of the separated flow. 


Introduction 


HAS BEEN KNOWN for a long time that the heat 
transfer from a body to a fluid flowing over it is increased if the 
turbulence level of the stream increases, and recently the experi- 
mental results demonstrating the effect have been augmented, 
particularly for over-all heat transfer from circular cylinders 
placed transverse to an air stream. For systems like this, the 
gain in the over-all heat transfer is due to a number of effects, 
such as an increase in the rate of heat transfer to the laminar layer, 
an earlier transition to turbulent flow with its higher heat trans- 
fer, and movement of the separation point and alteration of the 
character of the flow in the separated region. Giedt [1]! was the 
first to achieve a separation of these effects by measuring local 
heat-transfer coefficients on a circular cylinder placed normal to 
an air stream, and these results demonstrated in particular the 
increased rate of heat transfer to the laminar boundary layer on 
the front of the cylinder. Kestin [2] has recently presented ad- 
ditional results on the effect of screen-produced turbulence on the 
over-all heat transfer from circular cylinders in air stream, to show 
that slight increases in the turbulent intensity produce substantial 
increases in the heat transfer, with a slight continued augmenta- 
tion as the turbulent intensity is increased further. Sato [3] has 
presented results for the over-all heat transfer from spheres, which 
apply to higher turbulence levels and indicate a more sustained 
increase as the turbulence level rises. 

The purpose of this paper is the presentation of additional evi- 


1 Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division of THe AMERICAN 
Socrety or MecnanicaL ENGINEERS and presented at the ASME- 
AIChE Heat Transfer Conference, Storrs, Conn., August 9-12, 1959. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, April 
30, 1959. Paper No. 59—HT-3. 
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dence on the effect of screen-produced turbulence on the local 
heat transfer from cylinders. These results are not in as homo- 
geneous a form as might be expected, for they were obtained 
largely as an adjunct to other experimental programs, but they 
do have the co-ordinating element of all having been obtained in 
the same wind tunnel with similar screens. They are also so 
far unique in embracing air speeds of sufficient magnitude so that 
recovery factors could be determined. But these aspects imposed 
a disadvantage, not discovered in detail until near the end of the 
program, in that apparent low frequency oscillations in the wind 
tunnel rendered partly ambiguous the relative turbulence increase 
attained with the screens. While, therefore, no quantitative 
specification of the increase in the heat transfer relative to the tur- 
bulence level has been achieved, the results in one sense indicate 
that possibly none should be expected. 


Apparatus and Procedure 


The models for which results have been obtained were a 1:3 
elliptical cylinder, with 6-in. major axis, a 1.25-in. diameter circu- 
lar cylinder, and a 1.87-in. diameter circular cylinder fitted with a 
downstream splitter plate. All were made of bakelite, with 
Nichrome ribbons on their surfaces for heating by electric dis- 
sipation. The elliptic cylinder described previously [4] con- 
tained three longitudinal ribbons, 0.002 in. thick and 1 in. wide; 
the 1.87-in. circular cylinder was wound helically with such a 
ribbon, and the 1.25-in. cylinder was wound helically with a 
similar ribbon 0.001 in. thick and 0.50 in. wide. Thermocouples 
were located beneath the ribbons and pressure taps in small spaces 
between them, and on the circular cylinders all instrumentation 
was situated on one cylinder generator and local results were ob- 
tained by altering this location by rotation of the cylinder. All 
the models were mounted centrally in and spanned the 6-in. 
width of the 6 X 9-in. rectangular test section of the wind tunnel. 

Data were obtained at air velocities from 150 to 350 ft/sec, with 
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stagnation temperatures near 90 deg F and atmospheric stagna- 
tion pressure, with heating sufficient to give temperature dif- 
ferences of from 15 to 20 deg F between ribbon surface and stag- 
nation temperature. The results are presented conventionally as 
pressure coefficients or local free-stream velocity deduced from 
the surface pressure observations, and heat-transfer coefficients 
determined from the heat dissipation in the ribbons and the effec- 
tive temperature difference. These were evaluated directly from 
the data without correction and the errors due to the omission of 
the effects of internal conduction and surface radiation are small. 
They vary with the model and the conditions of observation but 
as a categorical specification it can be stated that all recovery 
factors shown are accurate to 2 per cent and heat-transfer co- 
efficients to 5 per cent. 

The screens were hardware mesh having the dimensions given 
in Table 1, where M is the distance in inches between wire centers 
and 6 the wire diameter. 


Table 1 

Mesh /in., 

nominal MM (in.) 5 (in.) M/s 
2 1/2 0.108 4.63 
3 11/32 0.062 5.37 
4 1/4 0.050 5.00 
8 1/8 0.028 4.47 
16 1/16 0.013 4.82 


These screens were mounted at the end of the nozzle contrac- 
tion, which was the beginning of the test section, and the leading 
edges of the models were located at the following distances down- 
stream of the screens: Elliptic cylinder, 9 in.; 1.25-in. cylinder, 
8*/, in.; 1.87-in. cylinder, 8 in. 

Turbulence determinations were made with a Flow Corpora- 
tion HWB hot-wire system, using a 0.35-mil diameter wire. Total 
output was determined on a Ballantine Model 320 rms meter and 
a Hewlett-Packard Model 300A wave analyzer was used for spec- 
tral determinations. 


Results With the Elliptic Cylinder 


Fig. 1 shows the effect of turbulence produced by the insertion 
of the 2 and 4-mesh screens on the heat transfer from this model, 
with the magnitudes for clear tunnel operation being shown as the 
curve which correlated these results in the earlier report [4]. 
There is apparent from the figure a negligible effect on the free- 
stream velocity distribution and this is so because of the relatively 
small change in the character of the small separated region that 
existed at the end of this model. 

Two effects appear in the results for the local heat-transfer co- 
efficient. The value at the stagnation point is increased in the 
order of 30 per cent and an increase is maintained in diminishing 
magnitude with increasing distance from the stagnation point. 
Also, the increased turbulence promoted earlier transition to tur- 
bulent flow in the boundary layer, after which it can be seen from 
the figure that the heat-transfer coefficients are predicted satis- 
factorily by the Colburn equation, hs/k = 0.027 (us/v)*. In 
the case of the elliptic cylinder, it is this change in the point of 
transition that produces the major effect of the turbulence on the 
over-all heat transfer. With the turbulent separation, occurring 
farther downstream than the laminar separation evidenced with 
clear tunnel operation, the heat-transfer coefficients at the rear of 
the cylinder are reduced slightly. 

Recovery factors are not shown in the figure; insertion of the 
screens gave values typical of turbulent flow after the indicated 
transition points, while in the upstream laminar region the magni- 
tudes were unaffected by the insertion of the screens. 

Other results of the same general nature, presented elsewhere 
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[5] were obtained with the cylinder at +6 deg angle of attack. 
These support the evidence as presented in Fig. 1. 


Results With the 1.25-In. Circular Cylinder 


The results for the circular cylinder are presented in terms of 
the pressure coefficient, recovery factor, and correlating group for 
the heat-transfer coefficient. Here the recovery factor differs, 
being referred to the upstream rather than to the local velocity. 
Fig. 2 shows the performance with a clear tunnel, and makes the 
comparison to values predicted by the method of Seban [6] in 
the laminar flow region. The correspondence between predicted 
and measured values is good. Beyond the separation point the 
rather uncorrelated behavior of the local heat-transfer coefficient 
is typical for this range of Reynolds numbers, and there the re- 
covery factors decrease continually to a negative value at the rear 
stagnation point. 

Only the 4-mesh screen was used with this cylinder, to give 
the results shown in Fig. 3 which reveal for this bluff body a more 
drastic alteration than with the elliptic cylinder of flow and heat- 
transfer characteristics. The drag is reduced due to an apparent 
delay in separation until 120 deg and the velocity distribution 
over the front of the cylinder is slightly altered. In the region of 
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Fig. 1 Heat transfer on the elliptical cylinder 
Curves a, d, Performance with clear tunnel. The dashed line near "a" ap- 
plies for the higher Reynolds number. 
b, c, Predictions for turbulent boundary layer flow for the high and 
low Reynolds numbers. 


Run Mesh Uw, ft/sec 
1 2 8.14 XK 105 284 
4 4 8.14 284 
5 4 4.63 160 
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laminar flow there exists the increase in local heat-transfer co- 
efficient similar to that found on the elliptic cylinder, now more 
extensive in effect. Remarkably, the recovery factor is unaltered, 
and is predicted as satisfactorily as before for the laminar flow 
region, and this laminar region appears to persist at least to 90 deg. 
There the heat-transfer coefficient rises, but to values less than 
would be predicted by the Colburn equation, while the recovery 
factor diminishes to values below any that would be associated 
with known boundary-layer effects. The behavior here is sug- 
gestive of a laminar separation near 90 deg, a transition to turbu- 
lent flow and a subsequent reattachment near 110 deg, with a 
new separation in the region of 120 to 130 deg. The final sepa- 
rated flow is so altered that the low recovery factors typical of 
the operation in the clear tunnel are no longer obtained. 

Fig. 3 also contains curves which define the low speed result 
of Giedt [1] obtained with a larger cylinder in a different wind 
tunnel behind a rope net of the same relative geometry as the 
screen. Similarity is evident except for the difference in pressure 
coefficient behind the cylinder, an alteration due to wake blockage 
in the closed tunnel used with the 1'/,-in. cylinder. 


Results With the 1.87-In. Circular Cylinder 


The 1.87-in. cylinder was examined with an extended range of 
screen sizes but, because of defects in the model, the results in 
their totality are somewhat less acceptable than those so far ex- 
hibited. This was caused by surface irregularities connected with 
the relatively larger helix angle of the wider ribbon winding on 
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Fig. 2. Results for the 1.25-in. circular cylinder 


Reynolds numbers, v.d/v, are indicated in the figure. Here the recovery 
factor, r, is taken a3 (T, — T)/u..?/2c>. 
Curves a, Predicted value for laminar boundary layer heat transfer. 

b, Recovery factor predicted from a local value of 0.84. 
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this cylinder, which irregularities, though being less than 0.001 
in., were sufficient to disturb the flow near the separation line, 
particularly at high speeds. The use of a downstream splitter 
plate, 0.09 in. thick and extending downstream 12 in. from the 
rear stagnation point, ameliorated this difficulty in part, and 
the exigencies of the experimental program at that time dictated 
the use of this cylinder. The results are pertinent here in be- 
ing typical at least in the region of laminar flow, to demonstrate 
further the effect of screen size. Further details on the operation 
of this cylinder in the clear tunnel are given in Reference [8]. 
Fig. 4 shows results for this cylinder with the clear tunnel and 
with the 4-mesh screen, at velocities from 190 to 380 ft/sec and 
irregularity is demonstrated directly by the pressure coefficient 
at an angle of 70 deg with clear tunnel operation. Here there is 
scatter in the results which are themselves averages of a spanwise 
variation of 10 per cent per inch in the velocity, a variation which 
was this serious only in this locality. Some of this variation is 
shown also in the heat-transfer correlation, which in the laminar 
flow region fails to resolve the performance at different speeds, an 
effect basically consequent on changes in the distribution of the 
local free-stream velocity. The two higher speed runs do, how- 
ever, agree fairly with the prediction, made for the corresponding 
pressure coefficient. The recovery factors, shown here on a local 
basis, agree with the expected value of about 0.84. Because of 
the poor behavior in that region, separation is ill defined but ap- 
parently occurred at about 80 deg. Downstream of separation 
the heat-transfer coefficient varies considerably and the correla- 
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Fig.3 Resuits for the 1.25-in. cylinder with 4-mesh screen 


The representation is the same as that of Fig. 2. 
Curves o, b, Heat-transfer and pressure coefficients from Giedt [1], for 


cylinder downstream of rope net at a Reynolds number of 
170,000, 

¢, Predicted value of the recovery factor, from a local value of 
0.84 for @ < 90 deg and 0.89 for @ > 90 deg. 

d, Heat-transfer prediction. 
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Fig. 4 Results for the 1.87-in. cylinder with splitter plate with clear 
tunnel and 4-mesh screen 
Curves a, b, Predictions for laminar flow from the velocity distribution 
for 280 ft/sec, clear tunnel and with screen. 
¢, Prediction for turbulent flow, Run 37. 
Experimental results 


Run 55 30 35 44 37 47 
Mesh one 4 4 4 
Ue ft/sec 192 288 
vod 19-5 171 250 3.24 250 2.80 


tion fails. The behavior is, however, typical of this system, in 
which the separated flow is ultimately reattached to the splitter 
plate. 

Insertion of the sereen produced a marked stabilization of the 
flow over the forward part of the cylinder, with resultant im- 
provements in the correlation of the pressure coefficient and the 
group containing the heat-transfer coefficient. Separation was re- 
tarded and there was in consequence a sufficient change in the 
pressure coefficient to require a new prediction of the heat-transfer 
coefficient to enable a proper evaluation of the effect of turbulence 
in the laminar region. The increase is seen to be of the same order 
as experienced with the 1.25-in. cylinder, being of 30 to 40 per cent 
at the stagnation point and diminishing with downstream dis- 
tance. The recovery factors remain the same and they indicate 
laminar flow up to the 100 deg position. Thereafter the heat- 
transfer coefficient increases and the recovery factor decreases in 
a way that suggests flow separation rather than the occurrence of 
any turbulent flow. 

Fig. 5 shows additional results for the forward region of this 
cylinder, to illustrate primarily the effect of changing the screen 
size. There is little trend with screen size, except that the results 
for the 16-mesh screen are substantially closer to the clear tunnel 
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Fig.5 Heat transfer from 1.87-in. cylinder and plate with various screens 
Experimental results 


Run 50 44 51 53 55 

Mesh 2 4 8 16 ne 

Ua ft/sec 260 194 194 195 192 

vad» 19-5 2.25 1.70 1.70: 1.70 1.70 


performance with respect to the heat transfer, and for it the pres- 
sure coefficient is almost identical to that for operation in the 
clear tunnel. 


Turbulence Produced by the Screens 


To the obvious increases in heat transfer realized by the inser- 
tion of the screens there is immediately associated the question 
about the nature of the turbulence so produced, and this was de- 
termined partially in terms of the intensity of the longitudinal 
component and in terms of the spectral distribution of that quan- 
tity. Adequate measurements of even the first of these quantities 
was realized only comparatively late in the program, and these 
revealed that the turbulent intensities in the clear tunnel were un- 
expectedly high. In consequence, the introduction of the larger 
mesh screens provided very small increases in the observed turbu- 
lent intensity. Spectral determinations demonstrated some of 
the causes of this, and as shown in Figs. 6 and 7 in terms of nor- 
malized spectrum functions, they indicate unusual low frequency 
contributions with clear tunnel operation, together with a discrete 
input at a wave number near 17. The source of these exceptional 
turbulent contributions has not been established. The spectra in 
Fig. 7 show that at least the discrete input at wave number 17 
was eliminated by the screens, though the low wave number 
energy in part remained. 

The intensities presented in Table 2 are to be viewed in this re- 
gard, where as a specific example the introduction of the 4-mesh 
screen led to turbulent intensities of the same magnitude as given 
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Fig.6 Normalized power spectra for the clear tunnel 
The line shown is a reference line having a slope of — 5/3. The ordinate is 
f= where T. is the intensity per unit wave number. 


by the clear tunnel. As a contrast to emphasize the difference in 
intensities at higher wave numbers, Table 2 contains also a few 
values giving the intensity above an arbitrarily selected wave 
number of 40 and these show the substantial and expected effect 
of the screens in this spectral region. 


Table 2 


Intensity 100 


1 


2 


1. 
1. 
3. 


While the intensities existing in the clear tunnel are subject to 
no absolute comparison, those produced by the screens can be 
appraised in terms of the known results obtained at lower speeds 
behind rod grids of similar relative spacing. This is done in Fig. 
8, to show that at the lower speeds the expected intensities were 
attained, while at the highest speeds the values attained were 
much too low. It is estimated that this latter difficulty was 
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Fig. 7 Normalized power spectra for various screens 


associated with possibly inadequate compensation, whereby the 
energy at the higher wave numbers was inadequately sampled. 

Another appraisal of the screen data is possible through com- 
parison of the scale with the mesh size, with the scale determined 
as L = (4/2)F(0). Because of the low frequency energy that was 
still existent with the screens, the “zero’’ value is estimated from 
the flat portion of the spectra, with the results shown in Table 3. 
The values of scale are of the order of the mesh size. 


Relation of the Heat Transfer to the Turbulence 


It is evident from the results that the heat transfer is affected 
by the free stream turbulence in a number of ways. There is a 
definite effect upon the heat transfer to the laminar boundary 
layer, upon the point of transition and extent of the subsequent 
flow, possibly upon the point of separation and certainly upon the 
character of the separated flow and the region of the body adjacent 
hereto. The cumulative effect of each of these depends upon their 
relative importance on the particular body form considered, as 
can be assessed, for instance, by comparing the present results for 
the elliptical and circular cylinders. Of the individual effects the 
increase in the heat-transfer coefficient and invariability of the 
recovery factor of the laminar boundary-layer flow are phenome- 
nologically the most definite in the present results and the rela- 
tively well-known behavior of the heat transfer in this region 
enables a more specific appraisal of the effect of increase in the 
turbulence level. 


Table 3 


Screen mesh 
M inches 
F(“0O" )it/rad 
L inches 
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10¢ 2 4 686 2 «668 Toh, 
10° 10! 10? 10° 10° 
q Vo ua 
Total for > 40 
150 250 340 190 340 
Clear 5 1.0 0.90 0.10 0.07 
BS 4 mesh 6 1.1 1.0 0.94 a 
3 1.6 1.4 1.0 
es 2 3 2.6 1.0 1.0 355 
a 
2 3 4 8 16 Clear 
0.50 0.34 0.25 0.12 0.07 
ats. 0.017 0.10 0.008 0.005 0.002 > 0.20 
0.32 0.20 0.15 0.10 0.04 > 3.8 
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Fig. 8 Turbulence levels for various speeds and positions 


The results are mostly for three speeds, of the magnitudes 145, 190, and 
340 ft/sec. At any position the results for the intermediate speed are 
distinguished by solid points and the results for the higher and lower speed 
are open points, respectively, above and below the solid points. Here the 
distance x is the distance downstream of the screen position. The line is 
Batchelor's [7] specification for M/S = 5.5. 


A uniform effect of the introduction of the screens was an in- 
crease in the heat transfer in the region of laminar flow. This was 
& maximum at the stagnation point and diminished with increas- 
ing distance downstream, to suggest a possible strong dependence 
upon pressure gradient. There is a variable sensitivity to turbu- 
lence level, with a slight increase over the clear tunnel value pro- 
viding an increase of the order of 30 per cent in the coefficient at 
the stagnation point, with additional increases in turbulence pro- 
viding only slight additional increases in the heat transfer. The 
failure to make this relation quantitative is an unfortunate aspect 
of the present result which is due, as already noted, to the ab- 
normally high total intensity of clear tunnel operation. This 
appears to be due primarily to low frequency components which 
apparently exercised no influence on the laminar flow, so that in- 
ferentially some special spectral dependence is also anticipated 
for the ultimate rationalization of the increase in the heat-transfer 
coefficient which is realized by the insertion of the screens. 

Finally, there arises the question of theoretical explanation of 
the increase, for which two general arguments might equally 
apply. It has been suggested that Gértler type instability might 
occur near the stagnation point where the onstream flow turns in 
& concave way to pass over the body surface. The wave length 
of such an instability would be of the order of 10 momentum 
thicknesses, which for the operating conditions embraced by these 
results would be from 0.04 to 0.10 in., at least of the order of what 
might be expected as a scale produced by the screens, though cer- 
tainly not by the clear tunnel. In this sense such oscillations are 
possibility. 
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Alternatively, there is the possibility of applying a theory such 
as that of Lin [8] and accounting in it both for temporal and 
spatial dependence of a periodic perturbation of the free stream 
velocity. This indicates that a periodic oscillation gives rise to an 
effective momentum transport increasing in importance with 
wave number and also with intensity. Presumably the heat trans- 
fer would be augmented also, but considerations along this line 
really indicate only an increased influence of higher frequency 
oscillations, an effect which is, of course, inferred also from the 
nature of the present results. 

Finally, it must be noted that the as yet unavailable theory 
which would be descriptive of the presently exhibited effects on 
the laminar boundary layer would need to predict also the inde- 
pendence of the recovery factor from the free-stream turbulence, 
which would require a special kind of similarity between the in- 
creases in the heat-transfer and the friction coefficients. 


Summary 


This exhibit of the effect of screen-produced turbulence on the 
local heat transfer from three different cylindrical bodies reveals 
varied effects, due to alterations in the viscid flow which produced 
changes both in its local character and in the flow field over the 
body. The over-all effect of increase in the turbulence is an in- 
crease in the over-all heat-transfer coefficient, due to the totality 
of the following specific changes. 


1 The coefficient for the laminar boundary layer is increased, 
with a maximum gain at the point of largest pressure gradient and 
a@ minimum at the smallest. The local recovery factor, however, 
is unaffected by the alteration of the turbulent intensity. 

2 As expected, the transition Reynolds number is reduced to 
lower values, but the results for the elliptical cylinder indicate no 
effect of the turbulence on the heat transfer to the region of tur- 
bulent boundary layer flow. 

3 The results from the circular cylinders reveal no effect on 
the separation point per se but the alteration of the wake flow is 
sufficient to alter substantially the position of separation by 
changes in the free-stream velocity distribution. 

4 If the flow remains laminar at separation, as it apparently 
did on the 1.25-in. circular cylinder, then the increase in free 
stream intensity promotes earlier transition in the separated re- 
gion, with probable reattachment before a final separation. This, 
of course, reduces the drag in the classic effect whereby bodies of 
this type have been used as turbulence indicators. The great 
attendant alteration of flow conditions in the rear of the cylinder 
alters completely the distribution of local heat-transfer and re- 
covery factor in such a complicated way as to be discerned only 
by observation of the results as pictured here. Of the many 
effects a noteworthy one is the complete elimination at the rear of 
the cylinder of the low recovery factors that are characteristic of 
subcritical fiow. 


A major portion of this work was accomplished under the spon- 
sorship of the United States Air Force, under AF 33(616)-348 and 
because of the duration of the effort there were many participants. 
Among them special mention is due to A. Levy, A. Emery, P. 
Gildea, and D. Doughty. 
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This paper represents an additional contribution to the grow- 
ing background of information concerning the influence of grid- 
induced turbulence upon the local thermal transfer from cylinders. 
Apparently an important difficulty in obtaining satisfactory agree- 
ment between different investigators is the lack of understanding 
of the details of the turbulent degradations taking place in the 
wake of a grid system. The author appears to recognize this, 
and the information portrayed in Fig. 8 illustrates some of the 
difficulties in relating turbulence level as measured with a hot 
wire anemometer to the geometry of the screen or grid for a range 
of flow conditions. The agreement with the work of Giedt [1] is 
as good as should be expected, and the scatter of the experimental 
points in the after portion of the cylinder is typical for such 
measurements. 

The author is to be commended for his frankness in regard to 
the inadequacy of the data on the 1.87-in. cylinder. This labora- 
tory also has found that small surface irregularities cause un- 
usually large variations in the thermal transfer, particularly in the 


? Chemical Engineering Laboratory, California Institute of Tech- 
nology, Pasadena, Calif. 


after part of a cylinder. The increased stability of the flow pat- 
tern around the cylinder in the wake of a grid as compared to 
open channel shear flow appears typical of the behavior of 
cylinders and spheres. 

The indicated increase in thermal transfer in the forward seg- 
ment of the cylinder as a result of the increase in turbulence in 
the wake of a grid is not greatly different in total effect from that 
found by other investigators for both cylinders and spheres [I, 
2, 3]. However, this laboratory found a progressive increase in 
the rate of change of local thermal transfer with increase in the 
apparent level of turbulence. The markedly differing trends thus 
indicated as to the effect of apparent level of turbulence upon the 
thermal transfer in the forward segment of the cylinder and in 
spheres [3] may well be related to the microscopic aspects of the 
degradation phenomena in the wake of grids which at the present 
time do not appear to be well understood. 

In closing, the author is to be complimented on bringing to an 
effective conclusion a difficult experimental program. 


Author's Closure 


Professor Sage correctly indicates the question that remains 
about the relation between the increase in the local heat transfer 
in the region of laminar flow and the turbulent intensity of the 
stream. In his work, results have been obtained which indicate a 
simple functional relation between the increased heat transfer 
and the turbulent intensity. A survey of those given here re- 
veals no such relation. 

For example, if results at the stagnation point are chosen and 
these are represented as a function of the intensity, there is an 
indication of a trend of increasing heat-transfer coefficient with 
increased turbulence level. Considerable scatter exists and the 
question of a relation is confounded by the abnormally high tur- 
bulence level of the “clear’’ tunnel. Removal of this contradic- 
tion by using intensities over a wave number of 40, as indicated by 
the last columns in Table 2, does not minimize the scatter of those 
results representing increases of the order of 30 per cent in the heat- 
transfer coefficient. Hence, as indicated in the paper, the present 
results establish no functional relationship in this regard and ad- 
ditional experimental results are needed to clarify this aspect of 
the problem. 
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arbitrary circumferential heat flux. 


Heat Transfer to Fully Developed Laminar 
Flow in a Circular Tube With Arbitrary 
Circumferential Heat Flux 


A solution is presented for fully developed laminar heat transfer in a circular tube with 
Examples included indicate that the influence of 


circumferential heat-flux variation on wall temperatures can be quite significant and 


provide some insight into the nature of the effects. 


wall temperatures for any arbitrary peripheral heat-flux variation for fully developed 


Introduction 


HE CIRCULAR TUBE has been the subject of a con- 
siderable amount of heat-transfer analysis and experimentation. 
The theories are now sufficiently advanced that it is possible to 
predict the heat-transfer performance in either laminar or tur- 
bulent flow, with arbitrary axial wall temperature [1, 2] or heat 
flux (3, 4]. The essential restrictions of these theories are: 1 
That the velocity profile be fully established at the entrance of 
the tube; 2 that the fluid properties be constant; and 3 that 
the teranerature field is radially symmetric, i.e., the heat flux and 
wall temperature do not vary around the circumference of the 
tube. For the special case of constant wall temperature, the 
effects of a developing velocity profile have been investigated 
(5, 6). The effect of fluid property variation has also been 
studied [7]. The present paper treats fully established laminar 
flow in a circular tube with constant axial heat input but variable 
circumferential heat flux. Although restricted to a particularly 
simple case, the present analysis provides some general insight 
into the nature and magnitude of the effects of circumferential 
heat-flux variation. 

The method of analysis is as follows: First, a solution is ob- 
tained for the case of a tube with constant heat flux over a por- 
tion of its circumference, insulated over the remainder; reduction 
of this solution to the limiting case of heating over a differential 
segment, of circumference leads to a superposition kernel which in 
turn allows the arbitrary circumferential flux problem to be 
handled. Throughout, both the temperature and velocity profiles 


! Numbers in brackets designate References at end of paper. 
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laminar flow under the restriction of constant axtal heat input. 


are assumed to be fully established, the heat flux is taken as con- 
stant in the axial direction, fluid properties are assumed constant, 
and viscous dissipation is neglected. An example illustrates how 
circumferential heat-flux variations might influence the wall tem- 
peratures in a nuclear reactor. 


Analysis 


Under the previously mentioned assumptions, application of an 
energy balance and suitable rate equations to a differential ele- 
ment of the fluid leads to the following governing differential 
equation for the fluid temperature [8]: 

dt 1 O% Vv 


(1) 
dr? r Or r? a or 


The velocity profile may be shown to be [8] 
V = 2V,(1 — r**) (2) 


If the temperature profile is fully established, then the rate of 
change of temperature with length must be the same at every 
point, and therefore 

dt at 

— = — = const 3 

(3) 
If the heat input is constant in the axial direction, then an 
energy balance dictates that, for the tube heated as shown by 
Fig. 1, 


= 2roBq" (4) 


Again using the fact that the temperature profile is fully estab- 
lished, one can introduce the temperature difference as de- 
pendent variable, 


At(r, 6) U(r, 6, 2) t,,(z) (5) 


flow area, ft? 


= radius from tube axis, ft mean fluid velocity, ft/sec 
c, = specific heat at constant pressure, ro = radius of tube, ft xz = distance along tube, ft 
Btu/lb deg F r* = r/r% a = thermal diffusivity of fluid, 2a = 
C, = Fourier expansion coefficients k/pe,, ft?/hr 
h = local convective heat-transfer co- - deg F as : B = half angle of heated segment, ra- 
efficient, Btu/hr ft* deg F , = local wall temperature, deg F dians 
k& = thermal conductivity of fluid, At = t —t,, deg F An = eigenvalues 
Btu/hr ft deg F At, = t, — t,, deg F p = fluid density, lb/ft? 
Nu = Nusselt number, Nu = 2hro/k V = fluid velocity, ft/sec 6, w = angle in polar co-ordinate system 
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Fig. 1 Tube with heated segment centered at 0 = O deg 


Equation (1) may now be written as follows: 
orAt 1 1 


r* or 
Here 
a = 4Brq" (7) 


For the problem of Fig. 1, Equation (6) is to be solved subject 
to the boundary conditions 


ot 
k— =q’ 


ot 


The analysis is simplified by the change of dependent variable, 


+8 (8a) 


B<0<2r (8b) 


r*2 


at = - (9) 


Substituting (9) in (6) one sees that the variable F must satisfy 
Laplace’s equation, 


oF 1 OF, 1 
Ore 


oF 
r*2 26? 
F must also satisfy the boundary conditions, 


oF 


=0 (10) 


(lle) 


= —(q’ro/k)(B/m) B<@0<2r (110) 


Equation (10) may be solved by separation of variables. Let 
F(r*, 6) = R(r*)-T(@). Under this assumption, the functions 
R and T must satisfy the following equations: 

R’ + (1/r*)R’ - MR =0 (12a) 


T’+NrT =0 
Equations (12) have the solutions 


(126) 
R = + 


T = B, cos + sin XO 
Because of symmetry, 
oAt 
=@ 
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(13a) 
(138) 


(14a) 


oAt 


These symmetry conditions lead to 
B, = 0 (15a) 
sin \r = 0 (156) 
Therefore the eigenvalues are simply integers, 
A (16) 
Moreover, since the temperature must be finite at the tube center, 
A; = 0 (17) 


Equation (10) is linear, and therefore the linear combination of an 
infinite number of solutions may be taken in order to satisfy the 
boundary conditions (1). The solution for F then takes the form 


C,r** cos nO (18) 


The boundary conditions are to be applied to the derivative and 
lead to the evaluation of the expansion coefficients C, in the usual 
Fourier manner. One finds the following: 


The solution now takes the form 


2q"ro 
+Q+ > sin nB-r**-cosn@ (20) 


The single remaining unknown, C), must be such that the mixed 
mean temperature is 


1 
21 
te (21) 


Equation (21) may be written as 
f, (1 — r**)r*Atd@dr* = 0 


When (20) is used to evaluate the integral of (22), the cosine 
terms all drop out in the integration over 0. The remain- 
ing terms may be readily integrated, and one finds for Cy 

Co = —148/489 (23) 


The wall-temperature difference, At,(0) = At(ro, 6), may now be 
expressed as (for problem of Fig. 1) 


(22) 


At,(9)/(q"2ro/k) = + — sin nB-cosn@ (24) 
T 


Note that this is simply the inverse of the local Nusselt number, 
and that for 8 = 7, or constant circumferential heat flux, the 
Nusselt number has the familiar value of 48/11 = 4.364. This 
completes the first part of the analysis. 

We next consider the temperature distribution produced in a 
tube heated as shown by Fig. 2. The tube is insulated around its 
periphery, except for @ between w and w + Aw, where the tube is 
heated with a constant flux q’(w). The solution for this problem 
may be obtained by a simple shift of axes. We will denote the 
solution to this problem by 


bAt(r, 0; w) 
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‘ 
Fig. 2 Tube with heated segment at 0 = w 
which denotes a small temperature difference induced by heating 
between w and w + Aw at a rate g"(w). The wall temperature 
for this case is, from Equation (24), 


48 


+ Fo 


k 


This solution can be used to obtain the solution for the 
temperature distribution in the tube heated at different rates 
over several segments of its circumference, as shown by Fig. 3. 


q’ 
A 
4 
-7 6, 


Fig. 3 Tube with different heat fluxes over different segments 


Imagine that only the m™ segment was heated at q’(w,,), and the 


remainder of the tube was insulated. The temperature distribu- 
tion would be 


6,,U(r, 9, 2; w,,) 


Since the governing differential equation (1) is linear in tempera- 
ture, the sum of M such solutions is also a solution, and hence 


M 
u(r, 0,2) = 0, 2; w,,) (26) 


m=1 


is also a solution of the differential equation. Moreover, such a 
solution satisfies the boundary conditions that 


ot 
k = q" 27 
q"() (27) 
where q’(@) represents the heat flux at any point @ prescribed as 
shown in Fig. 3. Since (26) satisfies both the differential equa- 
tion and the imposed boundary conditions, it must represent the 
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aw) 


Fig. 4 Tube with arbitrary heat flux 


solution to the problem of Fig. 3. The wall temperature is there- 
fore given by 


M M 
At(0) = 6, w,,) = | Aw,, 
m=1 k 48 


The problem of completely arbitrary circumferential heat flux, 
shown in Fig. 4, can be solved simply by letting Aw be infinitesi- 
mal, dw. The summation over the elemental segments dw in 
Equation (28) is then replaced by an integral around the circum- 
ference and the solution reduces to 


q"(w)2ro i 
The infinite series appearing in Equation (29) has the sum [9] 
7 n 2 


We now denote an integration kernel, g(@; w), by 


ll = cos n(@ — w) 
— + (31a) 


n=1 


(31b) 


It then follows that the wall-temperature distribution is given by 


= cos n(@ — w) | dw 
|z (29) 


or alternatively, using (30) 


11 
g(9; w) = In 


dw 
= (2ro/k) f w) (32) 
0 


Equation (32) represents the general solution to the arbitrary 
circumferential flux problem. Knowing the heat-flux variation 
around the circumference, g’(@), one need only perform the in- 
tegration of (32) to evaluate the wall-temperature distribution. 
In some cases the series form of the integration kernel (31a) will 
be easier to use. The closed form expression of the kernel (316) 
makes numerical integration extremely easy in problems where 
exact integration is not possible. 

It should be noted that the kernel has a singularity at the 
point @— w. At this point the series of (31a) behaves like 1/n, 
and therefore diverges. However, the integrated series behaves 
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Fig. 5 Illustration of effect of circumferential heat-filux variation on fully 
established laminar heat transfer 


like 1/n?, which converges, and therefore the wall temperature, as 
expressed by (32), is finite. At all other points the convergence is 
absolute. 


Examples 


Constant Heat Flux. For the special case of constant heat flux, 
q°(9) = q” = const, and the integration of (32) yields 
At (9) = (2reqo”/k)(11/48) = const (33) 

The Nusselt number for constant heat flux is then equal to [7] 
Nu = 48/11 = 4.364 (constant heat input) (34) 


Cosine Circumferential Flux Distribution. Consider the heat-flux dis- 
tribution shown at the top of Fig. 6, 


q"(8) = qo"(1 + b cos 8) (35) 


This might represent the sort of circumferential flux distribution 
found in a nuclear reactor. For this example the most con- 
venient form of the integrating kernel is the series form (31a). 
From (32) the wall-temperature distribution is given by 


At(@) = f qo"(1 + cos w)(2re/k) 
0 


48 n 2r 
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PARTLY HEATED, 
PARTLY COOLED ™ 
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Performing the indicated integrations one finds 


(37) 
The local Nusselt number can be calculated from the local heat- 
flux and temperature difference, 


Nu(8) (1+ /( +08) 


The results of this example for b = 1 are shown by Fig. 5. 
Note that the local Nusselt number varies substantially from the 
constant heat-input value previously determined as 4.364. This 
indicates a substantial variation of the heat-transfer coefficient 
around the tube periphery. The Nusselt number is infinite at the 
point where the wall temperature equals the mean fluid tempera- 
ture, and is negative where the wall temperature is less than the 
mean temperature. The ratio of the local wall-temperature dif- 
ference to the average temperature difference Af,., is shown in 
Fig. 5. Note that the local temperature difference is more than 
three times the average difference for 6 = 0! This indicates the 
profound influence of the circumferential heat-flux variation. It 
might be thought that the local heat flux, together with the heat- 
transfer coefficient from the constant heat-input solution, could 
be used to estimate the wall-temperature variation. But as the 
comparison of Fig. 5 shows, this is not very satisfactory because 
the flux variation markedly alters the heat-transfer coefficient. 

Fig. 6 shows the extreme temperature differences which occur 


6 
At(8) = (qo"2ro/k) +b 


(38) 
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at 6 = 0 deg and 180 deg as functions of the parameter bAlso 
shown in this figure are the values of qmin”/qmax” corresponding 
to given values of b. Note that the effect of heat-flux variation is 
surprisingly pronounced, even with relatively small circumferen- 
tial flux variations. For example, if b = 0.1 (corresponding to 
Qmin”/Qmax”" = 0.817), then Atmss/Atsyg = 1.22, which means 
that the maximum temperature difference exceeds the average 
temperature difference by 22 per cent! If one were to compute 
the maximum temperature difference using the actual heat flux 
at that point (@ = 0 deg) and the heat-transfer coefficient for 
constant circumferential flux, he would obtain Atus:/Atw, = 
1.10, and would be therefore in error by about 12 per cent. This 
example serves to illustrate how severe the circumferential 
flux variation must be before its influence becomes important, 
and provides one with some feel for the magnitude of error in- 
troduced by using the constant circumferential flux value of the 
heat-transfer coefficient to estimate peripheral temperature 
distributions. 


References 


1. 8S. N. Singh, “Heat Transfer by Laminar Flow in a Cylindrical 
Tube,” Applied Scientific Research, sect. A, vol. 7, 1958, p. 325. 

2 C. A. Sleicher and M. Tribus, “Heat Transfer in a Pipe With 
Turbulent Flow and Arbitrary Wall-Temperature Distribution,” 1956 
Heat Transfer and Fluid Mechanics Institute. 

3 R. Siegel, E. M. Sparrow, and T. M. Hallman, “Steady Laminar 
Heat Transfer in a Circular Tube With Prescribed Wall Heat Flux,” 
Applied Scientific Research, sect. A, vol. 7, 1958, p. 386. 

4 E. M. Sparrow, T. M. Hallman, and R. Siegel, ‘‘Turbulent 
Heat Transfer in the Entrance Region of a Pipe With Uniform Heat 
Flux,” Applied Scientific Research, sect. A, vol. 7, 1958, p. 37. 

5 W. M. Kays, “Numerical Solutions for Laminar Flow Heat 
Transfer in Circular Tubes,”” TR-No. 20, Contract N6onr 251, Task 
Order 6, Stanford University, October 15, 1953. 

6 R.G. Deissler, “Analysis of Turbulent Heat Transfer and Flow 


in the Entrance Region of Smooth Passages,""”’ NACA TN 3016, 
October, 1953. 

7 W.M. Kays and A. L. London, “Compact Heat Exchangers,” 
National Press, Palo Alto, California, 1950. 

8 M. Jacob, “Heat Transfer,”’ vol. 1, John Wiley & Sons, Inc., 
New York, N. Y., 1949, p. 451. 


9 K. Knopp, “Theory and Application of Infinite Series,” Blackie 
and Son, Ltd., London, England, 1928, p. 378. 


112 / may 1960 


This paper provides an interesting addition to the heat-transfer 
literature on pipe flows. The illustrative example shows that 
circumferential heat-flux variations can be important. However, 
if the flow were turbulent rather than laminar as assumed in the 
paper, then the effects of the imposed circumferential variation 
would very likely be smaller. 

The analysis was carried out by first finding the response to an 
isolated heating at 6 = w and then generalizing this result by 
superposition to provide the response to an arbitrary circumferen- 
tial heat flux q’(@). This is certainly a correct and elegant pro- 
cedure. However, if this discusser had been faced with solving 
Laplace’s equation (10) for arbitrary q"(@), his first impulse 
would have been to call on his experience with similar heat con- 
duction problems and utilize Fourier series methods directly. 
The solution would have had the same form as equations (13a) 
and (130), and the coefficients B,, and Bz, would be found by the 
Fourier theorems by integrating q’(@). It would be appreciated 
if the author would comment on the validity of this alternate ap- 
proach and also discuss the relative advantages of the two 
methods. 


Author's Closure 


Professor Sparrow’s suggestion is of particular interest to the 
author, who is currently studying this problem in turbulent 
flow. For this case it appears that the resulting eigenvalue prob- 
lem is somewhat simpler to handle for a cosine flux distribution 
than for the pulse treated in the present paper. Both methods 
are of course identical in the end result. The difference is merely 
in the way the convolution notion is employed to reduce the super- 
position integral to the form of equations (31) and (32). 
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From Momentum Model 


Theoretical equations governing slip effects in forced circulation of boiling water are 
derived. The equations indicate that steam slip is dependent upon channel geometry, 


inlet water velocity, and rate of heat addition. A simplified momentum model is 
postulated which leads to equal friction and head losses of two phases. The model gives 
good agreement with available experimental results in horizontal and vertical test sections 
with and without heat addition at pressures from 12 to 2000 psia. Discussion of the model 
in terms of nonquasi steady-state unbalances of friction and head losses of the 
two phases explains experimental deviations from the predictions and the previously 
noted effects of water inlet velocity. It also gives trends for the effects of channel geometry 
and rate of heat addition. A pplication of the simplified model to calculating two-phase 
pressure drops is included. 


Introduction 


AM AND SATURATED WATER have different proper- . 


ties and cannot be expected to flow at the same velocity in forced 
circulation boiling systems. Unequal density and viscosity of the 
two phases usually force the steam to attain a relative velocity 
or “slip” with respect to the water. Knowledge of this slip 
velocity is of utmost importance because it determines the 
density, i.e., driving head, of natural circulation boilers and void 
distribution or neutron economy of boiling water reactors. 

In recent years, considerable information has been accumulated 
about the density of two phase mixtures. Flinn and Lottes [1],* 
Cook [3], and Marchaterre [2] have measured the density of boil- 
ing water in a rectangular channel at pressures between 115 and 
600 psia. Egen, Dingee, and Chastain [4] have investigated the 
same problem at 2000 psia. Additional data have been obtained 
by other investigators for steam-water mixtures in circular pipes. 
They include the work of Schwarz [5], Behringer [6], Dengler 
[7], Schurig [8], Smith and Hoe [9], Isbin, Sher, and Eddy [10], 
and Larson [11]. Similar information has also been compiled for 
air-water and air-oil mixtures by Petrick [12], Lockhart and Marti- 


'The author wishes to acknowledge the assistance of Dr. R. O. 
Niemi and E. E. Polomik. Critical discussions with them helped 
the author formulate changes in the basic equations, 
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nelli [13], Alves [14], Johnson and Abou-Sabe [15], Zmola and 
Bailey {16], Bergelin and Gazley [17], and Baker [18]. 

Most of the work to date has, however, been of an experimental 
nature. Correlation of data has always been empirical and even 
the widely used predictions of Martinelli and Nelson [19] were 
obtained by introducing an empirical pressure dependency in the 
relations proposed in reference [13]. It is proposed herein to 
theoretically analyze the problems of steam slip in forced circula- 
tion of boiling water. The basic equations will first be derived 
and the most important variables established. Next, a simplified 
momentum model will be postulated to obtain a theoretical pre- 
diction of steam slip and two-phase pressure drop. 


Basic Equations 


Derivation of Equations. In a boiling water system with heat 
addition in the direction of flow, each fluid phase satisfies a 
momentum equation of the Bernouilli type. The equations are 
derived in the Appendix for an element length dy in the flow 
direction. They are: 


) dy — pz (sin @)dy = (1) 
LTP 


Vv 
dP + — d(AgpgVq*) + — d(A,p,V;,z) 
gAg gAg 


dP 
= (+) dy — pg(sin (2) 


Equations (1) and (2) assume that the static pressure drop for 


A = cross-sectional area of flow, ft? 
C = heat generation per unit volume of 
coolant, Btu/hr ft® 
D pipe diameter or equivalent hy- 
draulic diameter, ft flow, ft 
friction factor, nondimensional 
mass flow rate per unit area, 
lb/hr ft? 
gravitational constant, ft/hr? 
heat of vaporization, Btu/Ib 
nondimensional grouping 
pressure lb/ft? lb/ft? 
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ens 


ference 


WR 


velocity, ft/hr 

mass flow rate, lb/hr 

steam quality, nondimensional 
distance measured in direction of 


= void content, nondimensional 
nondimensional pressure drop dif- 


angle of inclination of test section 
density, lb/ft* 

shear stress at wall, lb/ft? 
frictional plus static pressure, 


exit end of heated section 

initial value corresponding to 
z=0 

steam or gas 

steam in two-phase flow 

liquid or fluid 

water alone 

water in two-phase flow 

corresponding to 100 per cent 
water flow 

two-phase mixture 
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the liquid and gas phase are equal and constant at every cross 
section and that the boiling system is inclined at an angle @ from 
the horizontal direction. 


Multiplication of Equation (1) by A, and Equation (2) by Ag 
and addition yields 


1 
dP + gA + Arp,Vz*) 


a 
(0.4 + po “2) sin 0 dy (3) 


d 

The term ( ~) represents the two phase frictional drop and 
dy /TP 

is equal to 


dy A dy /.re A dy Jarre 

Equation (3) is the basic equation proposed by Martinelli and 
Nelson to calculate steam-water pressure drop. It states that 
the total pressure drop is made up of a momentum loss (repre- 
sented by the second term on the left side of the equation) plus 
a frictional loss and a head loss (given by the last term on the 
right). Similarly, another important relation can be obtained 
from Equations (1) and (2) by subtracting one from the other. 


1 
{d[AgpgVg? + ArprV — dV 


d 
= Ag [(*) - ( + — Pg) sin a (5) 
dy dy / ure 


The advantage of writing a momentum equation for each of the 
phases is now apparent. In addition to the accepted pressure drop 
relation of Equation (3), there results a second expression, Equa- 
tion (5), which can be used to predict steam slip. 

Equation (5) can be rewritten by substituting the following: 


p:ALV, = — 2) (6) 
poAgVe = GAz (7) 


ain and 45 (8) 


A 
Equation (5) becomes for p,; assumed constant 


PL (1 — a) a Pe 2 (1 — a)? 


d 
(*) + — pa) sin a| dy (9) 
dy /ate dy / ure 


It is possible to relate the two variables y and z from a heat 
balance. Assuming thermal equilibrium at any cross section and 
a heat generation C per unit volume of coolant, one obtains 


Cdy = (10) 


In order to make Equation (9) nondimensional let us introduce 
the liquid phase pressure drop for the initial flow of 100 per cent 


water or 
dP SoG? 
— 11 
), 29Dp1, 


Equation (9) becomes 


2 
“La-a@ @ po 2 (1—a)? 
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dP dP 


(12) 


Equation (12) in its most general form is the steam slip equation. 
Discussion of the Equations. For a given heat generation, mass 
flow rate per unit area, system pressure, and hydraulic diameter, 
Equation (12) is a differential equation which relates a and z 
(i.e., defines steam slip with respect to water) if (dP/dy)arp/ 
(dP/dy)o and (dP/dy)ure/(dP/dy)» could be expressed in terms 
of z and @ and the system parameters G, D, C, and pressure. 
Even before considering such expressions, several important con- 
clusions can be drawn about steam slip from Equation (12). 

1 Steam slip depends not only upon pressure, but also upon the 
inlet velocity, hydraulic diameter, and distribution of heat 
production along the flow direction. The grouping K = Ghigfo/ 
DC can be expected to control. 

2 Steam slip will vary with flow pattern if the functions 
(dP/dy)are/(dP/dy)o and (dP/dy)ure/(dP/dy)o vary with flow 
pattern. 

3 The functions and (dP/dy)are/(dP/ 
dy)o have not been found susceptible to accurate correlation. 
Since these functions are expected to exhibit a scatter of 10 to 30 
per cent steam slip may be expected to exhibit the same order of 
scatter. 

A relation between @ and z can be derived from Equation (11) 
by two means. First, one could utilize experimental measure- 
ments of and (dP/dy)are/(dP/dy)o. Second, 
a theoretical model could be postulated to define all these func- 
tions. The first method, a priori, appears to be the most accurate 
one. However, upon close examination of test results it is found 
that the important ratios (dP/dy)ire/(dP/dy)o and (dP/dy)arr/ 
(dP/dy)»o cannot be measured experimentally and the only re- 
course is to a theoretical model. Such a theoretical model, re- 
ferred to as the “momentum exchange model” will be discussed 
next. It has the important advantage of being physically ac- 
ceptable without requiring expressions for (dP/dy)ure/(dP/dy)o 
and (dP/dy)are/(dP/dy)o. 


Momentum Exchange Model 


Definition of Model. Equation (9) can be rewritten 


@ pe 2 (l—a)? 


where y is used to represent the combined frictional and head 
losses. For instance, 


Let us first consider the simplified case of no heat addition and 
of no other source of variation of z, a, or p,;/pg. Equation (13) 


reduces to 
ay +) 
( dy a. ( dy /ure a4) 


which states that the frictional and head loss of the liquid must 
equal the frictional and head loss of the gas. This postulate was 
the basis of the semiempirical relations proposed by Lockhart 
and Martinelli [13]... Consider next small localized disturbances 
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of z, a, or pz/pg. Following these local disturbances the system 
will once again satisfy Equation (14) even though z and a@ are now 
slightly modified. If the changes occur slowly enough, one can 
assume that Equation (14) is satisfied at all times. This con- 
stitutes the basis of the momentum exchange model. In other 
words, momentum is exchanged between water and steam every 
time z, a, or p,/pg vary and the exchange tends to maintain the 
equality of frictional and head losses of the two phases. 

It is important to note that the model was tacitly utilized by 
Martinelli and Nelson when they extrapolated the unheated test 
data of reference [13] at atmospheric pressure to an evaporating 
system where z changes rapidly. The model also gives a possible 
explanation for the agreement of two-phase pressure drop data 
taken with and without heat addition in horizontal and inclined 
test sections. 

Derivation of Equations. Mathematically speaking, the model 
corresponds to writing 


d (1 — z)? z* pr, 1 (1 — 2)? 
| (15) 


Note that if Equation (14) holds Equation (15) is valid whether 
the flow is horizontal or vertical. 

Integration of Equation (15) with zero voids at zero quality 
gives 


Equation (16) can be rewritten to give 


and the quality z would vary as pg/p, instead of V po/ Pz 28 ob- 
tained from Equation (19). 

The slip ratio according to Equation (20) is not constant but 
increases with void content. This is in agreement with the pre- 
diction of Martinelli and Nelson. Machaterre’s plot of this ratio 
at 115 psia along a heated channel verifies the predicted trend of 
Equation (20). The slip ratio starts below one, increases along 
the test section, and appears to reach a constant value. Over the 
exit half of the test section the per cent change in steam volume is 
small and this partly accounts for the approximate constancy of 
the slip ratio. Equation (20) suggests a further note of caution. 
Empirical correlations of slip ratios versus inlet velocity are best 
done when the void content or quality, at which those slip ratios 
are plotted, are held constant. The effect of quality has been 
recognized and noted by Petrick. 

Comparison With Test Data. The validity of the model is best 
demonstrated by comparison with experimental data. The pre- 
diction and test results of Larson are shown in Fig. 1 for a hori- 
zontal pipe at 1000 psia without heat addition. The agreement 
is good in the region of high steam quality and void content where 
the experimental errors are small. The comparison with the 
Martinelli and Nelson prediction is also satisfactory, with the sim- 
plified model giving higher slip values particularly at low qualities. 
Similar comparisons with the data of Sher [10] and Schwarz [5] are 
given in Fig. 2 for unheated vertical and horizontal pipes at pres- 
sures of 14.7 and 1180 psia. The momentum model again yields 
lower than measured steam voids, the deviation ranging from 
about 20 to 30 per cent. 

Figs. 3, 4, 5, and 6 show the predictions and test results at 


a(l 2a) + +a} — a)? +a(l 


2a) | 


z= 
— a)? + all — 2a) 
Pa 


In the event that a is not equal to zero at z = 0 but equal to a; 
(as often reported in test results), Equation (16) becomes 


(1 — a) Pg 2 (1 — a)® 


1 
i-a’ 2a 


1 
— =0 (18) 


and it can be solved to yield a solution equivalent to Equation 
(17). 

It is important to note that Equation (17) gives a = 1 when 
xz = 1 and reduces to the accepted relation z = a@ at the critical 
pressure where p, = Pg. At low pressures where the quantity 


(1 — z) may be set equal to one even at high void contents, 
Equation (17) reduces to 


_@ Pa 
(1 — @)* py 


fad (19) 


and the ratio of liquid to gas velocity is 


Vo _ V/ 20 (20) 
Vi Pa 


The slip ratio, according to the momentum model may, there- 
fore, be expected at low pressures to vary with the square root of 
the ratio p,/pg, a dependence once proposed by Untermeyer [20]. 
On the other hand, the homogeneous model where V,, equals Vg 


would yield 
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various pressures in a vertical system with heat addition. Fig. 3 
deals with the data of Dengler and the agreement of theory and 
tests is good. Fig. 4 gives a check of the model at high pressure. 
The momentum model curve falls about 20 per cent below the 
results of Egen, Dingee, and Chastain[4].* Fig. 5 compares 
Marchaterre’s measurements at 115 psia with the theoretical anal- 
ysis. The predicted void volumes at 115 psia are about 25 to 40 
per cent too low, but the deviation decreases for the higher pres- 
sure runs of Marchaterre. The most comprehensive check of the 
model is obtained from the data of Cook. Fig. 6 shows predic- 
tions and experimental measurements when the parameter K = 
Ghigfo/CD is below 15. The agreement is good. Here again, the 
momentum model falls below the experimental measurements 
and appears to be a lower limit for the test data. 

This trend is understandable since the derivation assumed that 
Equation (15) is satisfied along the entire test section, while in 
reality departures may be expected when the steam quality varies. 
Unbalances of the steam and water pressure losses in Equation 
(15) account for the deviations as discussed in the next section. 


Momentum Model for Nonquasi Steady-State Conditions 


The preceding model assumed that the change in quality and 
fluid properties were slow enough for the water and steam pres- 
sure losses to remain equal. This is particularly true of hori- 
zontal and vertical flow in unheated channel with pressure drops 
small enough not to cause sizable changes in steam density or 
water flashing. (Note, however, a major difference between verti- 


* Assuming a different from zero at zero quality does not lower 
the deviation at high qualities between tests and prediction. The 
test results confirm that the steam voids are not highly dependent 
upon the initial void value at zero quality. 
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Fig. 1 Comparison of momentum exchange theory with Larson's data 
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cal and horizontal flow. In the vertical channel (dP/dy)arp is 
much larger rather than equal to (dP/dy)ure due to hydrostatic 
head differences. Even though the void content is the same in 
both systems at a given quality z, the flow pattern must be differ- 
ent to produce the higher gaseous frictional losses in vertical 
flow.) 

In most practical cases the rate of change of steam quality is 
rapid enough to produce different pressure losses of steam and 
water. If the steam quality were continuously changed the un- 
balance could be preserved and during this period 


(dy/dy)are — 
(dP/dy)o 


(22) 


Equation (12) can now be written 


(1 — Pe 2 (I — 
aGhifo 

B= KB (23) 

Equation (23) gives lower slip values than Equation (15) if 8 
is positive and higher ones if 8 is negative. Nonquasi steady-state 
unbalances of the two-phase pressure losses may, therefore, be ex- 
pected to give different void content than the momentum model. 
The role of the grouping K is now also apparent. Its effects upon 
slip cannot, however, be quantitatively described until an ex- 
pression for 8 becomes available. Tentative trends can still be 
obtained by examining limiting conditions of the parameter K. 
First, consider the case where heat addition is small or C — 0, 
i.e, K — ©. We assumed that as 8 approaches zero that the 
left-hand side of Equation (23) approaches zero. This was the 
single major assumption made in deriving the preceding momen- 
tum exchange model. Next consider the case of very high heat 
addition or C K 0. The term approaches a con- 
stant corresponding to the pressure drop ratio of 100 per cent 
steam flow to 100 per cent water flow. The left-hand side of 
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Equation (23) once again approaches zero. One can, therefore, 
presume that as K increases from zero to infinity the slip values 
will first agree with those of the momentum model, then deviate 
from it, and finally agree again with it at very high values of K. 

Cook’s data offer the only major opportunity of determining 
the role of K. Fig. 6 showed the experimental measurements 
where the parameter K was below 15. Similar plots are given in 
Fig. 7 for values of K between 15 and 20 and in Fig. 8 for K be- 
tween 20 and 30. The increased departure from the momentum 
exchange model is noticeable as K is raised and in particular at 
the higher values of a as predicted by Equation (23). The role of 
K appears to have leveled off in Fig. 9 where Cook’s data for K 
above 30 are shown. The decrease of slip with initial increase of 
K means that raising the inlet velocity decreased slip while in- 
creasing power density and hydraulic diameter would raise the 
slip. The trend of inlet velocity has already been noted at the 
Argonne National Laboratory but the effects of channel geometry 
and rate of heat addition have not been previously mentioned. 
Quantitative predictions of these effects must, however, await 
theoretical or experimental determinations of the function 8. 

While expressions for 8 are beyond the scope of this paper, 
further theoretical and experimental work suggests itself. Sim- 
plified theoretical models which neglect the frictional pressure 
drop of one phase with respect to the other could be utilized. 
These and other simplified models are presently being considered. 
A test program which would consist of applying heat to a steam- 
water mixture over a short distance and measuring voids before, 
during, and after heating would help to determine the function 8 
and establish the proposed momentum exchange model. Cook’s 
measurements of voids in the riser beyond the heated section were 
a first step in that direction. For those runs where the slip ratio 
appears stabilized an increase in slip was noted with the absence 
of heat production. This may be construed as a partial verifica- 
tion of the assumption that as K — = the momentum model is 
again satisfied. Unfortunately, the riser flow area was twice that 
of the heated section and the conclusion cannot be final due to 
the change in the parameter G. 
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— Pressure Drop Prediction From the Momentum 
odel 


Once the void content is specified in terms of the quality z, the 
acceleration and static losses in Equation (3) can be readily com- 
puted. The remaining term (dP/dy)rp can be obtained if the 
frictional losses are expressed in terms of the void content and 
other known system parameters. A simplified relation consists 
of writing* 

(dP/dy)rp “ 1 
(dP/dy)er ~ (1 — a)? 


(24) 


where (dP/dy)isp represents the pressure drop per unit length of 
the liquid phase if it were assumed to flow alone in the system. 


( aP ) — 2)* 
dy /.sr prD2g 
‘Equation (24) is an approximate relation because it does not 


account for the flow dependency reported in several two-phase in- 
vestigations. 


(25) 


dp 
dy/ LSP 
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A relation similar to Equation (24) was previously derived for 
an annular pattern by Levy [21] and has been proposed by Flinn 
and Lottes. It can be derived if one assumes that most of the 
liquid phase is in contact with the channel walls. The two-phase 
frictional loss for equal shear forces at steam-water interface can 


be expressed as 
(26) 


where the shear stress at the wall r,, is obtained from water flow 
conditions 

V. 2 

T. = Prfure 


2g 


If one uses the average water flow velocity and assumes that the 

wetted perimeter is confined to water-channel interface, then 
Gi — 
(1 — a) 


(27) 
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Fig. 10 Two-phase pressure drop data from Moen's thesis 
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Dure = (1 — a)D (29) 
Since the friction factor f depends only upon Reynolds number 


Sise = fire (30) 


Substitution of Equations (27), (28), (29), and (30) into (26) 
gives the relation of Equation (24). 

The validity of Equation (24) has been verified experimentally. 
A plot in reference [22] has been reproduced in Fig. 10. It is 
seen that the data of Martinelli and Lockhart and Martinelli and 
Nelson are in good agreement with it. Further, the data of 
Moen, Wickey, and Mosher at high flow rate give a further check 
of the model if the steam void contents are evaluated from the 
Martinelli and Nelson correlations. 

The two-phase pressure gradient can be obtained by integrating 
Equation (24) from start of heating where z = 0 to end of the 
heated section where z = zr, and 
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Equation (31) can be rewritten [19] 


(AP/ALye C 


= (32) 
(dP/dy)o z, 

and substitution of Equation (17) means that the ratio 
(AP/AL)re/(dP/dy)o is only a function of a and pressure, or 
again the mean void fraction & obtained from 


a 
— dz 


= - (33) 


Equations (17) and (24) were used to predict two-phase pressure 
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Fig. 11 Comparison of predicted pressure drop with data of Moen and Larson 
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Fig. 12 Comparison of predicted two-phase pressure drop with data of Isbin, Sher, and Eddy 


drop in two instances. Fig. 11 shows the two-phase pressure drop 
for the tests of Moen and Larson. Values for quality above 0.8 
were purposely not plotted as these are more readily computed 
from 100 per cent steam flow. The comparison between theory 
and tests is satisfactory and the flow effects have been somewhat 
masked. A similar comparison is shown in Fig. 12 with the data 
of Sher and the prediction for frictional plus head losses fall 50 
per cent above and 30 per cent below the experimental measure- 
ments. The deviations are partly due to the large void content 
measured by Sher and partly to scatter in the test results. The 
scatter in test results is probably caused by the effects of flow rate 
which were neglected in the analysis. 
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APPENDIX 
Derivation of Basic Equations 


The net pressure forces are obtained by considering an ele- 
mental length of fluid stream tube. The net force is aceording to 
the sketch below 


d 
aP 


dy 


PA wy) aa 
2 dy 


dP dP 
iy) (A +dA)= —-A dy (34) 


The momentum change for each fluid has been given in ref- 
erence [23]. For the liquid phase flowing at a rate W, lb/hr one 
has 


124 / may 1960 


1 
> + dW, Vz'+ dV.) — WiV, — 
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Equation (35) assumes that dW, lb/hr of fluid are converted to 
gas. The corresponding gaseous momentum change is 


1 
(We + dWeXVe + — WeVe + VidW,z) 
1 
= [d(WeVe) + VidW,) (36) 


The gravity forces are computed as follows. For the liquid 
phase they are equal to 


1 
Pi (4 +> ia) dy = pyAdy (37) 


The frictional forces are expressed by means of expression (34) 
so that for the liquid phase they are 


dP 
A({— d 
( dy 
Combination of Equations (34) to (38) yields Equations (1) and 
(2). 


(38) 
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Pressure Drop and Heat Transfer in a 
Duct With Triangular Cross Section 


Friction factors have been measured for a duct whose cross section has the shape of an 
isosceles triangle with a side ratio 5 to 1 in the fully developed flow region for laminar, 


transitional, and turbulent conditions. In addition, local and average heat-transfer co- 
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efficients and the temperature field in the duct wall have been determined for the condition 
of constant heat generation per unit volume of the duct walls. 

Friction factors in laminar flow agreed well with analytical predictions. 
turbulent flow range they were by 20 per cent lower than values calculated from relations 
for a round tube with the use of the “hydraulic diameter.” 


In the 


Heat-transfer coefficients 


averaged over the circumference of the duct were only half as large as values calculated 
from round tube relations in the Reynolds number range from 4300 to 24,000. The 
measurements also revealed that thermal starting lengths were in excess of 100 diameters. 
In round tubes a length of 10 to 20 diameters has been found sufficient to develop the 


temperature field. 


Introduction 


ECAUSE of their wide technical applications, pres- 
sure drop and heat transfer in duct flow have been studied ex- 
tensively in the past. The majority of these investigations have 
been concerned with ducts having a circular cross section. 

In many present-day cooling problems, such as in nuclear reac- 
tors, in high performance heat exchangers, or in rocket power 
plants, design considerations dictate ducts with noncircular cross 
sections. Papers in the literature which investigated the pressure 
drop and heat transfer in ducts with noncircular cross sections 
(ref. [1])* indicate that turbulent pressure drop can be calculated 
with good accuracy from relations for circular ducts when the 
diameter in these relations is replaced by the “hydraulic diame- 
ter” of the particular cross section. This rule has also been used 


1 Publication of the Heat Transfer Laboratory, University of 
Minnesota. 

2 Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division of THe American 
Society or Mecnanicat ENGIneEeRs and presented at the ASME- 
AIChE Heat Transfer Conference, Storrs, Conn., August 9-12, 1959. 

Nore: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those 
of the Society. Manuscript received at ASME Headquarters, May 
11, 1959. Paper No. 59—HT-10. 


to calculate average heat-transfer coefficients in noncircular 
ducts. It has, however, not been verified as well by experiments. 

Engineering applications which utilize cooling ducts at a high 
temperature level make it necessary to predict not only the 
average heat-transfer condition, but also the local temperature 
distributions in the duct walls. Theoretical studies [2—5] as well 
as experimental investigations [6 and 7] which have been made 
in this connection cast some doubt on the general applicability of 
the hydraulic diameter concept. The authors of this paper have, 
in addition, demonstrated [9] that for laminar flow, not only the 
local heat-transfer conditions, but also the average heat-transfer 
coefficients depend very strongly on the prescribed local tem- 
perature distribution of the duct wall in lengthwise and especially 
in peripheral direction. The analysis for a particular cross section 
demonstrated that the Nusselt number averaged over the 
periphery may differ by an order of magnitude depending on 
whether a constant wall temperature or a constant heat flux from 
the duct wall into the fluid is prescribed peripherally. 

It appeared important to investigate whether the heat-transfer 
coefficient in noncircular ducts is as sensitive to the boundary 
conditions in turbulent flow. The present paper presents the re- 
sults of an experimental study on the flow characteristics and heat 
transfer for a duct whose cross section has the shape of a narrow 
isosceles triangle with an apex angle of 11.46 degrees, correspond- 


Nomenclature 
C = perimeter of duct Pr = Prandtl number €,, = turbulent diffusivity for momen- 
A = cross-sectional area of duct de = local heat flow from inner duct tum 
b = duct wall thickness wall surface to air stream €, = turbulent diffusivity for heat 
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L ¢, = specific heat of sir at constant qe = heat generated electrically per 7 = Cartesian co-ordinate in plane of 
pressure unit time and duct surface area duct cross section 
minus heat loss to outside = dynamic viscosity 
Re = Reynolds number, see Equation v = kinematic viscosity 
J = friction factor, see Equation (3) (10) p = density of fluid 
h = convective heat-transfer co- 7 = air temperature T. = wall shear stress 
efficient T, = air bulk temperature & = Cartesian co-ordinate in plane of 
k = heat conductivity T.. = duct wall temperature duct crose section 
k, = heat conductivity of duct wall velocity in z-direction Indexes 
L = distance between pressure taps ae soartinate along cut wall im * = reference state 
circumferential direction — = average over circumference or 
m = mass flow through duct 2 = co-ordinate in direction of duct cross section of duct 
Nu = Nusselt number, see Equation (8) axis B = bulk value 
p = static pressure a = thermal diffusivity w = duct wall 
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ing to a side ratio 5 to 1. This shape was chosen because it was 
expected that with such a strong asymmetry of the cross section 
and with the sharp corner, conditions would differ in a pro- 
nounced way from those in a circular tube. In addition, coolant 
passages with triangular cross sections are used in various of the 
applications that have been mentioned. 


Experimental Apparatus 


The experimental apparatus consisted of a duct with triangular 
cross section, partially wood and partially metal, with an over-all 
length of 209 hydraulic diameters. As illustrated in Fig. 1, under 
the influence of a downstream blower, air flows through an as- 
sembly (a) consisting of fine mesh screens to remove disturbances 
from the stream and enters the duct through a rounded entrance. 


Fig. 1 Diagram of experimental flow circuit 

Next, a wooden duct section (b) 93 hydraulic diameters in length, 
allows the flow to become fully developed hydrodynamically before 
entering the heated section. This wooden section was removed 
in some experiments. The heat-transfer section (c) which is made 
from a stainless-steel plate, 0.039 in. thick, is 116 hydraulic 
diameters long, has 119 thermocouples to measure the duct wall 
temperatures, and is surrounded by a wooden box filled with 
“Santocel A” insulation. (Thermal conductivity = 0.0142 
Btu/hr ft F at a temperature of 150 F.) The box runs the whole 
length of the metal duct section with overlaps at both ends. The 
minimum thickness of the insulation was 5 in. 

Downstream from the duct are a mixing chamber (d) for meas- 
uring the exit bulk temperature of the air with thermocouples 
(e), flow measuring instruments (either rotameters or orifices) and 
the blower which returns the heated air to the room. Located at 
the ends of the wooden and metal sections are pressure taps (f) 
to measure the static pressure level at those respective cross sec- 
tions. 


Fig. 2 Construction details (c) and final dimensions (b) of heat-transfer 
duct 


The metal heat-transfer duct was formed from a single sheet of 
stainless steel 0.039 in. thick. The construction details and final 
dimensions of the metal duct are shown in Fig. 2. Grooves were 
milled along lines representing the corners which were then bent 
to the correct angle. The grooves served to simplify bending and 
eliminate excess corner material. The single seam, which occurs 
in the center of the triangle base, is a tongue and groove connec- 
tion held together with glass fiber tape painted with “‘glyptol,”’ 
and tested for leaks. The hydraulic diameter of the duct is 0.904 
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Fig. 3 Thermocouple locations in heat-transfer duct 


in. and the apex angle is 11.46 deg. The duct was arranged in 
such a way that the height of the triangular cross section and the 
axis were in a horizontal plane. 

As mentioned previously, 119 thermocouples were installed in 
the metal duct at locations shown in Fig. 3. For convenience in 
data representation, the locations of the measuring stations are 
given in the figure to the nearest hydraulic diameter. The ther- 
mocouple installation was performed by drilling a small hole 
slightly larger than the thermocouple junction part way through 
the duct wall, filling the hole with copper oxide cement, and then 
carefully holding the junction in place away from contact with 
the duct until the cement hardened. The junctions were thus 
not in electrical contact with the duct, avoiding spurious voltages, 
and were located approximately 0.010 in. from the inner duct 
wall. 

The metal test section was heated under the condition of con-. 
stant heat generation per unit wall volume by passing an alternat- 
ing current through the stainless-steel wall. Power was supplied 
and controlled by means of a transformer and ‘“Variac’’ combina- 
tion with the current entering and leaving the duct through a 
large number of fine wires attached at either end around the 
periphery of the duct and passing the electric current in parallel. 
The purpose of these wires was to insure a constant electrical po- 
tential around the end cross sections and simultaneously to pre- 
vent or at least reduce conduction losses from the duct ends. 


Accuracy of Measurements 


The important measurements in these experiments included 
wall temperatures, entering and leaving air temperatures, air 
mass flows, static pressures, and heat-generation rates. In the 
following paragraphs the accuracy of these measurements is 
examined. 

Thirty-gage iron-constantan thermocouples were used to make 
wall temperature measurements. The following calibration pro- 
cedure was followed. In order to avoid the necessity of individually 
calibrating 119 thermocouples, special thermocouples of uniform 
length were made from each of five rolls of iron-constantan wires. 
These thermocouples were then calibrated against mercury-in- 
glass thermometers which had been certified by the National 
Bureau of Standards to be accurate to 0.05 F. The thermocouple 
emf’s were read with a Brown self balancing potentiometer whose 
error is believed to be less than the equivalent of 0.1 F. When all 
calibrations were compared, it was found that the maximum dif- 
ference between thermocouples was 0.3 F. This latter figure 
was then taken as the uncertainty in the temperature measure- 
ments. 
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The temperatures measured in the duct wall are not the inner 
surface values, because the thermocouple junctions were posi- 
tioned within the stainless-steel duct wall. A calculation in- 
dicated that even under the most unfavorable conditions this dif- 
ference between the measured and the surface temperature should 
be less than 0.3 F or the calibration uncertainty and thus no 
corrections were made for this effect. Every effort was made to 
minimize conduction errors by running the thermocouple wires 
along the outer surface of the duct in flow direction (where the 
temperature gradients are least) before taking them to the 
potentiometer. Under these conditions it is felt that 0.3 F is a 
reasonable uncertainty figure for measurement of the wall inner 
surface temperature. 

The temperature of the air entering the heated section was 
measured by a thermocouple two inches upstream from the end of 
the wooden section. The thermocouple was shielded from any 
radiation coming from the metal section and since the air tem- 
perature measured in this way always agreed with the room 
temperature measured in the vicinity of the entrance with a mer- 
cury thermometer, it is felt that these data are reliable. 

The temperature of the air leaving the heated section was meas- 
ured by a seven-junction averaging thermocouple after the air 
had been thoroughly mixed by a series of baffles. Radiation and 
conduction errors in these measurements were shown by calculs- 
tion to be negligible. 

Mass flows of the air were measured by either an orifice, built to 
standard ASME specifications, or a rotameter depending upon the 
flow range. In the region of overlap between the two methods, 
the measurements agreed to within one per cent. 

Static pressures were measured through wall pressure taps 
having a diameter of 0.020 in. with no inside ridges to disturb 
the flow. Pressure differences were determined using a micro- 
manometer (Meriam Type A-750) which can be read within 0.002 
in. of water. In the worst case, when laminar pressure drops were 
being measured, the accuracy was estimated to be within two 
per cent. 

The heat-generation rate was determined by reading the total 
power dissipated in the duct. The current through the duct was 
measured by a calibrated laboratory ammeter (Westinghouse 
Type PY 5) and a current transformer. The meter accuracy was 
0.25 per cent of full scale or 0.0125 ampere. The electric potential 
difference across the duct was measured by a Weston model 904 
voltmeter with a guaranteed accuracy of '/: per cent of full scale 
of 0.015 volts. In these experiments, the uncertainty in the 
power input varied between one and two per cent. 

The amount of heat transferred to the air could be determined 
from the mass flow of the air and its entering and leaving tem- 
perature. The difference between the electric power input and 
the energy absorbed by the air was between 2 and 10 per cent, 
This difference is interpreted as external heat loss of the apparatus 
and is listed in Table 1. The losses determined in this way were 
later applied to the data as correction factors. Relaxation calcu- 
lations of the heat loss through the insulation surrounding the 
duct had previously established the validity of the correction 
procedure. 

An investigation was also made of the possibility of unequal 
rates of heat generation in various sections of the duct due to the 
change in duct temperature and to the temperature dependence 
of the electric resistance of stainless steel. A rather lengthy 
numerical calculation showed this effect to be negligible. 


Dimensionless Parameters and Their Determination 


The results of the present investigation are appropriately pre- 
sented in the form of dimensionless parameters. The type and 
number of the parameters, which describe heat transfer under the 
conditions as outlined before, will be determined by a dimensional 
analysis. 
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The flow field for fully established turbulent flow is described by 
the following differential equation 


ow ow 1 dp 
+2 lore (1) 


The symbols appearing in this and the following equations are 
contained in the Nomenclature. The fluid properties related to 
the present investigation will be assumed as constants. If Equa- 
tion (1) is made dimensionless in the customary way (see for in- 
stance ref. [10]), then it can be anticipated that the solution to 
this equation must have the following form 


w 9 
(Re, 2) (2) 


In developed flow it can be assumed that the turbulence level is 
determined by the internal flow process and that it can not be in- 
fluenced by changes in the inlet condition. As long as the duct 
walls are hydraulically smooth, the ratio ¢,,/v can therefore not 
be regarded as a free parameter. The solution of Equation (1) 
also leads to a relation between the pressure gradient and the 
mean velocity. 

The pressure drop in the duct may be described by a triction 
factor defined by the following equation 


2 
Sp (4) (3) 
L pw? L m 

The temperature field in the duct wall can be described by the 
following differential equation, assuming that temperature varia- 


tions through the duct wall normal to its surface are negligibly 


small. 
2 /, or 


The left-hand term in Equation (4) describes the net heat flux into 
a volume element of the wall by heat conduction in the peripheral 
direction z. A corresponding term for the direction z of the duct 
axis was omitted because it was found to be small as compared 
with the contribution of conduction in z-direction. 

The temperature field in the fluid flowing through the duct is 
described by the equation 


re) oT oT oT 
| +2 fata (5) 


Heat transport by conduction and turbulent mixing is again 
neglected in z-direction. The form of the solution to the two 
equations (4) and (5), which describes the duct wall and the fluid 
temperatures can again be anticipated by making these equations 
dimensionless. This process indicates that a relation of the 
following form must exist for the wall temperature: 


kb s 
(i. x) @ 


The Reynolds number appears in this equation because the rela- 
tion (2) for the flow velocity w has to be introduced into Equation 
(5). The wall temperature is indicated in this equation as a 
parameter AT, because Equations (4) and (5) are obviously in- 
dependent of temperature level as long as the properties are con- 
stant. The temperature can then be measured from any reference 
value and the fluid bulk temperature is used for this purpose in 
the following presentations. 

A heat-transfer coefficient A averaged over the periphery of a 
particular cross section will be defined in the following way 


Ge = KT, = Ts) (7) 
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The average Nusselt number, as a dimensionless expression for 
this heat-transfer coefficient, is then determined by an equation 
of the following form 


AD, ( kb 


Nu ED,’ 


2 

Re, Pr, (8) 
The parameter on the left side of Equation (6) is therefore equal 
to 1/Nu when AT, is interpreted as T,, — T'y. The relation (8) 
is presented later on. In addition, the dimensionless parameter 
(T, — — a8 representation of the local wall tem- 
perature will be shown as well as a local heat-transfer coefficient 
h which is defined by the equation 


= A(T. T's) (9) 


The parameters appearing in these relations are obtained from 
the measured values in the following way. The Reynolds number 
is defined as 

mD, 

Ap 
The properties are generally introduced into the dimensionless 
parameters at the local fluid bulk temperature. The average wall 
temperature is obtained by an integration of the measured 
wall temperatures around the periphery of the duct. The fluid 
bulk temperature rises linearly from the beginning to the end of 
the heated duct provided the heat addition to the fluid per unit 
length and the specific heat of the air is constant. The local bulk 
temperature was determined in this way from the measured en- 
trance and exit temperatures of the heated section. The average 
Nusselt number then was obtained from Equations (7) and (8). 
The heat conductivity of the duct wall material was taken from 
the literature as 9.27 Btu/hr ft deg F. This value, the length 
dimensions of the duct, and the heat conductivity of the air de- 
termine the parameter (k,b)/(kD,). Jt varied in the present 
experiments between 23.2 and 24.9. This variation is considered 
small and the value 24 has been assigned to the results presented 
later on in this paper. The Prandtl number of air varied in the 
experiments between 0.709 and 0.700. 

The calculation of the local Nusselt number is somewhat more 
involved. A heat balance was made for this purpose on an ele- 


Re (10) 


1 


Fig. 4 Control volume for calculating local heat-transfer coefficients 


ment with 1 in. X 1 in. dimensions located in the wall as indi- 
cated in Fig. 4. The heat balance on the wall element can then 
be approximated by the following equation when the heat con- 
duction in z-direction can be neglected. 


+ (Ta + Twa 2T w) (11) 
2x0 


h= 
Ts 
The temperatures 7’, T'1, and 72 are obtained from the meas- 
ured temperature profiles along the duct periphery; the measure- 
ment of the other quantities appearing in this equation has al- 
ready been described. 


Results of Pressure Drop Measurements 


The first series of measurements are concerned with the static 
pressure drop for the fully developed isothermal flow through 
the metal heat-transfer duct. The purpose was to obtain a partial 
check on the duct dimensions, sealing integrity, and flow transi- 
tion by comparison of the measurements with analytical laminar 
flow solutions and with previous measurements [7 and 8]. The 
pressure drop information which is conveniently presented as 
friction factor against Reynolds number with the hydraulic 
diameter taken as the characteristics length dimension is shown 
in Fig. 5. 
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Fig. 5 Fully developed friction factors in 11.46 degree duct. Solid points calculated by method of ref. [5]. 
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Table 1 


Total heating 
Total bulk power input Average 


Run Mass flow rate, Reynolds number temp rise, to duct, losses, 
no. lbm/sec Entrance Exit deg F Btu/hr per cent 
1 0.02097 7440 6840 62.5 1222 5.2 
2 0.01326 4725 4340 65.8 795 5.3 
3 0.03005 10580 9640 69.9 1817 2.0 
4 0.04381 15620 14380 61.2 2472 6.5 
5 0.04500 16000 14740 59.9 2463 5.5 
6 0.0618 21750 20330 51.3 2940 7.0 
7 0.0682 24100 22570 48.7 3070 6.7 
8 0.0379 13420 12510 58.3 1990 4.2 
9 0.01806 6400 5950 53.1 904 8.2 
10 0.0143 4990 4615 59.0 810 9.9 
ll 0.0210 7360 6830 56.1 1093 6.5 
12 0.02864 10160 9420 56.8 1054 6.6 
13 0.0450 15850 14820 50.2 2107 7.3 
14 0.0608 21660 20340 46.7 2650 7.3 
15 0.0673 23800 22200 51.2 3210 7.3 


In the laminar flow regime, the experimental friction factors 
are compared with the calculated friction factors for a circular 
sector duct with the same opening angle. Below a Reynolds 
number of 1000 the agreement appears to be satisfactory. It 
can be seen that at Reynolds numbers larger than 1000, the 
data gradually depart from the laminar solution. This is in agree- 
ment with previous measurements [7 and 8], in which it was 
found that the transition process in such a duct is a gradual one 
with the turbulent flow appearing first in the duct core near the 
triangle base and filling more and more of the duct as the Reynolds 
number increases. 

When the flow becomes fully turbulent, it is seen that the ex- 
perimental data fall approximately 20 per cent below the Blasius 
expression for a circular tube. This does not agree with either 
our previous experience on a trianguiar duct with 23-deg apex 
angle [8], or with the measurements of others on ducts with 
various cross sections [1] where it has been found that pressure 
drop can be represented by the circular duct relation if the hy- 
draulic diameter is used to define the friction factor and Reynolds 
number. A search of the literature, however, failed to reveal any 
measurements on a triangular duct with such a small apex angle. 

In an attempt to account for the departure of the experimental 
data from the circular tube solution, a calculation of the velocity 
field was made using the method of Deissler and Taylor [4 and 5}. 
They assume that on lines normal to the lines of constant velocity, 
the turbulent velocity distribution may be represented by the 
universal velocity profiles found to describe the flow in circular 
ducts. With this assumption, it is possible by an iterative proce- 
dure to calculate the velocity field and thus the friction factor in a 
duct such as the present one. The calculation was made for two 
Reynolds numbers and the results are shown as solid points in 
Fig. 5. The agreement between the measured and calculated 
values is very good. It was, however, found less satisfactory 
when the calculation was repeated for a triangular duct with 23- 
apex angle. The calculated friction factors were found to be only 
slightly higher than the values for the present duct, whereas the 
measured values reported in ref. [8] agreed well with the circular 
duct line in Fig. 5. 

On the basis of this figure it can be assumed that turbulent 
friction factors for triangular ducts with small apex angle are 
smaller than the values calculated by the hydraulic diameter con- 
cept and can be calculated with better accuracy by the method of 
Deissler and Taylor. An experimental investigation of pressure 
drop in even more extreme cross sections is the subject of a cur- 
rent program in our laboratory. 

Summarizing the pressure drop measurements, it may be said 
that the laminar flow pressure drop was found to be in good 
agreement with calculations and that the transition process was 
the same as previously observed. On the other hand, the turbu- 
lent pressure drop was approximately 20 per cent lower than pre- 
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dicted by Blasius’ circular tube relation. Under these circum- 
stances it would not be surprising if the fully developed turbulent 
heat transfer also deviated from the circular duct correlation. 
That this is the case, and that the deviation is many times 
greater, will be shown in the next section. 


Results of Heat- Transfer Measurements 


A series of runs was made, at a variety of duct Reynolds num- 
bers, in which the wall-temperature distribution and the appro- 
priate flow parameters were measured. In order to perform these 
measurements, the heating current was applied to the duct, the 
desired air flow was set, and the equipment was allowed to operate 
until a steady-state condition was reached. The existence of 
steady state was determined by the lack of change of the duct wall 
temperatures with time. Using the procedure given, a total of 15 
runs was made. A summary of important measurements and 
operating parameters is given in Table 1. Reference to this table 
indicates that all measurements were taken in a Reynolds number 
range between 4300 and 24,000. It can be seen from Fig. 5, that 
the flow is already turbulent at the minimum Reynolds number. 
Therefore the data discussed in this section will be for turbulent 
conditions only. 

It was originally intended to obtain data under conditions of 
laminar flow so that some comparisons could be made with the 
analytical solutions presented in ref. [9]. The experiments at low 
Reynolds numbers, however, indicated the existence of free con- 
vection flows which precluded such comparisons. These second- 
ary flows were perceived indirectly by noting a difference between 
the temperature distribution on the top long side of the duct and 
on the bottom side. It was only when the flow became turbulent 
that these temperature differences disappeared and the tempera- 
ture distribution along the two sides of the triangle became sym- 
metrical. This experience indicates that caution should be exer- 
cised in applying laminar flow heat-transfer solutions for non- 
circular ducts to actual, situations. The secondary currents en- 
gendered by the wall-temperature differences at a given cross 
section may produce free convection effects which are of the same 
order of magnitude as the forced convection heat transfer. 

The measured wall temperature distributions for three 
Reynolds numbers are shown in Fig. 6(a, 6, and c). In this 
figure, the difference between the local wall and the cross-section 
bulk temperature is plotted as a function of the dimensionless dis- 
tance from the corner. The Reynolds numbers shown are the 
average between the values at duct inlet and exit. When the 
temperatures are plotted in such a way, the existence of a fully 
developed temperature field should be indicated by the coinci- 
dence of the temperature profiles for all sections in the developed 
region. In Fig. 6(a), it will be noted that the temperature dis- 
tributions for Sections 7, 8, and 9 are nearly coincident. It 


may 1960 / 129 


ia 
a 
7 
= 
a 
| 
| 
| 
i 
4 


i = Re* 4530 
aN 
30 
020 040 060 080 100 
6(a) 


130 / may 1960 


Fig. 6(b) 


Fig. 6 Transverse wall-temperature distributions along heat-transfer duct 


may therefore be inferred that the temperature field is fully 
established past the location of Section 7. Reference to the table 
in the figure indicates that a total of 71 hydraulic diameters were 

in order to fully establish the turbulent temperature 
field. Fig. 6(b and c), however, shows that with increasing 
Reynolds number the temperature field does not become fully de- 
veloped in the investigated duct length because a significant 
change in the temperature profiles persists in all sections to the end 
of the duct. The excessively long thermal starting lengths indi- 
cated by Fig. 6 are somewhat surprising, since measurements as 
well as calculations [11 and 12] predict that, for turbulent flow in 
a circular tube, the thermal entrance length is generally not in 
excess of 20 diameters. 

The duct wall-temperature field was replotted in a manner 
shown in Fig. 7(a, b, and c). In this figure, the difference be- 
tween the local wall and the fluid bulk temperature in the same 
cross section is plotted against the distance down the duct made 
dimensionless with the hydraulic diameter. The individual curves 
are for a constant wall location away from the duct corner. This 
figure again is plotted for three average Reynolds numbers and 
corroborates the information presented in Fig. 6. For the 
Reynolds number 4530, the temperature profiles, except in the 
immediate vicinity of the narrow corner, become horizontal in- 
dicating that the field is developed. At the higher Reynolds 
numbers, however, the profiles never attain a constant value, im- 
plying that the thermal field does not become established over 
the entire duct length of 116 diameters. 

The region of a fully developed thermal field can also be defined 
as that region in which the heat-transfer coefficient is independent 
of axial location. This region can be determined from Fig. 8 in 
which the average Nusselt number, as defined by Equation (7), 
is plotted over the dimensionless distance z/D,. Fig. 8 illustrates 
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Fig. 6(c) 
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Fig. 9 Average cross section Nusselt numbers at various axial stations 
(hydrodynamic starting section removed) 


the same phenomenon as discussed previously. At the lower 
Reynolds numbers, the average Nusselt number drops from large 
values at the start of the heated section and approaches a con- 
stant value at a duct location in excess of 70 diameters. At 
Reynolds numbers in excess of 6000, the Nusselt number never 
becomes constant but is always decreasing over the range of 
z/D, investigated. 

Since these measured values of the thermal starting length in 
the triangular duct are so at variance with the circular duct in- 
formation, every effort was made to make sure that this behavior 
was pot due to some peculiarities of the particular setup. Such 
factors as duct leakage and duct distortion by the interior pres- 
sure at higher Reynolds numbers were investigated but no fact 
was found that could contribute to the observed behavior. 
Physical property changes which accompany the increase in 
temperature level as the fluid travels down the duct are too 
small to cause departure from the constant property analysis 
which predicts an asymptotic value of the heat-transfer co- 
efficient. It therefore appears that the extended starting length 
is a peculiarity of this duct geometry with the applied boundary 
conditions. A possible explanation will be presented later on. 

A series of runs was also made for which the wooden entrance 
section was removed and thus the temperature and velocity fields 
were required to develop simultaneously in the heated test section. 
Except in the region near the entrance, little difference was noted 
in the performance of the duct whether or not the velocity field 
was established at the beginning of the heated section. Fig. 9 
shows the longitudinal variation of the average Nusselt number 
with this sliarp inlet. Cross reference to Fig. 8 will show that the 
performance is essentially the same. 

The fact that the temperature field was not established within 
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Fig. 10 Average cross section Nusselt numbers at Section 9. Properties 
evaluated at local bulk temperature. 


the duct length which was used for the experiments causes some 
difficulties. It means that strictly the Nusselt number obtained 
in the present tests should not be compared with information for 
thermally established flow in a round tube. Nevertheless, such 
a comparison is presented in Fig. 10 for the cross section 9 which is 
farthest away from the start of the heated section. Values ob- 
tained with the hydraulic entrance section and with a rounded 
inlet are indicated by circles and values for simultaneous hy- 
draulic and thermal development by the triangles. It may be 
observed that the difference between both flow situations has a 
negligible effect on heat transfer at the distance of 114 hydraulic 
diameters. The experimental points lie well on a straight line 
which is represented by the equation 

Nu = 0.0325 Re (12) 
The physical properties were introduced into the dimensionless 
parameters at the local fluid bulk temperature. An evaluation 
was also carried out in which fluid properties were evaluated at 
a reference temperature halfway between the fluid bulk tempera- 
ture and the average wall temperature at the particular duct 
cross section were used. The experimental points evaluated in 
this way arranged on a straight line represented by the equation 


Nu* = 0.0318 Re**. (13) 


It may therefore be observed that the difference between both 
evaluations is quite small. For this reason the bulk temperature 
was used as reference temperature in the rest of the evaluations 
presented in this paper. This has the additional advantage that 
the bulk temperature is directly established by the measurements. 
Fig. 10 also contains the Dittus-Boelter relation for established 
turbulent heat transfer in a round pipe. It can be observed that 
Nusselt numbers obtained in the triangular duct, when based on 
the hydraulic diameter, are only approximately half as large as the 
Nusselt numbers for turbulent pipe flow. They also increase more 
slowly with Reynolds number than the pipe Nusselt numbers. In 
this connection it should be remembered that the temperature field 


in cross section 9 is further away from established state the larger 
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Fig. 11 Dimensionless wall-temperature distribution at Section 9 com- 
pared with laminar solution from reference (9) 
the Reynolds number. Therefore it is to be expected that, for 
completely established thermal situation, the increase of Nusselt 
number with Reynolds number will be even smaller. ~ 
Fig. 10 also contains a solid point which has been calculated by 
a method which is essentially a modification of Deissler and 
Taylor’s procedure [5] and which is described in more detail in 
the Appendix. This point was determined for the condition that 
the heat flux to the fluid is constant around the duct periphery. 
The present experiments have a peripheral boundary condition 
for the temperature field of the fluid which lies between a constant 
heat flux and a constant wall temperature. For this situation 
one has to expect a Nusselt number larger than the value shown 
in Fig. 10 because of the influence of wall conduction. The 
measured Nusselt numbers are, however, very much lower than 
this calculated value. An explanation for this fact is suggested 
by an inspection of Fig. 11 in which the difference between the 
local wall temperature and the fluid bulk temperature is plotted 
over the distance from the duct corner. The figure indicates that 
the temperature variation around the duct wall is up to two and 
one half times as large as the difference between the average 
wall temperature and the bulk temperature. This means that in 
the corner region the heat entering the fluid from the duct wall 
has to flow in a direction toward the base of the triangle before 
it is picked up and carried away by convection in the bulk of the 
fluid. It may well be that the turbulent diffusivity describing 
the turbulent heat exchange in this direction is smaller than the 
turbulent diffusivity in a direction normal to the duct walls. 
This directional characteristic of the turbulent diffusivity may be 
connected with the fact that the velocity gradient in the direction 
from the sharp corner toward the base is much smaller than the 
velocity gradient normal to the duct walls [7]. It is also suggested 
by Prandtl’s mixing length theory, according to which the turbu- 


d 
lent diffusivity ¢,, is = with / indicating the mixing 
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Fig. 12 Distribution of dimensionless local heat-transfer coefficient at 
Section 9 compared with isothermal calculation (dashed line) by method 
of reference [5|, Re = 14,300 


length. The diffusivity depends, according to this equation, on 
the velocity gradient. The assumption of a small turbulent dif- 
fusivity also explains the fact, discussed before, that an unusual 
duct length is required to establish the temperature field. One 
might attempt to account for the starting length by referring to 
the height of the duct rather than to the hydraulic diameter. 
The corresponding ratio is 5.43, and on this basis the length of our 
duct is 21.5 triangle heights. 

The measured temperature variations around the periphery of 
the duct wall are compared in Fig. 11 with a calculated tempera- 
ture variation for laminar flow. The calculation, however, is 
made for a peripherally constant heat flux q,, into the fluid which is 
produced in a duct with uniform heat generation when the parame- 
ter (k,b)/(kD,) is equal to zero. Since the corresponding 
parameter for the experiments has the finite value 24, a direct 
comparison is impossible. In a qualitative way, however, the 
two curves have a similar shape. The drop of the measured wall 
temperatures near the sharp corner is probably due to increased 
heat losses on the outer surface of the duct in this region. 

Fig. 12 presents the local heat-transfer coefficients for the same 
duct cross section as obtained by evaluation of the measurements. 
Theory indicates that the local heat-transfer coefficient should 
drop to the value zero in the duct corners. The fact that finite 
values have been obtained from the measurements is partially ex- 
plained by the finite differences used in Equation (11), partially 
by the larger external heat losses from the corners, partially by 
radiative heat exchange from the corners toward the cooler parts 
of the duct wall. For comparison purposes, turbulent local heat- 
transfer coefficients are plotted as a dashed line as they have 
been calculated by the method of ref. [5] with the boundary con- 
dition that the duct wall temperature is peripherally constant. 


Conclusions and Recommendations 
Experiments have been performed to establish friction factors 


may 1960 / 133 


38 
| 
| 
| 
4 |\ — = 
N 
0. 
Ate 
06 
| 
020}- + t— / © 22,600 02 
| -020 } 0 02 o4 06 
“040 
idka 
. 
ix 
Be 
j 


and heat-transfer coefficients for flow of air through a duct with 
the cross section of an isosceles triangle having a side ratio of 5 
to 1. Friction factors were measured for established laminar 
and turbulent flow. The laminar flow data agreed well with 
calculated values. The turbulent friction factors were approxi- 
mately 20 per cent lower than friction factors obtained from 
circular duct information with the hydraulic diameter concept. 
They agreed, however, well with values calculated by Deissler 
and Taylor’s method. 

Heat-transfer coefficients and temperature distributions in the 
duct walls could be obtained for turbulent flow only. The duct 
wall was heated electrically and the distribution of the heat flux 
around the periphery of the duct was determined by the uniform 
heat generation and by conduction within the duct walls. The 
parameter (k,b)/(kD,) describes this effect of heat conduction 
and it was established by analysis of heat transfer in laminar 
duct flow that it influences local and average heat transfer con- 
siderably. The present experiments were made for a fixed value 
of this parameter equal to 24. The measurements indicated that 
a considerable variation of the wall temperature existed in all 
cross sections of the duct. This temperature variation was from 
two to three times as large as the difference between the average 
wall temperature and the fluid bulk temperature at any cross 
section. An established temperature field existed in the down- 
stream portion of the duct at a Reynolds number 4530. For 
larger Reynolds numbers, duct lengths beyond z/D, = 114 are 
required to obtain an established temperature field. The average 
Nusselt number for the end cross section of the duct which comes 
closest to an established situation was found to be only half as 
large as Nusselt numbers calculated with the hydraulic diameter 
concept from circular duct relations. It also increases at a 
slower rate with Reynolds number than the heat transfer in 
circular ducts. It is suspected that the turbulent diffusivity on 
which the heat transport primarily depends is considerably 
smaller in the present duct than in a round pipe. 

The results mentioned indicate that the maximum wall tem- 
peratures of similar duct geometries will be in any engineering 
application considerably larger than the values obtained from 
circular duct information with the hydraulic diameter concept. 
Additional research is therefore obviously required to clarify the 
energy transport mechanism and to obtain information from 
which wall temperature variations and maximum temperatures in 
heat exchangers can be predicted. 
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APPENDIX 


R. Deissler and M. Taylor calculated in ref. [5] the velocity 
field in developed turbulent flow through a noncircular passage 
on the basis of two assumptions which will be discussed in con- 
nection with Fig. 13. Assume that the velocity field is known and 
described by lines of constant velocity as illustrated in Fig. 13(b). 


Fig. 13 Sketch of universal velocity and temperature profiles (c), 
velocity field (b), and constant peripheral heat rate temperature field (c) 


Then attention is concentrated on the hatched small area AA’, 
with the base Az at the duct wall and bounded by two lines which 
are st any point normal to the lines of constant velocity. It is 
then assumed that no net flux of momentum occurs into the 
fluid moving in axial direction within the space appearing as area 
AA’ in the plane of the drawing from the rest of the fluid. With 
this assumption, the following force balance exists 


= AA’ (14) 
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It is also assumed that the velocity profile along the dashed center 
line of AA’ is identical in shape with the universal profile describ- 
ing the velocities in a circular duct. In addition, it has, of course, 
to be required that proceeding from any location at the wall along 
a line normal to the curves of constant velocity one ends up with 
the same maximum velocity at the center point O’. The relation 
between the turbulent wall shear and the universal velocity pro- 
file in a circular duct gives thus a second possibility to calculate 
the wall shear stress in the noncircular duct at the location Az. 
Lines of constant velocity in the duct have to be established by 
trial and error in such a way that both relations for the wall shear 
are simultaneously fulfilled. 

To actually carry out this analysis, a diagram for the universal 
velocity profile in a circular tube is needed. Such a diagram is 
indicated in Fig. 13(a). It is presented as a dimensionless velocity 


(15) 
plotted over a dimensionless wail distance 


(16) 


yt 


For our later calculations, a universal temperature profile will be 
needed which is also indicated in Fig. 13(a) as a dimensionless 
temperature 


— T) 
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(17) 


plotted over the dimensionless wall distance. Such a curve exists 
for each Prandt! number provided the properties of the fluid can 
be assumed constant. Deissler and Taylor introduce, in addition, 
the following dimensionless parameters 


(18) 


(19) 


They extended their method to a calculation of the temperature 
field and heat transfer in a noncircular duct for established turbu- 
lent flow under the condition that the temperature gradient in the 
fluid in axial (z) direction is constant and that the wall tempera- 
ture is constant around the periphery or varies by an amount 
which is small compared to the difference between wall tempera- 
ture and bulk temperature. An inspection of Fig. 11 will show 
that this condition does not apply to our duct geometry. 

For this reason we extended the method for the condition that 
the heat flux from the wall into the fluid is constant around the 
periphery and that the wall temperature differences around 
the periphery are large. The procedure is again based on the as- 
sumption that no energy transport occurs in directions which lie 
on the surfaces of constant temperature and that the temperature 
varies along lines which are everywhere norinal to the surfaces of 
constant temperature in the way described by the universal tem- 
perature profile for a duct with circular cross section. Fig. 13(c) 
indicates, within the triangular duct cross section, lines of con- 
stant temperature. Ap area element AA may be bounded by 
the length Az along the duct wall and by the two curved lines 
which are everywhere normal to the isotherms. This area ele- 
ment, hatched in Fig. 13(c), has the dashed center line. The 
following energy balance holds for a cylindrical control volume 
with the bases AA and with height dz: 
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Azq, = Wm AA (20) 


The term w,, in this equation denotes the average velocity through 
the area element AA. A similar relation holds for the whole duct 
cross section A 


dT 


= pc, wA (21) 


From both relations (20) and (21) one obtains 
D, Az 

— —— grt 

4 AA 


w,** = 


(22) 


In this equation the velocities w have been additionally replaced 
by the dimensionless parameters defined in Equation (18). Equa- 
tion (22) gives a relationship by which an original assumption of 
the temperature field can be improved in an iteration process. 

The procedure which we used is as follows: One assumes the 
location of point O within the cross section at which the tem- 
perature has a maximum or minimum. Then a number of iso- 
therms are drawn as a first estimate and lines are inserted which 
start at the origin O and are everywhere normal to these iso- 
therms. In this way the area elements AA are established. 
For any of these area elements, for instance, for the one hawhed 
in Fig. 13(c), the length yo can be measured. With this value 
and with Equation (16), the parameter yo* is calculated and the 
temperature parameter 7+ is read from the universal tempera- 
ture profile in Fig. 13(a). The following parameter is used to 
describe the temperature field 


(23) 


It is connected with the parameter 7'* by the equation 
To* — T* 
AA 
Az D, 


= 


(24) 


This equation permits the calculation of the parameter ++ 
when the term 7'* is made equal to 0. In order to establish a 
new, improved isothermal curve, a value @++ is arbitrarily as- 
sumed. The corresponding temperature parameter 7'* is cal- 
culated from Equation (24), the corresponding value y* is read 
off Fig. 13(@), and the distance y is calculated from Equation (16). 
It is assumed in this connection that the velocity profile and the 
local wall shear stress 7,, have already been established by a 
previous calculation. A point of the isothermal curve is in this 
way obtained by plotting the length y along the dashed line in 
Fig. 13(c). A repetition of this procedure in other area elements 
AA results in other points of the same isotherm so that this iso- 
therm can now be inserted in Fig. 13(c). Other isotherms are 
calculated in the same way and can be considered as a second 
approximation to the true temperature field. By a repetition of 
the procedure this temperature field can be approximated to any 
degree provided the procedure converges. 

The dimensionless Nusselt number can now be obtained from 
the equation 


Pr D,*+* 


Nu = — 


(25) 


which is easily derived from Equations (9), (19), (23), and from 
the definition of a Nusselt number. The value @,*+ is the di- 
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mensionless parameter describing the fluid bulk temperature and 
can be obtained from the velocity field and temperature field re- 
sulting from the calculation just described. An equation cor- 
responding to Equation (25) describes the Nusselt number Nu 
when the parameter 6, ** describing the average wall temperature 
replaces the parameter 4)**. 

The actual calculation following this procedure turns out to be 
quite tedious because the location O as well as the shape of the 
isotherms has to be varied in the course of the iteration procedure. 
We were therefore satisfied to replace the curved boundaries of 
AA by straight lines. A variation of O along the duct center line 
leads to a position for which Equation (22) was fulfilled around 
the duct periphery within +20 per cent. A study of the con- 
vergence obtained to that point indicated that the result is not 
accurate enough for a description of the local Nusselt numbers, 
but that it should approximate, on the other hand, the average 
Nusselt number fairly well. The value of the Nusselt number 
obtained in this way is inserted as a full point in Fig. 10. 


DISCUSSION 
L. C. Hoagland* 


The authors should be commended for an excellent piece of ex- 
perimental work and a very significant contribution to the rather 
sparse literature on the subject of heat transfer in noncircular 
ducts. Research on heat transfer to fluids in noncircular ducts is 
fraught with many difficulties not encountered in work with 
circular ducts. Perhaps the most significant difficulty is the 
sensitivity of the average heat-transfer coefficient or Nusselt 
number to changes in the thermal boundary condition applied 
around the duct periphery. 

As stated by the authors, previous analytical work has shown 
that “the Nusselt number averaged over the periphery may dif- 
fer by an order of magnitude depending on whether a constant 
(uniform) wall temperature or a constant (uniform) heat flux 
from the duct wall into the fluid is prescribed peripherally.” 
These two extremes, of course, represent the limits of actual 
boundary conditions that may be achieved in practice. It ap- 
pears that the authors have obtained experimental results for a 
boundary condition approaching the uniform peripheral heat- 
flux case. (It would be interesting to know how the heat flux 
into the fluid varied peripherally in these tests.) 

Previous analytical work reported in the literature has indicated 
that for extremely noncircular shapes such as the 5 to 1 triangle, 
the average Nusselt numbers for uniform peripheral heat flux 
are considerably lower than those for uniform peripheral wall 
temperature. Therefore it is not too surprising that the average 
Nusselt numbers reported in Fig. 10 are considerably lower than 
one would predict from the Dittus-Boelter equation. It is hoped 
that the authors will continue their program to investigate other 
values of the parameter k,b/kD», and, in particular, the uniform 
peripheral wall-temperature case to determine whether the Nus- 
selt number for this case approaches the Dittus-Boelter correla- 
tion. This, of course, is easier said than done. 

The experiments reported here indicated that very long duct 
lengths (measured in terms of number of hydraulic diameters) 
are required to develop the temperature profile along the axis of 
symmetry (across the long duct dimension). The author has re- 
cently been doing some work at M.I.T. on the thermal entrance 
region for laminar slug flow in rectangular ducts. Although it is 
impossible to make a direct comparison between laminar slug 
flow and turbulent flow, some of the results obtained to date may 
help to explain this seemingly strange discovery. Two thermal 
boundary conditions have been investigated, namely, (a) uniform 
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Groetz number (Gz) fe Pr 
L/D, 


short side (axis) long side 


a Uniform heat flux everywhere 


short side long side 


b Uniform wall temperature everywhere 


Fig. 14 Laminar slug flow in 4 to 1 rectangular duct—development of 
temperature profiles across long and short axes 


heat flux to the fluid both peripherally and axially, and (b) uni- 
form duct wall temperature both peripherally and axially. Fig. 
14 shows for a 4:1 aspect ratio duct typical temperature distribu- 
tions along the long and short axes of symmetry at various 
positions from the duct entrance [measured by the conventional 
Graetz number Gz = RePr/(L/D,)|. These curves indicate 
that the temperature profiles along the long and short axes be- 
come fully developed at widely different values of the Graetz 
number (based on the conventional hydraulic diameter) or, in 
other words, at considerably different distances from the entrance 
to the heating section. This effect appears about the same quali- 
tatively for both boundary conditions. 

Of special interest is the manner in which the Graetz number 
required to develop the temperature distribution across the long 
and short axes varies with the duct aspect ratio. Fig. 15 shows 
this relationship using the criterion that the temperature parame- 
ter (see Fig. 14) must everywhere be within 5 per cent of its value 
when fully developed. The Graetz number required for developing 
the short side profile increases by a factor of 3 as the aspect ratio 
goes from one to ten whereas the Graetz number required to de- 
velop the long side profile simultaneously ‘‘decreases’’ by a factor 
of 30. Note that a lower Graetz number means a larger L/D, 
for fixed values of Re and Pr. For comparison, the Graetz num- 
ber required to develop the temperature profile in a circular duct 
is about 15 or 20. 

At this point it is interesting to philosophize a bit. Should the 
thermal entrance length for noncircular ducts be comparable to 
that for circular ducts? Should the hydraulic diameter be a 
meaningful parameter determining the thermal entrance length 
for noncircular ducts? Physically, in order to obtain a fully de- 
veloped temperature profile (say, across the long axis of sym- 
metry), the effect of the sudden change in the thermal boundary 
condition at the entrance of the heating section must propagate 
across the long axis from the wall to the center, ‘a distance a. 
The propagation phenomenon is accomplished physically by heat 
conduction and is inhibited by convection and thePmal storage 
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required to develop temperature profiles across long and short axes of 
symmetry (within 5 per cent) 


ratio - 


processes. Hence the thermal entrance length L required for the 
boundary effect to propagate across the distance a should be 
directly proportional to the fluid velocity u and the thermal 
capacity of the fluid through which the effect must propagate 
(pcyAn° a) and inversely proportional to the thermal conductance 
of the fluid (kA,,/a) where Aj is the area normal to long axis or a. 
Thus 


‘a 

This physical notion suggests that a Graetz number built with the 
duct dimension representative of the temperature profile in ques- 
tion should provide a better thermal entrance length parameter 
than the Graetz number built with the hydraulic diameter. When 
the results shown in Fig. 15 are converted to Graetz numbers 
based on the width (2a) and height (2b) of the duct, respectively, 
the Graetz number required to obtain fully developed temperature 
profiles is indeed a constant (for a given boundary condition) and 
is the same for both long and short sides and independent of aspect 
ratio. 

Although the solutions illustrated in Figs. 14 and 15 indicate 
large entrance lengths required to develop temperature profiles 
across the long duct dimension, another boundary condition of 
practical importance should be expected to have a much smaller 
thermal entrance length. Consider the boundary condition of 
uniform heat input per unit duct length with uniform peripheral 
wall temperature. For this case the fully developed temperature 
distribution across the long side is completely different in charac- 
ter from the “‘parabolic type’’ temperature profiles obtained with 
uniform wall temperature and uniform heat flux everywhere. 
Fig. 16 shows fully developed temperature profiles across the long 
axis of rectangular ducts with uniform heat input per unit length 
(uniform peripheral wall temperature). The temperature dis- 
tribution is nearly parabolic for a square duct, but is markedly 
flattened in the center for the larger aspect ratios.‘ Because here 


* This fact can also be visualized physically by noting that the mem- 


brane or soap film analogy (generally used for torsion of noncircular 
bars) is applicable to this problem. 
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Fig. 16 Laminar slug flow in rectangular ducts—fully developed tem- 
perature profiles across long axis of symmetry 


the ultimate effect of the boundary condition is localized, the 
thermal entrance length required to develop the long side tem- 
perature profile is probably considerably smaller for this boundary 
condition than for those mentioned previously. As yet, no 
analytical solution is available for the thermal entrance length 
under the influence of this boundary condition. It will be interest- 
ing to see if this notion can be verified by a solution. If so, it 
appears that a knowledge of the wall conductivity parameter 
(k,.b/kD,) is essential for determining the thermal entrance 
length when a uniform heat flux is generated in the duct wall. 

Finally, it is interesting to compare the duct length required to 
develop temperature profiles along the long and short sides based 
on the results shown in Fig. 15. Consider a rectangular duct with 
an aspect ratio of 5:1, a Reynolds number of 1000, and a Prandtl 
number of 1.0. For the short side, a Graetz number of 40 is re- 
quired. The corresponding length to hydraulic diameter ratio 
is 25. For the long side, the required Graetz number is 1.5 and 
the necessary length-diameter ratio becomes 667. The writer 
plans to finish this work soon and submit a paper with a complete 
description to the ASME. 


Authors’ Closure 


The preceding discussion of our paper by L. C. Hoagland 
focuses attention on the question: Which length should be used 
as basis for a dimensionless parameter characterizing the duct 
length required for the development of the temperature field? 
Mr. Hoagland was successful in finding such a parameter for 
laminar bulk flow through a rectangular duct the value of which 
is independent of the side ratio of the duct cross section. In a 
similar attempt, it was investigated in our paper whether the 
ratio of duct length L necessary for development of the tempera- 
ture field to duct height would have a value which agrees better 
with the ratio of L to diameter D for a circular tube than the 
ratio L/D,. The experiments showed, however, that the tem- 
perature field in turbulent flow through a triangular duct was 
not developed for a ratio of duct length to duct height equal 
21.5, whereas it has been found to be developed for a ratio L/D 
from 10 to 15 in a circular pipe. This led us to the conclusion 
that the turbulent flow process itself must be different in the 
triangular duct and in the circular tube. 

In answer to a suggestion by Mr. Hoagland, we are presenting 
in Fig. 17 a parameter which is proportional to the local heat 
flux g,, plotted over the distance z from the sharp corner of the 
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Fig. 17 


duct for a specific Reynolds number. It may be observed that 
the heat flux varies considerably over the duct periphery indicat- 
ing that the specific conditions applied in our experiments led to 
a situation approximately halfway between a constant heat 
flux and a constant wall temperature boundary condition. 

Friction and heat transfer in turbulent flow through a rectan- 
gular duct are presently studied in our Heat Transfer Laboratory. 
A report on the friction characteristics is presently being pre- 
pared. 
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Heat Transfer in Flow Through 
Rotating Ducts 


Experimentally measured heat-transfer coefficients and visually observed flow phe- 


nomena for water flowing through an electrically heated rotating pipe with various 


inserts are reported. 


Convection heat transfer in flow through a rotating pipe or a rotating annulus was 
found to be an exceedingly complex mechanism which can not be treated adequately in 
terms of conventional parameters, but requires a detailed understanding of wave and 
vortex phenomena which were visually observed and photographed. Heat-transfer co- 
efficients were measured in completely as well as partially filled flow channels. Three 
distinct flow regimes were observed in all cases; transition between regimes was found 
to depend on the rotational speed, but not on the axial flow velocity. Stationary vortexes 
and a counterrotational fluid motion were observed. 


OUGH CONVECTIVE HEAT TRANSFER in rotating 
systems is of importance in the thermal design of many types of 
industrial equipment, relatively little information on heat trans- 
fer in such systems is available in the literature, especially with 
fluids other than air. 

This study was originally undertaken under the sponsorship of 
the National Printing Ink Research Institute in order to obtain 
information on the factors which control the cooling of printing 
ink roll mills in which an excessive temperature rise in the 
products is undesirable. The results of the initial phases of the 
research program raised, however, more questions than they 
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answered. The scope of the study was therefore extended under 
a grant from the National Science Foundation to include a basic 
investigation of convection heat transfer to, and flow phenomena 
in, water flowing through a heated rotating pipe and annulus. 
Unfortunately, even the extended scope of the research program 
leaves many questions still unanswered and also raises new ones. 
It does, however, give an insight into some of the complexities of 
the transfer mechanism and can therefore serve as the basis for 
further study. 

Heat transfer and flow phenomena in systems in which the 
fluid motion is influenced by body forces and Coriolis forces have 
until recently been investigated only by meteorologists and astro- 
physicists. In recent years, however, practical problems in aero- 
nautics, atomic power, gas turbine cooling, printing and paper 
roll design, and chemical engineering have drawn the attention of 
engineers to problems associated with convection in fluids subject 
to various force fields. 9 

This paper presents experimental results on heat transfer and 
flow phenomena in a fluid flowing through a horizontal, elec- 
trically heated pipe and annulus which are rotating and by virtue 
of their rotation cause the fluid to rotate. The rotation vector in 
this system is horizontal, and the temperature gradient is radially 
orientated. Table 1 summarizes the information available in the 
literature on some other rotating diabatic systems of interest. 
Papers dealing with heat transfer from bodies rotating in an in- 
finite environment, such as cylinders and disks, are not included 


= heat capacity, Btu/lb deg F 


ID of test section, ft 

equivalent diameter, ID of outer 
tube minus OD of inner tube, ft 

separation distance between sur- 
faces rotating in opposite 
directions; D,/48, in. 

internal energy of water layer, 
Btu 

gravitational constant, ft/hr* 

convective heat-transfer coef- 
ficient, Btu/hr ft? deg F 

thermal conductivity, Btu/hr 
ft? deg F/ft 

depth of fluid in free convection 
studies, ft 
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rotational speed, rev/sec 

rate of heat transfer, Btu/hr 

radius of inner tube, in. 

inlet water temperature, deg F 

outlet water temperature, deg F 

average wall surface tempera- 
ture, deg F 

average bulk water tempera- 
ture, deg F 

average water velocity, ft/sec 

mass flow rate of water, lb/hr 

radial distance from inner sur- 
face of outer tube, ft 

Prandtl number, c,u/k 

flow-rate Reynolds number, 
DV/v 


Nre, = rotational Reynolds number, 
D*w/v 
= flow-rate Reynolds number based 
on equivalent diameter, D,V/v 
Nusselt number, hD/k 


d 
Ro d2rRyn/v 


coefficient of volumetric expan- 
sion, deg 
thickness of water film on upper 
inside surface of outer tube 
when tube is partly full, ft 
= viscosity, lb/ft hr 
= kinematic viscosity, ft?/sec 
= angular velocity, sec™ 
= density of water, lb/ft* 


Taylor number, 
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Table 1 Summary of types of flow investigated in rotating diabatic systems 


Reference no. 1 2 3 4 5 6 7 8 9 This paper 
Temperature gradi- 
ent R V V V R R R R R R 
Rotation vector V V V V ¥ V H H H H 
Geometry A HL HL M4&HL A A A A A P, A 4 
Type = E E E,8 T E 
Legend: 
H—horizontal P—pipe 
V—vertical HL—horizontal layer PIPING ARRANGEMENT 9 
R—radial E—experimental 
A—annulus T—theoretical " HALF FULL SYSTEM 
S—summary 


in Table 1 because in these systems no phenomena directly per- 
tinent to this study have been reported. Summaries of the in- 


formation available on heat transfer from bodies rotating in an 
infinite environment are presented in Refs. [8 and 12],* and on the 
mechanies of the flow in rotating systems without heat transfer in 
Refs. [7, 8, and 11]. 


Equipment and Procedure 


Heat Transfer. A schematic flow diagram of the equipment is 


shown in Fig. 1. Water was supplied to the constant head tank “ 
through valve A. The water flow rate was regulated by valve E : 
and measured by a rotameter, R, before entering the test sec- 


tion. 
Two piping arrangements were available at the test section PIPING ARRANGEMENT 
outlet. In the first arrangement, shown in Fig. 1(a), the outlet FILLED SYSTEM 


duct was horizontal. With this arrangement, corresponding to 
that employed in many ink and paint dispersion mills, the test 
section was only half filled with water, and the space above the 
water was filled with air. With the alternate arrangement shown 


in Fig. 1(6), the test section was completely filled with water. 
The air was first removed through a hole in the top of the pipe, 
and the hole was then plugged. A standpipe with a vent at the 
outlet insured full flowing conditions in the test section during the , 

«MERCURY JUNCTION 


The heat-transfer test section is shown in detail in Fig. 2. It 
consisted of a horizontal, 3-in-ID steel tube, 25'/, in. long with 
'/~in. thick walls. When assembled, the internal length was 24 


in. Heat was supplied to the section by radiation and convection Fig. 1 Schematic sketch and flow diagram of experimental equipment 5 
(a, arrangement for partially filled ducts; b, arrangement for full flowing | 
5 Numbers in brackets designate References at end of paper. duct) 
3° ID STEEL ROTATING 4 
ROLL 
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Fig. 2 Heat-transfer test section—construction details 
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from 32 stationary, electrically heated, No. 15 Nichrome, re- 
sistance wires, running parallel with the axis, near the outside 
surface of the tube. The test section was placed in an insulated 
cabinet to which the heating wires were attached. Electrical 
power input was regulated by a rheostat and measured with a 
wattmeter for establishment of operating conditions. 

Inside surface temperatures of the tube were measured by 
three thermocouples embedded in the wall in accordance with the 
Patton-Feagan method [20]. These thermocouples were con- 


24-0” 


nected to a potentiometer through mercury pool junctions as 
shown in Fig. 4. Inlet and outlet water temperatures were 
measured with thermocouples immersec in the water stream. 
The wall and water thermocouple locations are shown in Fig. 3. 
Readings from thermocouples 1, 2, 3, 6, and 9 were used in the 
heat-transfer calculations. 

Tests were conducted with two types of inserts in addition to 
the plain tube studies. These inserts are shown in Fig. 5. The 
first type, the paddles-insert, consisted of four baffles made of 
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Water Thermocouples 
Fig. 3 Location of thermocouples 
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Fig. 4 Construction of mercury junction for connecting rotating thermocouple leads to potentiometer 
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X 9/4 24-in. steel strips. The second type, the annulus- 
insert, consisted of a 2*/s-in-OD steel tube, 22'/2 in. long, sealed 
at both ends and supported by '/s X */4 X 3-in. crosspieces as 
shown in Fig. 5. Both types of inserts rotated along with the 
main tube. 

Independent variables were the water flow rate, rotational 
speed, and power input. Dependent variables were the tem- 
peratures of the duct wall and the temperature rise of the cooling 
water. Measurements of all of the variables were taken at each 
experimental condition for at least 15 minutes after equilibrium 
conditions had been established. 

The rate of heat transfer was determined by the equation 


we, (te come ti), (1) 


and the average convective heat-transfer coefficient h was calcu- 
lated from 


h = q/A(t, — t,); (2) 


where w is the mass flow rate of the cooling water in lb/hr, c, is 
the heat capacity of the water in Btu/lb deg F, 4, is the inlet tem- 
perature of the water in F, t; is the outlet temperature of the water 
in F, A is the total inside heat-transfer surface of the duct in sq 
ft, t, is the arithmetic average temperature of the three wall 
thermocouples, and ¢, is the mean between the inlet and outlet 
water temperatures. 

The convective heat-transfer coefficients measured in this study 
are estimated to be reliable within 28 per cent at the lowest heat 
flux, but within about 5 per cent in the higher heat flux range. 

Flow Visualization. The equipment used for the flow visualiza- 
tion study was essentially the same as that used for the heat- 
transfer tests except that the heater and insulation were removed, 
and a transparent tube was substituted for the steel tube. The 
transparent duct was made from 3-in-ID plexiglas tubing, 22'/, 
in. long with a wall thickness of */;.in. In order to study the 
annulus flow, a 19*/,-in. length of 2'/:-in-OD plexiglas tube, with 
plugged ends, was fastened in the center of the rotating duct. 
Crushed wheat particles were added to the system in order to aid 
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Fig. 5 Flow cross sections of test sections and sketches of inserts 


visual observation of the water flow patterns. Visual studies 
were made with and without axial flow, at rotational speeds in 
the range from 10 to 500 rpm using water at an average tem- 
perature of 68 F. The movements of the wheat particles were 
observed both visually and with photographs taken with a 35-mm 
camera, tri-X film, f-11 lens opening, and '/,-sec exposure. 


Experimental results for the geometrical configurations studies 
are shown in plots of Nusselt number, hD/k, versus rotational 
Reynolds number (Re, = D*w/v), in Figs. 6 through 10. The 
characteristic length dimension D is the inside diameter of the 
test section, 0.25 ft. The parameter employed in presenting dif- 
ferent flow-rate curves in each of these figures is the water flow rate 
in lb/hr rather than the flow-rate Reynolds number, because the 
latter varied with temperature while the flow rate was constant in 
each series of tests. The relationship between flow rates and aver- 
ages of the experimental flow-rate Reynolds numbers is shown in 
Table 2. The lack of strict proportionality between flow rate and 
Reynolds number is due to differences in average water tempera- 
ture in the different flow-rate ranges. 

The data of Fig. 6, for the partly full, simple rotating tube, 
yield straight lines when plotted as Nusselt number versus 
logarithm of the flow Reynolds number (Re, = DV/yv), at con- 
stant rotational Reynolds number. Lines for the various rota- 
tional speeds have a common intercept at a Reynolds number of 
24 on the log Reynolds number axis. It is therefore possible to 


Table 2 Correspondence between flow rate and experimental average 
flow-rate Reynolds number 


Flow rate, Reynolds number, 
lb/hr DV/» 


29 60 
60 120 
168 300 
317 550 
465 800 
530 950 
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Fig. 6 Nusselt number versus rotational Reynolds number for partially filled plain tube at various flow rates 
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Fig. 7 Nusselt number versus rotational Reynolds number for full flowing plain tube at various flow rates 
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Fig. 9 Nusselt number versus rotational Reynolds number for full flowing annulus at various flow rates 
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Fig. 10 Nusselt number versus rotational Reynolds number for particlly filled pipe with paddle insert at various flow rates 
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Fig. 11 Correlation of Nusselt number with flow rate and rotational Reynolds numbers for partially filled plain tube 
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Fig. 12 Calibration of experiment— pori of heat-transfer coefficient for flow through the annulus without rotation 


with data of Ref. [16] 


correlate the data of Fig. 6 in a single plot of log [Nu/(Re,/24)] 
versus log Re,. The latter relationship is shown in Fig. 11. The 
relationship between flow rate and heat t:ansfer was more com- 
plicated in the other geometrical systems, and the data presently 
available do not warrant more extensive correlation than shown 
in Figs. 7 to 10. 

A check of the experimental equipment was provided by ob- 
taining heat-transfer data without rotation in the completely 
filled annulus system. The results are compared with those of 
Chen, Hawkins, and Solberg [16] in Fig. 12. This comparison 
shows that the slopes of the two curves are similar, and that the 
experimental data agree within 8 per cent with those of Ref. [16]. 
The characteristic diameter used in the Reynolds number of Fig. 
12 is the equivalent diameter D, which is equal to the ID of the 
outer tube minus the OD of the inner tube. For the data of this 
paper, D was 0.0313 ft. 


Discussion 


An inspection of the experimental results shown in Figs. 6 
through 11 reveals the following significant points: 


1 There exist two discontinuities in the slopes of the curves of 
Nusselt number versus Reynolds number for all of the test sec- 
tions. The first discontinuity occurs at a rotational speed of 
about 30 rpm, and the second at about 150 rpm. 

2 The locations of these discontinuities are essentially in- 
dependent of the net mass flow rate except at the highest flow 
rates in the completely filled annulus shown in Fig. 9. 

3 At any given rotational speed and net flow rate the Nusselt 
number for a partially filled duct is equal to, and in some cases 
even slightly larger than, the Nusselt number for the same duct 
completely filled with water. 


In the following discussion, an attempt is made to present 
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plausible, but still tentative, explanations for these observations. 
The authors wish to emphasize the tentative nature of this dis- 
cussion, but feel that its presentation is justified by its agree- 
ment with the experimental observations. Additional tests will 
be required to validate these arguments under conditions different 
from those investigated here. 

At first sight, the conditions in a rotating annulus with fluid in 
the annular space heated from the outer surface resemble the 
conditions prevailing in free convection between two horizontal 
flat surfaces with heating from below [13]. In the rotating an- 
nulus the centrifugal body force is directed perpendicular and 
toward the heating surface, while in the free convection system 
the gravitational body force acts in a similar manner. For the 
gravitational free convection system it is well known that mo- 
tion begins when the vertical temperature gradient becomes suf- 
ficiently large, a condition corresponding approximately to a 


32 
Rayleigh number, SoeeT of 1700. Chandra [17] showed 


that the motion can either be “cellular” or “columnar,” the form 
of motion being dependent on the depth of fluid. When a steady 
flow is superimposed parallel to the heating surface, the fluid was 
found to align itself into a series of vortexes with axes transverse 
to the steady flow direction as soon as the steady flow velocity 
becomes sufficiently large to destroy the columnar flow pattern. 
For even larger steady velocities (the exact transition value is 
not known) the flow pattern changes into a series of vortexes with 
axes parallel with the flow direction, i.e., a series of longitudinal 
rolls. 

The situation in the rotating annulus or pipe is, however, con- 
siderably more complicated, because within the range of rotational 
speeds investigated the gravitational torce g is of the same order 
of magnitude as the centrifugal force. For example, at a rota- 
tional speed of about 150 rpm, the centrifugal force is just equal 
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Fig. 13 Heat flux versus temperature potential for partially filled rotating plain tube 


to one g. Since the component of the gravitational force per- 
pendicular to the heating surface varies with angular position, 
the analogy between the rotating system and horizontal free 
convection is only approached at high rotational speeds when the 
gravitational force becomes negligible compared with the centrif- 
ugal force. It will be shown in the following discussion that in 
the lower range of rotational speeds investigated here, the flow 
conditions appear to be somewhat analogous to those found by 
G. I. Taylor [18] and others in an annular space between rotating 
co-axial cylinders. H. Jeffreys [10], who investigated mathe- 
matically the stability of a layer of incompressible fluid with a 
decreasing temperature in the vertically upward direction with 
a steady crossflow superimposed parallel to the surface, found that 
the instability in this system is analogous to that observed by G. 
I. Taylor [18] for flow without heat transfer between rotating 
co-axial cylinders. It is therefore reasonable to expect that a 
transition from conditions where the gravitational force is sig- 
nificant to the conditions under which centrifugal effects predomi- 
nate would be subject to similar flow regimes and stability phe- 
nomena as have been found in each system separately. In order 
to determine experimentally whether or not pure free convec- 
tion phenomena analogous to those between two horizontal sur- 
faces are of importance for the rotating duct within the range 
of rotational speeds and temperature gradients investigated, the 
heat-transfer coefficient in the partially filled plain pipe was 
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measured at various wall-to-fluid temperature differences while 
holding all other variables constant. As can be seen from Fig. 
13, the heat-transfer coefficient is essentially independent of the 
temperature difference, since the rate of heat flow is directly pro- 
portional to the mean temperature difference. From these tests 
it was concluded that those flow phenomena which depend on 
the rotational speed and the mass flow rate predominate over 
those induced by pure free convection. This may, however, not 
be true at rotational speeds or temperature gradients considera- 
bly different from those used in this study. 

An inspection of Figs. 6 through 10 reveals that the variation of 
the Nusselt number with rotational speed changes abruptly at 
rotational Reynolds numbers of about 2 X 10‘ and 7 X 10‘, cor- 
responding to rotational speeds of about 30 and 150 rpm, respec- 
tively. In order to shed some light on the reason for the change 
in the functional dependence of Nyu on Nre,, & transparent test 
section having a geometrical configuration similar to that of 
the section used for the heat-transfer study was built of plexi- 
gias, and the flow was observed visually. It was observed that 
the general features of the flow pattern did not vary qualita- 
tively with flow rate, but depended primarily on the rotational 
speed. Since it was considerably easier to photograph the flow 
pattern without any net through flow, the pictures shown here 
were taken without net flow through the system. 

The series of photographs in Fig. 14 show the flow pattern in 
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i5 RPM FRONT 


33 RPM FRONT 


33 RPM REAR 


65 RPM FRONT 


500 RPM 


the plain pipe, half filled with water, at various rotational speeds. 
The dark areas in the photographs are the crushed wheat particles 
mentioned previously. 

At 15 rpm the particles are almost uniformly distributed along 
the tube. At 33 rpm the particle motion indicates the existence 
of ring-shaped vortexes spaced regularly along the pipe. Two pic- 
tures are shown for this speed. One, identified as ‘‘Rear,” was 
taken from 15 degrees below the horizontal, looking along the 
direction of rotation at the cylinder top. The picture marked 

- “Front” was taken at an angle of 45 degrees from the horizontal, 
looking down in the direction opposite to that of rotation. The 
wave length of these vortexes, i.e., the spacing between their 
axes, is about one half the pipe diameter. At 65 rpm the wave 
length has increased to about one pipe diameter, but the main 
vortex is beginning to break up into smaller disturbances. This 

break up continues at 110 and 138 rpm, and at about 160 rpm a 
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Fig. 14 Photographs of flow in partially filled plain pipe without net through flow 


110 RPM FRONT 


RPM 


FRONT 


ig 


198 RPM FRONT 


FRONT 


sheet of water of substantial thickness is carried along by the tube 
wall. It should be noted that at 150 rpm the centrifugal force 
just balances the gravitational force. At 200 rpm the flow pat- 
tern is irregular, but at 300 rpm a more regular vortex motion be- 
comes apparent again. The wave length is now only about one 
fifth of the tube diameter. At 360 rpm a new flow pattern ap- 
peared. At this speed there were 6 distinct bands of entrapped air 
observed along the inner circumference of the tube. At 450 rpm 
the number of bands increased to 9; while at 500 rpm, the last 
picture shown, 11 bands were counted. Seven of these bands 
can be seen in the portion of roll photographed. 

At speeds below 150 rpm the flow pattern in the bulk of the 
fluid, as seen from an end view, i.e., in a cross section perpendicu- 
lar to the axis of rotation, is sketched in Fig. 15. The fluid in 
the central portion of the pipe rotates in a direction opposite to 
the pipe. This somewhat unexpected flow pattern may be neces- 
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Fig. 15 Sketches of flow pattern in partially filled pipe and annulus in a cross section 


perpendicular to the axis of rotation 


sary to conserve the moment of momentum of the air and water 
in the test section. A similar reversal in the flow was also ob- 
served in the water in the rotating annulus [see Fig. 15()]; it 
may serve to explain the discontinuity at 30 rpm. 

In 1923, G. I. Taylor [18] investigated the stability of an in- 
compressible viscous fluid contained in a narrow annulus between 
rotating concentric cylinders of infinite length for the case of zero 
axial flow. In 1937, 8. Goldstein [19] extended the analysis to 
the case of a viscous flow through an annulus formed between a 
rotating inner cylinder and a stationary outer cylinder. These 
papers, as well as additional mathematical treatments and ex- 
perimental investigations are summarized in Refs. [7 and 8). 
Although Taylor’s analysis leads to the conclusion that no in- 
stabilities should occur when the inner and outer cylinder rotate 
at the same speed, and Goldstein's analysis is inconclusive on that 
point, one can expect instabilities to arise in the partially filled 
annulus and pipe as shown by the following considerations. 
According to the observed cross-sectional flow pattern, the water 
in the center of the annulus rotates in the direction opposite to 
that of the walls at low speeds. The flow system, therefore, 
roughly resembles two concentric annuli in which the inner and 
outer surfaces rotate in opposite directions. Meskyn’s 
[15] predicts for this case a critical Taylor number of about 85 
without net flow. For the idealized system above, the critical 
Taylor number at 30 rpm, the speed at which the first discon- 
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analysis 


tinuity in the slope of the curve of Nyu versus Npre occurs, is, 


Noe = V/d/R, d2xRon/v (3) 


- vm 0.094 X 6.28 X 1.5 X 0.5/(144 X 1X 107-5) = 82 
15 


The agreement may be fortuitous, but the same approach has 
been used successfully by Brewster and Nissan [14] on a similar 
problem with satisfactory results. 

A similar argument can be used to explain the origin of the 
vortex pattern in the partially filled pipe where the flow in the 
outer annulus resembles the flow between rotating co-axial 
cylinders. 

The preceding arguments admittedly do not apply to the full- 
flowing pipe and annulus. In these two cases it was very difficult 
to observe the flow with the available equipment, and although 
vortexes were observed and the curves of Nyw versus Nre, show 
similar discontinuities in slope, neither a reliable description of the 
flow nor a satisfactory explanation of the cause of the instabilities 
can be given at this time. 

Having shown that the discontinuities in slope of the curves of 
Nwu versus Ne, coincide, at least in the partially filled systems, 
with changes in the flow pattern, we now turn our attention to a 
comparison of the Nusselt number with rotation to the Nusselt 
number without rotation. For an annulus filled with air, Bjork- 
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Fig. 16 Ratio of Nusselt number with rotation to Nusselt number without rotation versus rotational 


Reynolds number 


lund [8] found that, at Taylor numbers less than critical, the 
Nusselt number in the annular space between rotating cylinders 
equals the Nusselt number for pure conduction. A similar com- 
parison for the rotating annulus is shown in Fig. 16, where the 
ratio of the Nusselt number with rotation to the Nusselt number 
without rotation at the same net flow rate is plotted as a function 
of the rotational Reynolds number. An inspection of this curve 
reveals that this ratio is approximately unity at rotational 
Reynolds numbers below those corresponding to the critical value 
of Taylor numbers for the idealized system. In other words, rota- 
tion does not increase the Nusselt number until vortexes begin to 
appear. 

The situation is qualitatively similar for the plain pipe, but in 
this system the Nusselt number was found to increase continu- 
ously with rotational speed, even at low rotational speeds. 

The tests with the paddles insert were conducted chiefly to de- 
termine whether or not this type of insert would improve the heat 
transfer at a given flow rate. Contrary to expectations the pad- 
dies did not improve the heat transfer, but actually reduced the 
Nusselt number at rotational speeds above 30 rpm. A reasonable 
explanation for this observation is that the motion induced by 
the vortexes is more effective in transferring heat by convection 
than the motion induced by the paddles which undoubtedly inter- 
fered with the formation of the ring-shaped vortex patterns. 

The third significant observation was that the heat-transfer co- 
efficients for a partially filled duct at a given rotational speed 
and water mass flow rate is essentially equal to that for a com- 
pletely filled duct at the same speed and flow rate. At first sight 
this seems unreasonable because in a partially filled duct only one 
half of the inside area is in contact with the main body of water, 
and consequently the convection heat-transfer surface area is 
much smaller. However, it was observed that under all condi- 
tions a film of water, approximately 0.025 in. thick, adhered to 
that portion of the inner pipe surface which was above the water 
level in partially filled flow. This thin layer of water will be 
heated by conduction from the pipe, but is effectively insulated 
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on the other side which is in contact with the air in the upper 
half of the pipe. Little or no mass transfer is expected to take 
place because the air is saturated in the steady state. 

If we now postulate that this layer of water which adheres to 
the upper pipe surface mixes effectively with the main stream 
after the rotational motion submerges it again in the bulk of the 
water, one can ask if this periodic heat flow mechanism can 
possibly be as effective as the transfer mechanism in a completely 
filled duct. 

Assume for the sake of this argument that the convection heat- 
transfer coefficients over that portion of the duct surface which is 
in contact with the main stream are equal for partially filled and 
full flowing conditions. Then the rate of heat transfer through 
the lower half of the pipe surface per unit length is 


q = h(wD/2)(t, — t,) (4) 


in a half full duct. In the upper half of the duct the internal 
energy of a water layer of thickness 6 increases during one revolu- 
tion by the amount 


where y is the radial distance from the surface 


t:(y) is the temperature profile of the film as it emerges from the 
main body of water 
t,(y) is the temperature distribution of the film as it submerges 
again 
Assuming that the temperature profile of the emerging film is 


linear and that the entire layer is heated up to t, by conduction, 
the internal energy change per revolution becomes 


where ¢; is the water temperature at the air interface. 
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If the water film mixes with the bulk upon being re-immersed, a 
heat balance per unit time yields 


q = 90wipc,D(t, — ts) + h(wD/2)(t, — ty) (7) 
for the partially filled duct, while 
q = hrD(t, — t,) (8) 


for the full flowing case. Combining Equations (7) and (8) and 
solving for 6 yields 


whit, — t,) 
1800.pe,(t, — ta) 


Assuming (t, — t;) = 2(¢, — ts), Equation (9) yields 


th 
900wpc, 


6 (10) 


Equation (10) can be recast into the more convenient dimension- 
less form, 

3) = 4r 

D Re, 
Substituting the experimental data from Figs. 6 and 9 into Equa- 
tion (11) shows that the film thickness 6 required to have equal 
rates of heat transfer in a partially and a completely filled duct 
is about 0.08 in. at Nre, of 4 & 10* and decreases to about 0.008 
in. at Nre, of 2 XK 10°. Although experimental observations indi- 
cate that the actual film thickness slightly increases with increas- 
ing speed, the preceding analysis is still valid because with in- 
creasing speed less time per revolution is available for heating the 
film. The change in internal energy, however, can still be more 
than sufficient to offset the reduction in surface area available for 
convection. This mechanism may have numerous practical 
applications in the thermal design of milling and calendering rolls. 


(11) 
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Comparison of Turbulent Heat-Transfer Results for Uniform 
Wall Heat Flux and Uniform Wall Temperature 


R. SIEGEL’ and E. M. SPARROW: 
Introduction 


Tae purpose of this note is to examine in a more precise way 
how the Nusselt numbers for turbulent heat transfer in both the 
fully developed and thermal entrance regions of a circular tube 
are affected by two different wall boundary conditions. The 
comparisons are made for: (a) Uniform wall temperature 
(UWT); and (5) uniform wall heat flux (UHF). Several papers 
which have been concerned with the turbulent thermal entrance 
region problem are given in references {1 to 4].2 Although these 
analyses have all utilized an eigenvalue formulation for the ther- 
mal entrance region (in contrast to reference [6] which used a 
boundary layer approach), there were differences in the choices of 
eddy diffusivity expressions, velocity distributions, and methods 
for carrying out the numerical solutions. These differences were 
also found in the fully developed analyses. Hence when making 
a comparison of the analytical results for uniform wall tempera- 
ture and uniform wall heat flux, it was not known if differences in 
the Nusselt numbers could be wholly attributed to the difference 
in wall boundary conditions, since all the analytical results were 
not obtained in a consistent way. To have results which could be 
directly compared, computations were carried out for the uniform 
wall temperature case, reference [4c], using the same eddy dif- 
fusivity, velocity distribution, and digital computer program em- 
ployed for uniform wall heat flux in references [4a and b]. In 
addition, the previous work was extended to a lower Reynolds 
number range so that comparisons could be made over a wide 
range of both Reynolds and Prandtl numbers. 

The analysis of heat transfer in the turbulent thermal en- 
trance and fully developed regions has already been thoroughly 
treated in references [4] and hence need not be presented here in 
detail. The calculations were carried out under the assumption 
of equal eddy diffusivities for heat and momentum. For the re- 
gion near the tube wall, the diffusivity was evaluated from Deiss- 
ler’s formula [5], while for the region away from the wall, the 
diffusivity was found by differentiating the logarithmic velocity 
expression and using a linear variation of shear stress. For con- 
venience, the Appendix gives a few of the final analytical expres- 
sions, and Table 2 provides numerical data for the cases which 
have not been previously tabulated in references [4]. 

In the next section, we then proceed directly to an examination 
of the Nusselt number results. 


Nusselt Number Comparisons 


Pr = 0.7. The results for Pr = 0.7 and four different Reynolds 
numbers are given in Table l(a) as a function of the length- 
diameter ratio (r/D) from the tube entrance. The per cent dif- 
ference is defined as (Nuvur — Nucwr)/Nuuwr. The largest 


1 NASA, Lewis Research Center, Cleveland, Ohio. Assoc. Mem. 
ASME. E. M. Sparrow is now Professor of Mechanical Engineering, 
University of Minnesota, Minneapolis, Minn. 

? Numbers in brackets designate References at end of paper. 

Contributed by the Heat Transfer Division of Tae AmeEerRIcAN 
Society or MecHANIcAL ENGINEERS. 


Manuscript received at ASME Headquarters, January 28, 1960. 


152 / may 1960 


After approval such contributions will be published as soon as possible, normally in the next issue of the Journal. 


differences are found in the low Reynolds number range, and for 
each case the differences become smaller at larger distances from 
the tube entrance. For the ranges of variables considered, the 


effect of the two different wall boundary conditions is never’ 


larger than 10 per cent. 

Pr = 10 and Pr = 100. Results for the higher Prandtl numbers 
are given in Tables 1(b) and (c). For this range there is essentially 
no influence of the two different boundary conditions. Since the 
per cent difference does not become significant even for small 
z/D, the details in the thermal entrance region are only given as 
an example for one Reynolds number. Some of the Nusselt num- 
ber results for uniform wall temperature are very slightly in ex- 
cess of those for uniform wall heat flux. This is a somewhat sur- 
prising finding, and it may be inferred that the differences of a 
few tenths of a per cent are probably due to accuracy limitations 
of the calculations. 

Conclusions. From these tabulations we may then conclude 
that for turbulent flow the heat-transfer mechanism in the ther- 
mal entrance and fully developed regions is quite insensitive to 
the two wall boundary conditions which were examined, at least 
for the range Pr > 0.7. This conclusion is in qualitative agree- 
ment with the findings of reference [6] for Pr = 0.73 (thermal 
entrance region calculations using a boundary layer model) and 
of reference [7] for Pr = 1.0 (fully developed region). The effect 
of other wall boundary conditions can be examined by the super- 
position techniques described in references [2 and 4b]. 
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APPENDIX 
Summary of Analytical Results 
Uniform Wall Temperature. The temperature distribution within 
the fluid is given by 
T-T, 


em 


n=1 
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Table 1(a); Pr 0.7 


os 


Table 1(b); Pr 10 


Table 2(a) Uniform wall temperature 
(1) Re = 10,000; Pr = 0.7 


(4@,/ax*) 


-0.036269 
-022654 

- .018059 

.016892 

.016687 

-017637 

- .018552 


50,000; Pr = 


-0 030622 
-017838 
- .013361 
.011209 
- 0098511 
-0090539 
- .0085049 


24167.0 


Table 1(c); Pr 100 


= 100,000 


Re = 100,000 


Table 2(b) Uniform wall heat flux 


diff. 


Re = 10,000; Pr = 0.7 


Nuyyp aiff. 


732. 0.2 
710. 
697. .07 
688. -03 
685. .O1 
683. 


n Ay 


-21646 
-11776 
-090430 
.073326 
-059857 


1556 
1545 
1539 
1534 
1532 
1532 


Fully developed Nusselt numbers 


Fully developed Nusselt 


~ .044084 
- .052717 


Re Nuyyr % aiff. Re 


0.3 


50,000 


500 , 000 


500 , 000 


For Re = 10,000; rf = 324 


where 7’, is the wall temperature, 7’) the entering fluid tempera- 
ture, C,, are constants, 6, and ®, the eigenvalues and eigenfunc- 
tions, z the distance from the tube entrance, and r the radial co- 
ordinate. The Reynolds number is defined by a@D/v, where @ is 
the mean fluid velocity. 

The Nusselt number is given by 


p> (=) ¢—48n%z/DRe 
=1 ro 


Nuvuwr = Pr — 
(4 e 48n*z/DRe 
2 
=i 


Uniform wall heat flux. The fluid temperature distribution is 
given by 
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48n*x/DRe 


n=1 


where q is the wall heat flux per unit area, 7> the tube radius, and 
G(r) is the fully developed temperature distribution. The Nusselt 
number is given by 


1 


Nuvar = 


2 
G(ro) 


n=1 


where A, = C,®,(ro)/G(ro). The ®,, C,, and 8, for uniform heat 
flux are different from those for uniform wall temperature. 

For a detailed description of the analytical method, the reader 
is referred to references [4]. 
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